题目是测算某种艾滋病指标细胞的个数，变量有血清转化开始计算的时间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   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