R で比較的単純な AR(1) 表現を書こうとしています。このコードには目立った問題は見つかりません。さらに、エラーを返しません。環境に書き込んでいないか、areone2 を関数。どんな提案でも大歓迎です。
areone2<-function(y,N,p,d){
yvec<-c(rep(y, times = N))
for(i in 1:N){
yvec[i+1]<-
((1+p*(yvec[i]-d))
+ d)
}
plot(yvec, type='l', xlab="N", ylab="yeild")
}
areone2(.3,10,.9,.2)