楼主: hopui2017
457 1

[SAS EM] SAS或者R解决非线性模型非等距2阶自相关消除CAR(2)的问题 [推广有奖]

  • 8关注
  • 0粉丝

大专生

80%

还不是VIP/贵宾

-

威望
0
论坛币
1044 个
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
150 点
帖子
30
精华
0
在线时间
73 小时
注册时间
2017-10-19
最后登录
2019-1-15

hopui2017 发表于 2018-12-9 21:21:37 |显示全部楼层
200论坛币
本帖最后由 hopui2017 于 2018-12-9 21:21 编辑

大家好,向大家请教一下如何使用R或者SAS消除非线性模型中非等距的2阶后滞CAR(2)。我的模型形式如下图所示 koz.JPG
其中d为因变量,D,H,h为因变量,其中h由于是同一个体非等距测量,所以才同一个个体中存在自相关性。通过查阅文献,消除自相关需要加入2阶后滞才能完全消除,残差公式表示如下。
zxg.JPG


我要想达到的效果是
zxgt.JPG



我的数据结构如下
sj.JPG

。其中bolt就是非线性模型中的h,tht是H,dbh是D.
我利用sas加入了一阶后滞CAR(1)后模型收敛,但是加入CAR(2)时模型参数p2不显著。
加入CAR(2)  SAS代码如下
PROC model data=fit_taper method=marquardt sur dw collin;
exogenous bolt tht dbh;
endogenous dob  ;
parms b0 0.9884 b1 0.9478 b2 0.0735 b3 0.4884 b4 0.2783
b5 0.5511 b6 -1.1  b7 -0.0389 b8 0.2579 p1 0.6 p2 0.001;/*表示赋初始参数*/
dob=b0*(dbh**b1)*(tht**b2)*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)
**(1/3)))**(b3*(bolt/tht)**4+b4*(1/exp(dbh/tht))+b5*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)**(1/3)))**0.1+b6*(1/dbh)+b7*tht**(1-(bolt/tht)**(1/3))+b8*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)**(1/3))))

+(J>1)*p1**(bolt-zlag1(bolt))*zlag1(dob-(b0*(dbh**b1)*(tht**b2)*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)
**(1/3)))**(b3*(bolt/tht)**4+b4*(1/exp(dbh/tht))
+b5*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)**(1/3)))**0.1
+b6*(1/dbh)+b7*tht**(1-(bolt/tht)**(1/3))
+b8*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)**(1/3))))))
+(J>2)*p2**(bolt-zlag2(bolt))*zlag2(dob-(b0*(dbh**b1)*(tht**b2)*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)
**(1/3)))**(b3*(bolt/tht)**4+b4*(1/exp(dbh/tht))+b5*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)**(1/3)))**0.1+b6*(1/dbh)+b7*tht**(1-(bolt/tht)**(1/3))+b8*((1-(bolt/tht)**(1/3))/(1-(1.3/tht)**(1/3)))))
;
fit dob  ;
run;    /*表示结束运算*/
其中红色部分表示原非线性模型,蓝色的依照消除自相关残差公式的2阶后滞。结果显示模型的p2不显著。所以有没有人会这方面的知识,能帮忙解决吗?谢谢





stata SPSS
hopui2017 发表于 2019-1-9 22:18:25 |显示全部楼层
还没解决!
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 我要注册

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2019-3-20 07:48