如何实现随机抽取

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我想实现的就是每一个缺失值用其已观测到的值中随机抽取一个值进行填补，比如变量var1 一共有6个观测（1  2  .  4 . 6）两个缺失值，其中每个缺失值 都在已观测值（1 2 4 6）中随机抽取一个进行填补，想请教各位高手怎么实现，我自己有一个很麻烦的程序，运行很慢，所以想请教大家怎样能够简单一点

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：如何实现 缺失值 怎样能 观测值 VaR 如何 程序

如果你的数据集观测值很多的话建议这样做

1. data one;
   
2. input v1;
   
3. cards;
   
4. 1
   
5. 2
   
6. .
   
7. 4
   
8. .
   
9. 6
   
10. ;;;
    
11. run;
    
12. data two;
    
13. set one(where=(v1 ne .));
    
14. run;
    
15. 
    
16. data three;
    
17. do i =1 to 4; /*这里的数字4 是数据集two里面剩下的非空数据的数量*/
    
18. set two ;
    
19. array choose {4};
    
20. choose{i}=v1;
    
21. drop i v1;
    
22. end;
    
23. 
    
24. run;
    
25. data final;
    
26. if _n_=1 then set three;
    
27. set one;
    
28. if v1 eq . then do;
    
29. v1=choosen(ceil(ranuni(0)*4),of choose1-choose4);
    
30. end;
    
31. keep v1;
    
32. run;

复制代码

如果数据量不是很多的话 可以这样

1. data final2;
   
2. set one;
   
3. array choose {4} (1 2 4 6);
   
4. if v1 eq . then do;
   
5. v1=choosen(ceil(ranuni(0)*4),of choose1-choose4);
   
6. end;
   
7. keep v1;
   
8. run;

复制代码

注意所有的choose后面的数字以及 array的数字 都是非空数据的数量

sas那个石家庄人 发表于 2015-10-1 02:41
如果你的数据集观测值很多的话建议这样做

choose 函数用得很好，第一次见到。不过感觉第一个程序不能工作啊。

我的设想：将非missing的数据输出为另一个dataset, 从中抽样2个数据，然后补回原来的数据中

teqel 发表于 2015-10-1 08:24
choose 函数用得很好，第一次见到。不过感觉第一个程序不能工作啊。

我的设想：将非missing的数据输出 ...

我第一个程序可能有些地方没有写好 不过和你的设想思路是一致的

data two就是包含非missing的数据集

data three就是把data two里面非空的4个值变成一个单独的observation （choose1 choose2 choose3 choose4）

最后一步就是把这一列 （choose1 choose2 choose3 choose4）合并到原始数据的每一个观测值的右边
如果是missing就从这四个里随机选一个

sas那个石家庄人 发表于 2015-10-1 02:41
如果你的数据集观测值很多的话建议这样做

试了一下，可以工作。但是如果数据很多，会很慢

1. data have;
   
2.   input x @@;
   
3. datalines;
   
4. 1 2 3 . 5 . 6 7 8 . 10
   
5. ; run;
   
6. 
   
7. proc sql outobs=1;
   
8.     create table random_x as
   
9.     select x as ran_x
   
10.     from have
    
11.     where ^missing(x)
    
12.     order by ranuni(12);
    
13. quit;
    
14. 
    
15. data want;
    
16.     if _n_ = 1 then set random_x;
    
17.      set have;
    
18.      if missing(x) then x=ran_x;
    
19.     drop ran_x;
    
20. run;

复制代码

1. data one;
   
2. input v1;
   
3. cards;
   
4. 1
   
5. 2
   
6. .
   
7. 4
   
8. .
   
9. 6
   
10. ;;;
    
11. run;
    
12. 
    
13. data two;
    
14. set one(where=(v1 ne .));
    
15. num=_n_;
    
16. run;
    
17. 
    
18. 
    
19. data final;
    
20. 
    
21. array choose {4} _temporary_;
    
22. if _n_=1 then do until (last);
    
23. set two end=last;
    
24. choose{num}=v1;
    
25. end;
    
26. 
    
27. set one;
    
28. if v1 eq . then v1=choose{ceil(ranuni(0)*4)};
    
29. drop num;
    
30. 
    
31. run;
    
32. 
    
33. 我改进了一了下 加入了个temporary的array 应该可以减短运行时间

复制代码