我正在尝试在 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)