[Stata]模拟标准正态分布

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

[Stata]模拟标准正态分布

1. clear
   
2. set obs 10000
   
3. gen u1=uniform()   //生成均匀分布随机变量u1
   
4. gen u2=uniform()   //生成均匀分布随机变量u2
   
5. gen r_square=-2*log(u1)    //生成参数为2的指数分布R^2
   
6. gen theta=2*_pi*u2         // 生成（0，2pi)上的均匀分布theta
   
7. gen x=sqrt(r_square)*cos(theta)      //转换极坐标系
   
8. gen y=sqrt(r_square)*sin(theta)      
   
9. histogram x,bin(50) norm       // 50个直条
   
10. histogram y,bin(50) norm
    
11. pwcorr x y     //计算相关系数

复制代码

*标准正态分布模拟方法二（快速）

1. clear
   
2. set obs 12000
   
3. gen r_square=sqrt(-2*log(uniform()))
   
4. gen v1=uniform()*2-1      //分别生成（-1，1）上的均匀随机变量
   
5. gen v2=uniform()*2-1
   
6. keep if v1^2+v2^2 <=1    //只保留v1 v2单位元之内的部分，不必再进行转化
   
7. keep if _n <=10000    //限制保留在10000内
   
8. gen x=sqrt(r_square)*v1/sqrt(v1^2+v2^2)    //  v1/sqrt(v1^2+v2^2)是cos(theta)
   
9. gen y=sqrt(r_square)*v2/sqrt(v1^2+v2^2)    //   v2/sqrt(v1^2+v2^2)是sin(theta)
   
10. histogram x,bin(50) norm   // 50个直条
    
11. histogram y,bin(50) norm
    
12. pwcorr x y     //计算相关系数

复制代码

*方法三

1. clear
   
2. set obs 100
   
3. 
   
4. generate x1 = rnormal(0,1)
   
5. generate x2 = rnormal(0.5,1.3)
   
6. 
   
7. summarize x1
   
8. local M1 = r(mean)
   
9. local SD1 = r(sd)
   
10. 
    
11. summarize x2
    
12. local M2 = r(mean)
    
13. local SD2 = r(sd)
    
14. 
    
15. twoway function x1 = normalden(x, `M1', `SD1') , range(-4 4) || ///
    
16. function x2 = normalden(x, `M2', `SD2'), range(-4 5)               

复制代码

扫码加我 拉你入群

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

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

分享0 收藏8 回帖

关键词：标准正态分布 Stata 正态分布 tata histogram 正态分布

顶一个赞

顶一个，谢谢分享。多发此类好贴。

很优秀的帖子，顶一个！！！！

方法三中 x3=rnotrmal（）括号中的0.5和1.3是怎么来的？