stata 生成随机数 按照组别

现在数据库有50个城市，每个城市内有200个个体的观察变量。
需要：随机赋予50个城市0或1，25个城市为1， 25个城市是0. ie: 如果一个城市被随机赋予0，则城市内的所有个体都是0.
求助需要如何实现命令？
谢谢！

1. clear
   
2. set obs 50
   
3. gen city=_n
   
4. 
   
5. set seed 123
   
6. gen temp=runiform()
   
7. sort temp
   
8. gen bool=(_n<=25)  
   
9. 
   
10. expand 200
    
11. sort city
    
12. by city: gen ob_num=_n
    
13. 
    
14. drop temp

复制代码