关于random effect自由度的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问线性mixed模型中，random effect的自由度应该如何计算？

比如下面程序生成一个简单的数据集，总共10个人，分为两组，每组5人，两组的平均数有差别，正态模型

data d1;
do i=1 to 10;
  group=int((i-0.000001)/5);
  y=group+rannor(0);
  output;
end;
run;

用PROC MIXED分析，
PROC MIXED DATA=d1;
  CLASS group i;
  MODEL y=group /noint s;
  RANDOM INTERCEPT /s SUB=i;
run;

得到以下结果：
Type 3 Tests of Fixed Effects
                                                    Num     Den
                                         Effect      DF      DF    F Value    Pr > F
                                         group        2       0       1.69

F检验的分母自由度为0，我理解这是因为RANDOM INTERCEPT/s SUB=i 这句为每个人都估计一个intercept，所以把自由度用光了。但是不太理解为什么会这样？
设这些intercept服从N(0,sigma)的分布，在fit model的时候， 应该是只有sigma被带入运算，也就是说，虽然intercept有10项之多， 但是实际上只贡献了一个参数而已。 一旦这个参数sigma定下来，具体每个人的intercept可以根据公式用sigma和y算出来。 那么为什么这些intercept会占用这么多的自由度呢？而不是只占一个自由度而已？

是我什么地方理解的不对么？请高手指教。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Effect random rand FECT 自由度 Effect 自由度 random