Stata 排列组合 完全枚举从k个数中选取c个数的分组情况和组合模式

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. /*
   
2. 问题需求：
   
3. 一个地区16个村，8个村为干预，8个村为对照，
   
4. 考察所有随机化分组的基线变量平衡情况
   
5. 从K个村里任选c个村，有k!/(c!*(k-c)!)种
   
6. 
   
7. 解决思路：
   
8. 我们只需要知道干预组和对照组，用0和1表示即可，这样可以不用计算16的阶乘
   
9. 只用考虑16个变量，每个变量包含2个元素1和0，这16个变量的组合，2的16次方
   
10. 然后再删除行合计不等于8的即可
    
11. 
    
12. 优势：
    
13. 1. 突破percom等命令包中因为排列计算过大对选取数c的限制
    
14. 2. 突破被选取数k太大，造成的内存不足限制
    
15. 3. 可以生成原始数值的组合或者指示变量的分组
    
16. */
    
17. 
    
18. clear
    
19. local k=6 // 从k个数中一次性选取c个数
    
20. local c=2
    
21. di "组合数:" exp(lnfactorial(`k'))/(exp(lnfactorial(`c'))*exp(lnfactorial(`k'-`c')))
    
22. 
    
23. set obs 2
    
24. forvalues i=1/`k'{
    
25.     gen v`i'=1 in 1
    
26.         replace v`i'=0 in 2
    
27. }
    
28. di "记录数:" 2^`k'
    
29. fillin v* // 2^k条记录数，受电脑内存限制
    
30. 
    
31. // 以下代码生成指示变量的分组
    
32. egen rt=rowtotal(v1-v`k')
    
33. drop if rt!=`c' // 删除不需要的组合数
    
34. drop _fillin rt
    
35. d,s
    
36. 
    
37. // 以下代码继续生成带原始数值的组合
    
38. forvalues i=1/`k' {
    
39.     gen t`i'="`i'" if v`i'==1
    
40. }
    
41. 
    
42. egen cv=concat(t1-t`k'),punct(" ")
    
43. replace cv=strtrim(stritrim(cv))
    
44. keep cv
    
45. split cv,generate(v) destring
    
46. drop cv
    
47. compress

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Stata tata Factorial forvalues forvalue