我首先进行了回归,然后根据系数和置信区间手动计算了三个拟合值。数据框如下所示。 yld
是使用系数计算的拟合值,yld.l
和yld.u
是使用上下置信区间计算的拟合值。
> head(m.m.m)
harv frez_j dd_j cupc_j sm7_j fitted_j max_spring_j sp_spring_j pc_spring_j lt qt yld
1 31793.35 320.0668 2437.82 13.0822 50.63953 46.94146 23.34035 0.0844567 3.260585 18.09834 409.1669 48.78323
2 31793.35 320.0668 2407.26 13.0822 50.63953 46.94146 23.34035 0.0844567 3.260585 18.09834 409.1669 48.73291
3 31793.35 320.0668 2638.01 13.0822 50.63953 46.94146 23.34035 0.0844567 3.260585 18.09834 409.1669 48.90514
4 31793.35 320.0668 2838.45 13.0822 50.63953 46.94146 23.34035 0.0844567 3.260585 18.09834 409.1669 48.72009
5 31793.35 320.0668 2944.05 13.0822 50.63953 46.94146 23.34035 0.0844567 3.260585 18.09834 409.1669 48.52131
6 31793.35 320.0668 2765.93 13.0822 50.63953 46.94146 23.34035 0.0844567 3.260585 18.09834 409.1669 48.81819
yld.l yld.u
1 7.006441e-09 339659325071
2 7.272086e-09 326577117107
3 5.525450e-09 432853975445
4 4.400996e-09 539343282695
5 3.918681e-09 600793423028
6 4.773276e-09 499283037958
我想画三条线,其中y为yld
,yld.l
,yld.u
,x为dd_j
。我可以用下面的代码成功地绘制一条线,却不知道如何绘制全部三条线。
ggline(data=m.m.m,x='dd_j',y='yld',color = 'red',numeric.x.axis = T)
有什么想法吗?
顺便说一句,我认为在ggplot和ggpubr中必须将宽格式表转换为长格式以画多行是非常不便的。