如何建立回归方程？？替换缺失值！！！

jingju11 发表于 2013-12-11 11:32
what does it mean?

对不起，这个说的有点乱！我重新回复了个贴。

1. %macro impute(VarList);       
   
2. data _null_;
   
3.         call symputx('nobs', nobs);
   
4.         if 0 then set have nobs =nobs;
   
5.         stop;
   
6.         run;
   
7.         %do i =1 %to %sysfunc(countw(&varList));
   
8.                 %let mvar =%scan(&varList, &i);
   
9.                 data have2;
   
10.                         do _n_=1 by 1 until(last.&mvar.);
    
11.                                 set have(keep=&mvar);
    
12.                                 by &mvar. notsorted;
    
13.                                 imp ++(&mvar. =0);
    
14.                                 end;
    
15.                         do _n_=1 by 1 until(last.&mvar.);
    
16.                                 set have(keep=&mvar);
    
17.                                 by &mvar. notsorted;
    
18.                                 &mvar.i =(imp >3);
    
19.                                 output;
    
20.                                 end;
    
21.                         imp =0;
    
22.                         keep &mvar. &mvar.i;
    
23.                         run;
    
24. 
    
25.                 data have3;
    
26.                         call streaminit(1234);
    
27.                         set have2 end=eof;
    
28.                         array xx[&nobs.,2] _temporary_;
    
29.                   xx[_n_,1] =&mvar.;
    
30.                         xx[_n_,2] =&mvar.i;
    
31.                   if eof then do;
    
32.                     if xx[1,1] =0 then do;
    
33.                                         %*Fetch 1st non-missing value;
    
34.                       do i =1 to dim1(xx) by 1;
    
35.                         if xx[i,1] ^=0 then leave;
    
36.                         end;
    
37.                       xx[1,1] =xx[i,1];
    
38.                       end;
    
39.                     if xx[&nobs.,1] =0 then do;
    
40.                                         %*Fetch last non-missing value;
    
41.                       do i =dim1(xx) to 1 by -1;
    
42.                         if xx[i,1] ^=0 then leave;
    
43.                         end;
    
44.                       xx[&nobs.,1] =xx[i,1];
    
45.                       end;
    
46.                           do i =1 to dim1(xx);
    
47.                             if xx[i,1] ^=0 then PreMissing =0;
    
48.                             if xx[i,1] =0 then do;
    
49.                               if PreMissing =0 then do;
    
50.                                       PreMissing =1;
    
51.                                       xL=xx[i-1,1];
    
52.                                                         %*Fetch next non-missing value;
    
53.                                       do j =i+1 to dim(xx);
    
54.                                         if xx[j,1] ^=0 then leave;
    
55.                                         end;
    
56.                                       xR =xx[j,1];
    
57.                                       xx[i,1] =ifn(rand('table', 0.5) =1, xL, xR);
    
58.                                       end;
    
59.                                                 %*impute missing when previous one imputed;
    
60.                               if PreMissing =1 then xx[i,1] =ifn(rand('table', 0.5) =1, xL, xR);
    
61.                               end;
    
62.                                         %*Wont impute those >3 consecutive missings;
    
63.                             &mvar._imp =ifn(xx[i,2]=1, 0, xx[i,1]);
    
64.                             output;
    
65.                             end;   
    
66.                                 end;
    
67.                         keep &mvar._imp;
    
68.                         run;
    
69.                 %if &i =1 %then%do;
    
70.                         data have_f;
    
71.                                 merge have;
    
72.                                 run;
    
73.                         %end;
    
74.                 data have_f;
    
75.                         merge have_f have3;
    
76.                         run;
    
77.         %end;
    
78. %mend impute;
    
79. 
    
80. 
    
81. %impute(A1 A2 A3 B1 B2 B3 C1 C2 C3 D1 D2 D3 E1 E2 E3)
    
82. ;

复制代码

jingju11 发表于 2013-12-13 23:46

太厉害了~~~谢谢！大侠，交个朋友吧！我Q：869470467

make a mark.