R语言笔记
••••设置种子123使用Y=rnorm(n,u,r)生成白噪声绘制yt（both方式）绘制其与y=0水平线的比较图fun1<-function(){set.seed(123)y<-rnorm(50,0,1)x<-c(1:50)plot(x,y,type=\"b\")abline(h=0)}••计算自相关系数•acf(yt)$ac检验该序列的自相关性•Box.test(yt)fun2<-function(){set.seed(123)y<-rnorm(50)*10+10a<-acf(y)a$ac#Box.test(y)}••••生成y=rnorm(50)*10+10生成y(t+1)=y(t)*2+rnorm(1)生成y(t+2)=y(t+1)+y(t)*2+rnorm(1)针对三个序列，分别•计算其自相关系数•检验该序列的自相关性fun3<-function(){set.seed(123)y<-c()y[1]<-1for(iin2:50){y[i]<-y[i-1]*2+rnorm(1)}a<-acf(y)a$acBox.test(y)}fun4<-function(){set.seed(123)y<-c()y[1]<-1y[2]<-1for(iin3:50){y[i]<-y[i-1]+y[i-2]*2+rnorm(1)}a<-acf(y)a$acBox.test(y)}••生成y=rnorm(50)*10+10生成y(t+1)=y(t)*20+rnorm(1)••生成y(t+2)=y(t+1)*10-y(t)*2+rnorm(1)针对三个序列，分别•计算其自回归fun51<-function(){y<-rnorm(50)*10+10plot(y)x11()pacf(y)LM<-lm(y~c(1:50))print(LM)}fun52<-function(n){y<-array(n)y[1]<-1for(iin2:n){y[i]<-20*y[i-1]+rnorm(1)}plot(y)x11()pacf(y)LM<-lm(y~c(1:50))print(LM)}fun53<-function(n){y<-array(n)y[1]<-1y[2]<-1for(iin3:n){y[i]<-10*y[i-1]-2*y[i-2]+rnorm(1)}plot(y)x11()pacf(y)LM<-lm(y~c(1:50))print(LM)}