2016-11-30 196 views
0

我已300个变量在10个时间点(行),在任何给定时间点对各变量I具有温度值A(列)和F.R线的曲线图,值的范围之外绘图区

附上的样品数据帧

structure(list(Timepoint = c(1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L, 
5L, 5L, 6L, 6L, 7L, 7L, 8L, 8L, 9L, 9L, 13L, 13L, 25L, 25L), 
Temperature = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 
1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), .Label = c("A", 
"F"), class = "factor"), Concentration.of.chylomicrons = c(1.29e-11, 
1.25e-11, 1.02e-11, 1.1e-11, 1.08e-11, 1.3e-11, 1.28e-11, 
1.26e-11, 1.06e-11, 1.32e-11, 8.85e-12, 1.21e-11, 8.83e-12, 
1.08e-11, 1.35e-11, 1.12e-11, 8.99e-12, 1.08e-11, 9.55e-12, 
1.04e-11, 0, 1.01e-11), Total.lipids = c(0.00268, 0.0026, 
0.00208, 0.00225, 0.00222, 0.0027, 0.00268, 0.0026, 0.00219, 
0.00273, 0.0018, 0.00247, 0.00179, 0.00221, 0.00276, 0.00229, 
0.00182, 0.00222, 0.00195, 0.00212, 0, 0.00204), Phospholipids = c(0.000224, 
0.000223, 0.000145, 0.00016, 0.000157, 0.000211, 0.00023, 
0.000211, 0.000165, 0.000224, 0.000109, 0.00018, 0.000113, 
0.000163, 0.000175, 0.000177, 0.000122, 0.000173, 0.000127, 
0.000156, 0, 0.000138)), .Names = c("Timepoint", "Temperature", 
"Concentration.of.chylomicrons", "Total.lipids", "Phospholipids" 
), class = "data.frame", row.names = c(NA, -22L)) 

我想绘制一个线图来展示每个变量如何随时间变化。在这个线条图上,我想要绘制A和F线。我已经成功地为此编写了循环代码。

# subset based on temperatures A and F 
a_df <- subset(df, Temperature == "A") 
f_df <- subset(df, Temperature == "F") 

# loop from columns 3:x 
for (i in 3:x) { 
    plot(a_df[, 1], 
     a_df[, i], 
     type = "l", 
     ylab = colnames(a_df[i]), 
     xlab = "Timepoint", 
     lwd = 2, 
     col = "blue") 
    lines(f_df[, 1], 
     f_df[, i], 
     type = "l", 
     lwd = 2, 
     col = "red") 
    legend("bottomleft", 
     col = c("blue", "red"), 
     legend = c("Temperature A", "Temperature F"), 
     lwd = 2, 
     y.intersp = 0.5, 
     bty = "n") 
} 

然而,对于某些变量,某些地方的绘图区之外,像下面附

Please click on this link for image 我怎样才能确保在这个循环中的命令,我可以有可见的所有点graghs。我确定有一个快速的方法来解决这个问题,任何人都可以帮忙吗?

我曾尝试下面的行,请建议 ylim = C(分钟(f_df [, - 1],最大值(f_df [, - 1]),

我得到以下错误消息

for (i in 3:229) { + plot(a_df[, 1], + a_df[, i], + type = "b", + ylim = c(min(f_df[,-1] ,max(f_df[,-1]), + ylab = colnames(f_df[i]), + main = colnames(f_df[i]), + xlab = "Timepoint", + lwd = 2, + col = "red") + lines(f_df[, 1], Error: unexpected symbol in: " col = "red") lines" f_df[, i], Error: unexpected ',' in " f_df[, i]," type = "b", Error: unexpected ',' in " type = "b"," lwd = 2, Error: unexpected ',' in " lwd = 2," col = "blue") Error: unexpected ')' in " col = "blue")" legend("bottomleft", + col = c("red", "blue"), + legend = c("Ambient", "Fridge"), + lwd = 2, + y.intersp = 0.5, + bty = "n") Error in strwidth(legend, units = "user", cex = cex, font = text.font) : plot.new has not been called yet } Error: unexpected '}' in "}"

Lakmal

+0

在你的绘图命令中尝试'ylim = c(min(f_df [, - 1],max(f_df [, - 1])? – dww

+0

在我写入该行的命令中,我得到错误没有任何地块 – NLM09

+0

它是一个想法将plot限制添加到for循环中的“plot”函数以匹配数据的范围吗?例如:'ylim = c(min(a_df [,i],f_df [ ,i]),max(a_df [,i],f_df [,i])),编辑:与dww相同的建议,但我认为你需要在两个数据集中指定范围作为最小值和最大值。把这个命令放到'plot'函数中就可以工作了。 – Niek

回答

0

要在回答中回顾一下。设置ylim解决问题

# loop from columns 3:x 
for (i in 3:x) { 
    plot(a_df[, 1], 
     a_df[, i], 
     type = "l", 
     ylab = colnames(a_df[i]), 
     xlab = "Timepoint", 
     ylim = c(min(df[,-1]) ,max(df[,-1])), 
     lwd = 2, 
     col = "blue") 
... 

将每个地块的地块边界设置为相等,如果您想比较地块,那么地块面积可能会比您的数据大得多。

# loop from columns 3:x 
for (i in 3:x) { 
    plot(a_df[, 1], 
     a_df[, i], 
     type = "l", 
     ylab = colnames(a_df[i]), 
     xlab = "Timepoint", 
     ylim = c(min(df[,i]) ,max(df[,i])), 
     lwd = 2, 
     col = "blue") 
... 

集对于每个小区,这是更差作比较,但减少了不必要的空积空间的新边界。我已将min(a_df[, i],f_df[, i])替换为min(df[,i]),因为它们应该是相同的。