R语言产生一阶自回归误差的线性回归模型

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

n=200;m=8
              for(i in 1:n)
              {
                library(MASS)
                set.seed(2)
                sigma<-matrix(c(rep(1,64)),8,8)
                sigma
                x<-mvrnorm(n=200,rep(0,8),sigma)
                x
                    for(t in 2:m)
                    {
         
                        e[i,t]<-0.3*e[i,t-1]+x[i,t]
                        y[i,t]<-0.5*y[i,t-1]+e[i,t]
                    }

              }
请问坛友这个程序哪里有问题啊，运行时总会出现下面的提示
错误: 找不到对象'e'

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

不仅找不到对象e，也找不到对象y。
另外library()写在了for循环内，是要加载MASS包200次吗？

swufeyumeng 发表于 2018-2-10 13:14
不仅找不到对象e，也找不到对象y。
另外library()写在了for循环内，是要加载MASS包200次吗？

额，library()写错地方了，除了这个问题，那为啥找不到e啊，是不是没有设定初始值啊，还是说有别的问题o(︶︿︶)o

数据小仙paper.Libra 发表于 2018-2-10 13:16
额，library()写错地方了，除了这个问题，那为啥找不到e啊，是不是没有设定初始值啊，还是说有别的问题o( ...

看这一条语句e[i,t]<-0.3*e[i,t-1]+x[i,t]，嵌套循环首先相当于执行e[1,2]<-0.3*e[1,1]+x[1,2]，现在x[1,2]是已知的，e[1,1]呢？总不能凭空产生吧。另外mvrnorm()也被执行了200次，如果不是一定要这么做，它也应该放在for循环外。

swufeyumeng 发表于 2018-2-10 18:41
看这一条语句e

好的，谢谢