仅绘制具有R中许多时间序列的data.frame的对角线条目

时间:2011-12-01 05:50:13

标签: r plot

我有几个时间序列,当我使用时,

plot(na.omit(d))

我得到以下图表:

enter image description here

而不是整个情节矩阵,我只想要对角线的那些。换句话说,我不想绘制不同时间序列的组合。

以下非常混乱的代码实现了我想要的,但我想知道是否有更简单,更紧凑的方法来实现这一目标:

png('./Pictures/acf1.png')
plot(acf(na.omit(f$C1)))
png('./Pictures/acf2.png')
plot(acf(na.omit(f$C2)))
png('./Pictures/acf3.png')
plot(acf(na.omit(f$C3)))
png('./Pictures/acf4.png')
plot(acf(na.omit(f$C4)))
png('./Pictures/acf5.png')
plot(acf(na.omit(f$C5)))
png('./Pictures/acfCO2.png')
plot(acf(na.omit(f$MLCO2)))

png('./Pictures/pacf1.png')
plot(pacf(na.omit(f$C1)))
png('./Pictures/pacf2.png')
plot(pacf(na.omit(f$C2)))
png('./Pictures/pacf3.png')
plot(pacf(na.omit(f$C3)))
png('./Pictures/pacf4.png')
plot(pacf(na.omit(f$C4)))
png('./Pictures/pacf5.png')
plot(pacf(na.omit(f$C5)))
png('./Pictures/pacfCO2.png')
plot(pacf(na.omit(f$MLCO2)))

也许类似的东西?

for(i in d) 
{
    png(PACF$i.png)
    plot(pacf(na.omit(d$i)))
}

1 个答案:

答案 0 :(得分:4)

尝试:

 for(i in 1:ncol(d)){
png(paste('c:/zzz/',i,'.png',sep =''))
plot(pacf(na.omit(d[,i])))
dev.off()
}

对于'./Pictures/pacfCO2尝试类似:

png(paste('./Pictures/pacf',i,'.png',sep=''))