data bq2;
input p1 p2 p3 p4 p5@@;
datalines;
68 73 85 88 86
92 69 74 65 83
88 76 76 70 80
81 73 84 98 94
83 79 95 83 98
84 67 86 56 66
76 76 68 64 86
53 96 65 95 94
. 97 76 87 64
66 93 80 90 73
85 95 81 81 69
78 66 99 90 71
58 86 72 63 81
94 84 70 78 86
94 81 80 66 92
93 66 91 74 97
63 74 90 63 92
91 79 83 85 84
94 95 64 96 95
56 67 91 97 56
61 80 79 70 69
86 96 79 84 75
69 90 65 65 76
92 85 82 66 68
68 . 65 84 87
71 66 61 75 94
61 74 76 87 78
63 80 69 76 84
86 68 95 71 84
64 83 61 90 96
60 85 96 67 87
82 84 97 78 60
88 92 66 59 95
60 91 78 78 81
59 97 75 76 88
65 87 86 64 96
84 78 83 61 85
65 93 62 99 83
92 99 79 86 90
84 82 92 95 76
94 90 65 66 84
90 79 85 81 58
67 89 84 75 93
63 82 65 69 66
85 97 83 84 70
86 76 64 87 69
88 88 96 80 87
62 98 74 93 62
80 93 85 82 72
87 84 80 93 64
94 85 94 74 93
55 75 93 84 60
90 68 88 92 83
59 95 69 75 74
98 63 80 63 84
93 55 66 84 96
75 64 65 94 63
63 94 . 82 76
71 82 61 57 61
55 72 95 85 64
86 55 67 62 80
51 65 78 94 80
81 94 73 63 95
90 63 95 91 87
60 83 64 79 83
74 94 96 89 76
63 74 91 94 83
58 63 84 84 72
68 93 91 82 91
70 83 75 96 76
86 73 73 75 94
97 83 97 64 68
78 81 87 78 69
63 71 92 86 68
67 82 87 63 86
91 73 90 79 74
63 93 97 90 76
87 83 65 91 68
65 84 73 87 98
78 64 82 85 90
81 92 65 77 82
90 82 92 66 90
64 73 84 58 76
78 94 77 67 95
61 84 75 69 72
90 93 72 94 73
93 73 83 90 90
69 72 88 94 74
88 63 88 76 66
76 56 72 75 82
82 74 94 89 87
60 65 84 85 73
75 84 66 70 75
79 74 78 63 85
74 64 91 94 79
70 55 95 83 69
93 94 74 73 85
85 83 79 95 71
81 63 70 79 95
86 85 92 87 74
92 78 85 70 93
;
proc mi data=bq2 seed=1000 out=putexp;
mcmc chain=multiple displayinit initial=em(itprint);
var p1 p2 p3 p4 p5;
run;
proc mi data =bq2 seed=6413 norprint nimpute=2;
mcmc acfplot(mean(p1)/symbol=dot lref=2);
var p1 p2 p3 p4 p5;
run;
proc corr data=putexp cov nocorr noprint out=outcov(type=cov);
var p1 p2 p3 p4 p5;
by_Imputation_;
run;
proc print data=outcov(obs=12);
title 'CORR Means and Covariance Matrices'
'(First Two Imputations)';
run;
proc mianalyze data=outcov edf=32 wcov bcov tcov mult;
modeleffects p1 p2 p3 p4 p5;
run;
有大神,可以帮忙解答下吗?程序的结果一直出不来也不报错