在 R 中的 1 个图中绘制 6 个系列答案

【问题标题】：Plotting 6 series in 1 plot in R在 R 中的 1 个图中绘制 6 个系列
【发布时间】：2013-10-16 14:10:25
【问题描述】：

我有以下脚本：

myfunctionSD2 <- function(mydata) { return(sd(mydata,na.rm=TRUE))}
SDmeas <- tapply(datIn$Measurement,list(as.factor(datIn$Measurement.location),as.factor(datIn$Tube.number)),myfunctionSD2)
SDMeasurement <- SDmeas[order(as.numeric(rownames(SDmeas))), ]
dput(SDMeasurement)
structure(c(0.6555264303958, 0.605687762634627, NA, 0.683435292388886, 
NA, 0.615207645374612, NA, 0.739018717912613, NA, 0.79341715750565, 
NA, 0.769473075819677, NA, NA, NA, NA, NA, 0.6555264303958, 0.605687762634627, 
0.576875822240553, NA, 0.502849025908516, NA, 0.516028792109233, 
NA, 0.486023134629369, NA, 0.489281948774758, NA, 0.464851913610455, 
NA, 0.499546896146173, NA, NA, 0.6555264303958, 0.605687762634627, 
NA, 0.683435292388886, NA, 0.615207645374612, NA, 0.739018717912613, 
NA, 0.79341715750565, NA, 0.769473075819677, NA, NA, NA, NA, 
NA, 0.6555264303958, 0.605687762634627, NA, 0.683435292388886, 
NA, 0.615207645374612, NA, 0.739018717912613, NA, 0.79341715750565, 
NA, 0.769473075819677, NA, 0.773437871034822, NA, NA, NA, 0.6555264303958, 
0.605687762634627, NA, 0.683435292388886, NA, 0.615207645374612, 
NA, 0.739018717912613, NA, 0.79341715750565, NA, 0.769473075819677, 
NA, 0.773437871034822, NA, 0.760395526989779, 0.607323723612999, 
0.6555264303958, 0.605687762634627, NA, 0.683435292388886, NA, 
0.615207645374612, NA, 0.739018717912613, NA, 0.79341715750565, 
NA, 0.769473075819677, NA, 0.773437871034822, NA, 0.760395526989779, 
NA), .Dim = c(17L, 6L), .Dimnames = list(c("1", "3", "4", "5", 
"6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", 
"17", "19"), c("1", "2", "3", "4", "5", "6")))

表格的清晰格式如下所示：

           1         2         3         4         5         6
1  0.6555264 0.6555264 0.6555264 0.6555264 0.6555264 0.6555264
3  0.6056878 0.6056878 0.6056878 0.6056878 0.6056878 0.6056878
4         NA 0.5768758        NA        NA        NA        NA
5  0.6834353        NA 0.6834353 0.6834353 0.6834353 0.6834353
6         NA 0.5028490        NA        NA        NA        NA
7  0.6152076        NA 0.6152076 0.6152076 0.6152076 0.6152076
8         NA 0.5160288        NA        NA        NA        NA
9  0.7390187        NA 0.7390187 0.7390187 0.7390187 0.7390187
10        NA 0.4860231        NA        NA        NA        NA
11 0.7934172        NA 0.7934172 0.7934172 0.7934172 0.7934172
12        NA 0.4892819        NA        NA        NA        NA
13 0.7694731        NA 0.7694731 0.7694731 0.7694731 0.7694731
14        NA 0.4648519        NA        NA        NA        NA
15        NA        NA        NA 0.7734379 0.7734379 0.7734379
16        NA 0.4995469        NA        NA        NA        NA
17        NA        NA        NA        NA 0.7603955 0.7603955
19        NA        NA        NA        NA 0.6073237        NA

这里的列标题是 Tube.numbers，行标题是 Measurement.locations。

我想将每个管号（所以 6 条线）绘制成一张图，x 轴是测量位置，y 轴是标准偏差（表格中的计算值）

我尝试使用 matplot，但我似乎无法绕过必须具有相同数量的行和列的必要性。我更喜欢循环执行。

谁能帮帮我？

非常感谢！

【问题讨论】：

通过调用dput(SDmeasurement) 并将其输出复制到您的问题中来提供您获得的exact output 对象可能会有所帮助。
... 或者提供一些代码来创建一个玩具示例供我们使用。此外，在 R 中，循环主要用作最后的手段。这里绝对不需要。
对不起各位，这还是个初学者，这对解决我的问题有帮助吗？感谢您的回复！
请注意，您的大多数 Tube 编号对于不同的测量位置具有相同的 SD 值，并且相应的线将被绘制在彼此的顶部。只有添加的最后一行在最终图中可见。

标签： r plot multiple-columns

【解决方案1】：

# Make an empty plot
plot(0, 0, type="n", xlim=c(1, 19), ylim=range(SDMeasurement, na.rm=TRUE))

# Add the lines
for(i in seq(df)){
    y <- SDMeasurement[[i]]
    lines(which(!is.na(y)), na.omit(y), col=i)
}

【讨论】：

Hey Backlin 感谢您的回复，但是鉴于我提供的额外数据，现在有没有更简单的方法来制定这个？我不太明白你的脚本...
现在好点了吗？第一段只是解析数据。重点仍然存在：先用plot(..., type="n") 制作一个空图，然后逐一添加不包括NA 值的行。

【解决方案2】：

另一个选项是使用ggplot2 来获得一个额外的自动图例（因为您应该将它添加到带有基本图的手动图例中）。首先，您将数据以长格式放置。我在这里使用来自reshape2 的melt。然后你应该删除缺失值。一种选择是使用complete.cases。这是我的代码：

library(reshape2)

dat.m <- melt(dat)
library(ggplot2)
ggplot(dat.m[complete.cases(dat.m),]) +
  geom_line(aes(x=Var1,y=value,color=factor(Var2))) +
  xlab("location")+ylab("Measurement")+
  guides(color = guide_legend(title = "Tube numbers"))

【讨论】：