题目是测算某种艾滋病指标细胞的个数，变量有血清转化开始计算的时间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，下列数据只是一部分，完整数据在附件

）：
-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 518 10029
-0.7419581102 2.64 0 1 518 10029
-0.251882 801 2.64 0 1 538 10029
0.251882 824 2.64 0 1 5 7 10029
0.769336 866 2.64 0 1 515 10029
1.412731 704 2.64 0 1 521 10029
1.806982 757 2.64 0 1 525 10029
2.420260 726 2.64 0 1 529 10029
-1.3935661277 11.28 3 1-4-7 10039
-0.7200551132 11.28 3 0-2-5 10039
-0.2600961454 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