我有一个如下所示的数据集:
> 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,control.0)到(30,control.30)和一条绘制(0,surgery.0)到(30, surgery.30)。我还有200种化合物,所以我希望能够将程序放在列表中并为每个化合物吐出一个图表,而无需我手动进入并更改行号。我该怎么做呢?
答案 0 :(得分:0)
对于每个图表的两个线段,请使用:
with(subset(averages, compound=="A"), plot(c(0,30,NA,0,30),c(control.0,control.30,NA,surgery.0,surgery.30), type="l"))
然后更改子集条件以创建循环。