R软件：ar()函数识别已经模拟生成的AR（1）问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

#程序
AR1_phi=-0.2
n=50;
set.seed(123456)
for(i in 1:10)                        #重复10遍
{
n0=1000                             #这1000个之前的是要删掉的，我们要的是后面n=50数据
n1=n0+n
z=rnorm(n1,mean=0,sd=1)          #生成1050个服从标准正态的误差z,但我们只要从1001到1050的后面50个就是下面的zt
x=filter(z,filter=phi,"recursive")      #生成AR(1)序列x(x里有1050个数，但我们只要从1000到1050的后面51个就是下面的xt)
zt=z[(n0+1):n1]                  
xt=x[n0:n1]
b=ar(xt, order.max=1,AIC=F,method= 'yw')      #再用ar函数识别刚才生成的xt
a=b[[2]]                                                              #a 使用Y-W方法估计真实的ph（i=-0.2）的估计量
print(a)      
}
#结果

[1] -0.5441331
[1] -0.2585425
[1] -0.3503617
[1] -0.2825224
numeric(0)
[1] -0.4107991
[1] -0.3618471
[1] -0.2454404
numeric(0)
[1] -0.2150401
#为什么会出现两个numeric(0)？求大神，我要崩溃了。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：r软件 Recursive numeric filter Method filter method 软件 程序

对不起大家，有个地方：把AR1_phi=-0.2改成phi=-0.2啊