我有几个时间序列,当我使用时,
plot(na.omit(d))
我得到以下图表:
而不是整个情节矩阵,我只想要对角线的那些。换句话说,我不想绘制不同时间序列的组合。
以下非常混乱的代码实现了我想要的,但我想知道是否有更简单,更紧凑的方法来实现这一目标:
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)))
}
答案 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=''))