Постройте только диагональные записи с данными. Кадр из многих временных рядов в R - PullRequest
0 голосов
/ 01 декабря 2011

У меня есть несколько временных рядов, которые, когда я строю график,

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 Ответ

4 голосов
/ 01 декабря 2011

Попробуйте:

 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=''))
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...