我有一个看起来像这样的数据集:
> averages
compound control.0 control.30 surgery.0 surgery.30
1. A 3.609958 3.578200086 3.556325 3.669107598
2. B 4.984090 4.798330495 4.965342 4.812247664
我只想为化合物 A 绘制两条线的图表 - 一条连接(0,控制.0)到(30,控制.30),另一条绘制(0,手术.0)到(30,手术.30)。我也有 200 种化合物,因此理想情况下,我希望能够让程序进入列表并为每个化合物吐出一个图表,而无需我手动进入并更改行号。我将如何做到这一点?
对于每个图形的两个线段,请使用以下命令:
with(subset(averages, compound=="A"), plot(c(0,30,NA,0,30),c(control.0,control.30,NA,surgery.0,surgery.30), type="l"))
然后更改子集条件以创建循环。