簡體   English   中英

R中的情節循環

[英]Loop of plots in R

我有一個N C T數組

我希望為每個行(即數組中的每個N)繪制一個時間序列acf和pacf

我可以在一個例子中做到這一點

par(mfrow=c(3,1))
x <- x_whole[N,C,(1:T)]
plot(x, t='l')
abline(h=mean(x), col='red')
acf(x, lag.max=T-10)
pacf(x, lag.max=T-10)

現在,我想創建一個循環,對每個N

另外一些N為0,如果是這種情況,我希望將其忽略/跳過

我正在使用R軟件

我最大的努力是

for(i in 1:N){
        if(i != 0){
             par(mfrow=c(3,1))
            plot(x_whole[i,1,(1:T)], t='l')
        abline(h=mean(x_whole[i,1,(1:T)]), col='red')
        acf(x_whole[i,1,(1:T)], lag.max=190)
        pacf(x_whole[i,1,(1:T)], lag.max=190)

    }
    else{
        return 0 
        }
} 

數據看起來像

[,,1]
3.45    4    5    6    7
0       2    3    6    6
2.34    4    5    6    5

[,,2] 
4.56    4    5    6    7
4.43    2    3    6    6
0       4    5    6    5 

等等

我希望有一個3x1的圖,以[1,1,(1:T)],[2,1,(1:T)] [3,1,(1:T)]的圖作為輸出

同樣重要的是要注意第3、4、5列中的數據已排序

嘗試

for(i in 1:length(N)){
  if(i==0)
    {return 0} else
      {par(mfrow=c(3,1))
        plot(x_whole[i,1,(1:T)], t='l')
    abline(h=mean(x_whole[i,1,(1:T)]), col='red')
    acf(x_whole[i,1,(1:T)], lag.max=190)
    pacf(x_whole[i,1,(1:T)], lag.max=190)}}

如果沒有您的數據樣本,我將無能為力。 能否以data<-c(...)data<-matrix(...,ncol=...,nrow=...)的形式作為矢量/矩陣提供它,謝謝:)

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM