[英]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.