a better simulation?

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在学长的帮助之下，完成对矫正破送回归的模拟验证和比较。以下是模拟数据。随后将尽快更新验证只结果。
正如我以前强调，利用SQL总结数据，然后结合某些过程里WEIGHT/FREQ STATEMENT的，使得很多过程的速度大为提高。
NOTE: variable ID was created in SQL procedure because GEE model needs ID variable.
京剧

1. data have;
   
2. call streaminit(12345);
   
3. array p[3,2,2,2,1] _temporary_( 0.40 0.40 0.2 0.2
   
4. 0.60 0.60 0.2 0.2
   
5. 0.80 0.40 0.4 0.2
   
6. 0.60 0.30 0.4 0.2
   
7. 0.90 0.30 0.6 0.2
   
8. 0.75 0.25 0.6 0.2);
   
9. do rr =1, 2, 3; /*---------true relative risk---*/
   
10. do r =1, 2; /*-------------------allocation---*/
    
11. do sampsize =100, 200, 500; /*sample size---*/
    
12. do rep =1 to 1000; /*--simulation times---*/
    
13. do i =1 to sampsize;
    
14. exposure =rand('bernoulli', 0.5);
    
15. if exposure =1 then strata =rand('tabled', 0.6);
    
16. else if exposure =0 then strata =rand('tabled', 0.4);
    
17. if strata =1 then do;
    
18. if exposure =1 then event =rand('bernoulli', p[rr,r,1,1,1]);
    
19. else if exposure =0 then event =rand('bernoulli', p[rr,r,1,2,1]);
    
20. end;
    
21. else if strata =2 then do;
    
22. if exposure =1 then event =rand('bernoulli', p[rr,r,2,1,1]);
    
23. else if exposure =0 then event =rand('bernoulli', p[rr,r,2,2,1]);
    
24. end;
    
25. output;
    
26. end;
    
27. end;
    
28. end;
    
29. end;
    
30. end;
    
31. run；
    
32. proc sql;
    
33. create table have2 as
    
34. select rr, r, sampsize, rep, strata, exposure, event, count(*) as count, monotonic() as id
    
35. from have
    
36. group by rr, r, sampsize, rep, strata, exposure, event
    
37. order by rr, r, sampsize, rep, strata, exposure desc, event desc;
    
38. quit;

复制代码

;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Simulation ulation better ATION ATI better

已经顺利完成simulation.
详细内容如下.
京剧

http://blog.sina.com.cn/s/blog_a3a926360102uyw1.html

jingju11 发表于 2014-8-12 05:43
已经顺利完成simulation.
详细内容如下.
京剧

更加详细的更新。京剧

http://blog.sina.com.cn/u/2745771574