文档介绍:#——1999年中国年纱产量序列时序图〔〕”,header=T)#如果有标题,用T;没有标题用Fplot(,type='o')=[,2])x=rep(0,1500)x[1]=x1x[2]=x1-*x0+rnorm(1)for(iin3:length(x)){x[i]=x[i-1]-*x[i-2]+rnorm(1)}plot(x[1:100],type=T)acf(x)pacf(x)x0=runif(1)x1=runif(1)x=rep(0,1500)x[1]=x1x[2]=x1+*x0+rnorm(1)for(iin3:length(x)){x[i]=x[i-1]+*x[i-2]+rnorm(1)}plot(x[1:100],type=T)acf(x)pacf(x)又一个式子x0=runif(1)x1=runif(1)x=rep(0,1500)x[1]=x1x[2]=-x1-*x0+rnorm(1)for(iin3:length(x)){x[i]=-x[i-1]-*x[i-2]+rnorm(1)}plot(x[1:100],type=T)acf(x)pacf(x)#均值和方差smu=mean(x)svar=var(x)mu=0mvar=1/(1-)#书上51页#总体均值方差cat("populationmeanandvarare”,c(mu,mvar),"\n")#样本均值方差cat("samplemeanandvarare",c(mu,mvar),"\n")svar=(1+)/((1-)*(1-1-)*(1+1-))#:
x=(n=1000,list(ma=-2))(x,type='l')acf(x)pacf(x)法二x=rep(0:1000)for(iin1:1000){x[i]=rnorm[i]-2*rnorm[i-1]}plot(x,type='l')acf(x)pacf(x)法一:
x=(n=1000,list(ma=-))(x,type='l')acf(x)pacf(x)法二x=rep(0:1000)for(iin1:1000){x[i]=rnorm[i]-*rnorm[i-1]}plot(x,type='l')acf(x)pacf(x)##错误于rnorm[i]:类别为'closure'的对象不可以取子法一:
x=(n=1000,list(ma=c(-4/5,16/25)))(x,type=T)acf(x)pacf(x)法二:
x=rep(0:1000)for(iin1:1000){x[i]=rnorm[i]-4/5*rnorm[i-1]+16/25*rnorm[i-2]}plot(x,type='l')acf(x)pacf(x)##错误于x[i]=rnorm[i]-4/5*rnorm[i-1]+16/25*rnorm[i-2]:
##更换参数长度为零##例3.