如何求混合模型中随机参数值，非常感谢！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是某参考文献中附录的编程，是求混合模型中的随机参数值，其中Za1和Za2不知是什么，也不知道怎么求的，求好心人帮帮忙，万分感谢！！！！

%let n=2; /*n=number of tress used for calibration*/

Data example;

Input d h ddom hdom weight;/*weight=242.87/(area of the subplot), in the case of unequal selection probabilities;weight=1, in other case*/

H=1.3+(ht-1.3)*(exp(-7.914*d**(-1.444+0.02259*(ht-1.3)))/exp(7.914*dg**(-1.444+0.02259*(ht-1.3))));

Cards;

9 11 18 18 3.092

9.5 12 18 18 3.092

10 . 18 183.092

13 . 18 18 0.7731

15 . 18 18 0.7731

18 . 18 180.7731

;

Run;

Data calibration;

Set example;

Res=h-hhat;

B0=-7.914*d**(-1.444+0.02259*(ht-1.3));/*auxiliaryvariable*/

B1=-7.914*dg**(-1.444+0.02259*(ht-1.3));/*auxiliary variable*/

Za1=(ht-1.3)*exp(b0)*(b0*log(d)-b1*log(dg))/exp(b1);

Za2=(ht-1.3)**2*exp(b0)*(b0*log(d)-b1*log(dg))/exp(b1);

W=(1/weight)**0.5;

Where h is not null;

Run;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：非常感谢 随机参数 混合模型 calibration Selection 模型 如何

自己先顶一个，不要沉

Za是random effect的design matrix

yuerqieqie 发表于 2013-8-27 20:36
Za是random effect的design matrix

那如果是换一个公式，这个design matrix该如何设置，换个说法就是这个design matrix是怎么得来的呢

dongyfei 发表于 2013-8-28 09:18
那如果是换一个公式，这个design matrix该如何设置，换个说法就是这个design matrix是怎么得来的呢

楼主最好详细的说说这段程序是在干什么，这段code也不能运行，很难进一步帮你。

yuerqieqie 发表于 2013-8-29 03:36
可能前面说的有歧义。更正一下Za看上去像是random effect的design matrix前的系数，就像B（fix effect的d ...

这是完整的程序编程：
%let n=2; /*n=number of tress used for calibration*/
Data example;
Input d h ddom hdom weight;
/*weight =242.87/(area of the subplot), in the case of unequal selection probabilities; weight=1, in other case*/
Hhat=1.3+(hdom-1.3)*(exp(-7.914*d**(-1.444+0.02259*(hdom-1.3)))/exp(7.914*ddom**(-1.444+0.02259*(hdom-1.3))));
Cards;
9 11 18 18 3.092
9.5 12 18 18 3.092
10 . 18 18 3.092
13 . 18 18 0.7731
15 . 18 18 0.7731
18 . 18 18 0.7731
;
Run;
Data calibration;
Set example;
Res=h-hhat;
B0=-7.914*d**(-1.444+0.02259*(hdom-1.3));/*auxiliary variable*/
B1=-7.914*ddom**(-1.444+0.02259*(hdom-1.3));/*auxiliary variable*/
Za1=(hdom-1.3)*exp(b0)*(b0*log(d)-b1*log(ddom))/exp(b1);
Za2=(hdom-1.3)**2*exp(b0)*(b0*log(d)-b1*log(ddom))/exp(b1);
W=(1/weight)**0.5;
Where h is not null;
Run;
Proc iml;
Use calibration;
Read all var {za1 za2} into z;
Read all var {res} into res;
D={0.1512 -0.004188, -0.004188 0.0001244};
Read all var {w} into w;
G=w#i(&n);
R=2.561*g*i(&n)*g;
B=d*z`*inv(z*d*z`+r)*res;
Print b;
就是求混合模型中两个随机效应参数u和v的值，

dongyfei 发表于 2013-8-29 09:31
这是完整的程序编程：
%let n=2; /*n=number of tress used for calibration*/
Data example;

Za的确是random effect的design matrix。一般来说，design matrix对建模的人来说是已知的。它取决于你的假设条件，比如说我想对学生的成绩建模，假设以班级为单位有random effect。现在有两个班，每班2人，这个design matrix就设定为[1 0， 1 0， 0 1， 0 1]

yuerqieqie 发表于 2013-8-29 21:02
Za的确是random effect的design matrix。一般来说，design matrix对建模的人来说是已知的。它取决于你的假 ...

啊？是已知的吗，可是一般论文上都没有明确说明唉，混合模型建好了，不知道怎么检验