大神来帮忙看下GLMM模型，用PROC MIXED来做，有问题啊，求指点：

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

题目是测算某种艾滋病指标细胞的个数，变量有血清转化开始计算的时间T，和细胞的假定年龄AGE，一共有2376个观测值，来自369个病人，模型是这样的：
Y(T)= β0+β1T+β2T²+β3AGE+b0+b1T+ε
b0+b1T是随即项，
我的代码是这样的：
DATA HW;
drop x4_x7;
infile 'D:\MACS.DATA';
input time count age x4-x7 id;
time2=time**2;
run;
proc mixed data=HW ;;
class id;
model count = time time2 age/s;
random time;
run;

由于T在FIXED项中 又在RANDOM项中 不知道为什么我算出来T的自由度为0，总感觉不对，跟同班的师姐对了一下（师姐用R自己编的），估计参数不一样。



数据结构大概是这样的（第一列为T，第二列为细胞个数，第三列为AGE，中间几列没用，最后一列为病人ID，下列数据只是一部分，完整数据在附件 data.txt (120.66 KB)

2012-10-24 10:49:42 上传

）：
   -0.741958   548    6.57   0   0   5   8   10002
   -0.246407   893    6.57   0   1   5   2   10002
    0.243669   657    6.57   0   1   5  -1   10002
   -2.729637   464    6.95   0   1   5   4   10005
   -2.250513   845    6.95   0   1   5  -4   10005
   -0.221766   752    6.95   0   1   5  -5   10005
    0.221766   459    6.95   0   1   5   2   10005
    0.774812   181    6.95   0   1   5  -3   10005
    1.256673   434    6.95   0   1   5  -7   10005
   -1.240246   846    2.64   0   1   5  18   10029
   -0.741958  1102    2.64   0   1   5  18   10029
   -0.251882   801    2.64   0   1   5  38   10029
    0.251882   824    2.64   0   1   5   7   10029
    0.769336   866    2.64   0   1   5  15   10029
    1.412731   704    2.64   0   1   5  21   10029
    1.806982   757    2.64   0   1   5  25   10029
    2.420260   726    2.64   0   1   5  29   10029
   -1.393566  1277   11.28   3   1  -4  -7   10039
   -0.720055  1132   11.28   3   0  -2  -5   10039
   -0.260096  1454   11.28   3   1  -3  -6   10039
    0.260096   738   11.28   3   0  -4  -7   10039
   -0.306639   994   17.99   0   1   5  -7   10048
    0.306639   486   17.99   0   1   5  -7   10048
    0.813142   605   17.99   0   1   5  -5   10048
    1.095140   880   17.99   0   1   5   7   10048
    1.593429   352   17.99   0   1   0  -7   10048
    2.094456   376   17.99   0   0  -4  -6   10048
    2.625599   212   17.99   0   0  -4   1   10048

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：proc mixed Mixed MM模型 GLMM ROC 模型

from your model specification, i would do the following:
random intercept time/subject = id;

I may be thinking of including time2 into random part as well.

jingju