2011-12-01 11 views
-2

J'ai plusieurs séries chronologiques, qui quand je conspire en utilisant,Plot que les entrées en diagonale avec data.frame de plusieurs séries chronologiques en R

plot(na.omit(d)) 

je reçois le graphique suivant:

enter image description here

Au lieu de cette matrice entière de graphiques, je ne veux que ceux de la diagonale. En d'autres termes, je ne veux pas tracer les combinaisons de séries temporelles différentes.

Le code suivant très salissant atteint ce que je veux, mais je me demandais s'il y avait plus simple, plus compacte façon d'y arriver:

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))) 

Peut-être que quelque chose de semblable à cela?

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

Je ne peux pas vraiment ni queue ni tête de cette question. Pouvez-vous essayer de l'éditer pour rendre votre problème plus clair? – joran

+0

@joran: c'est le bazar, laissez le meilleur résoudre le puzzle: D AndresT a déformé cela. Résolu. Pour moi, il est difficile de clarifier cette question, pouvez-vous maintenant voir ce que je voulais dire? – hhh

+1

Oui, maintenant je vois ce que tu voulais dire. Voir mes modifications pour un exemple de comment écrire une question claire. – joran

Répondre

4

Essayez:

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

Pour » ./Pictures/pacfCO2 essayer quelque chose comme:

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

Merci, heureusement il y a encore des utilisateurs lettrés: P Exactement. – hhh

Questions connexes