多重插值-SAS中的MI步

yuanyuan1011

我写了一个程序，SAS运行的时候只占4000多字节的内存，不知道有没有跑，过了几天都没出结果，但编辑器里一直显示running.
这个程序是这样的：（copy可直接运行）
data mydata;
input a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11 a12;
datalines;

a =
       1.0438      0.54232       1.4816      0.63658      0.78846      0.25464     -0.16252       1.1151       1.9974      0.25464        1.581      0.22314
       1.8825            .       1.9315       1.5261       1.7281       1.1249      0.85442       1.1151       2.2006       1.1249       1.9796       1.3762
       2.3514       2.1175       2.1041       2.1518       2.4336            .       1.6134       1.9892            .       1.8625       2.0794       2.0567
            .       2.1199       2.0541       2.1258       2.2332       1.8437        1.884       2.0229       1.8453       1.8437       1.9184            .
       2.1494       1.9242       2.0149       2.0202       2.1679       1.7102       1.9587       1.8856       1.7334       1.7102       1.7699       2.2773
            .        1.805       1.8245       1.9286       2.0229       1.5974       1.8961       1.7716       1.7352       1.5974       1.6525       2.1483
       2.0109       1.6864       1.6677       1.7544       1.9587       1.3913       1.6582       1.5539       1.4446       1.3913       1.4929       1.8856
       1.9301       1.5151       1.5892            .       1.775        1.2413       1.4793       1.5195       1.4134       1.2413       1.2865       1.8099
       1.7817       1.1019       1.3083       1.4327       1.4748       1.0225       1.2613       1.0986       1.1506            .      0.98954       1.5195
       1.1878     -0.10536      0.04879      0.13976      0.45108    -0.083382      0.13976      0.22314      0.11333    -0.083382     -0.15082        0.157
;

PROC MI  
        DATA=mydata  OUT=impute_theophylline
        SEED=1000  NIMPUTE=5;


其实是依照一篇论文里的程序做的，论文里的源代码是：
data example;
input oxygen time rate @@;
datalines;
43.509 10.27 175 44.313 10.14 182
55.018 8.73 163 59.571 . .
48.982 9.34 . 45.823 12.14 175
. 12.04 177 . 11.02 .
39.501 12.97 175 60.247 9.01 171
50.724 . . 37.413 13.92 183
45.011 11.34 177 46.924 . .
52.145 10.48 167 48.877 9.04 178
41.241 11.02 171 456.721 10.13 .
46.835 10.37 . 50.117 9.89 165
40.216 12.71 175 46.136 11.34 158
45.532 9.71 165 . 9.71 .
45.357 11.19 . 40.014 13.24 167
46.121 10.68 187 50.623 9.89 151
48.714 9.52 186 48.073 11.65 171
47.542 10.61 172
;
proc mi data=example seed=1000 out=outexp;
var oxygen time rate;
run;
也是拷贝直接运行，但这个源程序可以直接出结果。我的错在哪了呢。。。

       

hy1011

为什么没有高手理我呢。。