我运行以下程序
data a;
input TRIAL VD VWD NVD NVWD LATITUDE YEAR ALLOC;
LN_OR=log((VD/VWD)/(NVD/NVWD));
EST=1/VD+1/VWD+1/NVD+1/NVWD;
datalines;
1 4 119 11 128 44 48 1
2 6 300 29 274 55 49 1
3 3 228 11 209 42 60 1
4 62 13536 248 12619 52 77 1
5 33 5036 47 5761 13 73 2
6 180 1361 372 1079 44 53 2
7 8 2537 10 619 19 73 1
8 505 87886 499 87892 13 80 1
9 29 7470 45 7232 27 68 1
10 17 1699 65 1600 42 61 3
11 186 50448 141 27197 18 74 3
12 5 2493 3 2338 33 69 3
13 27 16886 29 17825 33 76 3
;
proc print;run;
Proc mixed method =ml data=a;
class trial;
model ln_or=/ s cl;
repeated /group=trial;
parms / parmsdata=a eqcons=1 to 13;
run;
data b;
set a;
keep est;
run;
Proc mixed cl method =ml data=a;
class trial;
model ln_or=/ s cl;
repeated /group=trial;
random int/ subject=trial s;
parms (0.01 to 2.00 by 0.01)(0.35712)
(0.20813)(0.43341)(0.02031)(0.05195)
(0.00991)(0.22701)(0.00401)(0.05698)
(0.07542)(0.01253)(0.53416)(0.07164)
/eqcons=2 to 14;
run;
得到结果,我如何才能把我最想要的 “Solution for Fixed Effects”这一行储存于一个数据集中呢?要添加什么语句呢?多谢高人指点!!!