全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
1317 7
2013-07-08
问题:生成一个数据集,该数据集中大约30%是服从N(450,150),其余部分服从N(300,150),请问如何操作?


二维码

扫码加我 拉你入群

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

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

全部回复
2013-7-8 23:27:35
**原本想说以下的话
楼主您应当对于您的服从N(450,150)进一步的说明,特别是后面的那个数字,
因为在 R 与 Stata 的习惯是不一样的。

**我发现现在的人都不喜欢讲清楚,好像问的人是高手,回答的人也是高手

**我的演练,楼主高兴就参考,可以当废物的,没关系

set se 888
set ob 1000
g x=0.3*rnormal(450,150)+0.7*rnormal(300,150)

hist x,norm

g x1=rnormal(450,150) in 1/300
g x2=rnormal(300,150) in 301/1000
replace x1=x2 if x1==.
drop x2
g y=(_n<=300)

hist x1,norm
tw hist x1,by(y, total)
二维码

扫码加我 拉你入群

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

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

2013-7-8 23:39:45
h3327156 发表于 2013-7-8 23:27
**原本想说以下的话
楼主您应当对于您的服从N(450,150)进一步的说明,特别是后面的那个数字,
因为在 R ...
=_=

蓝色版主, 我只是想表达 "提问的人要讲清楚些啦! 可以减少想回答的人去猜"

所以我演练的回答,也就不废话了! 换提问的人去猜吧!
二维码

扫码加我 拉你入群

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

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

2013-7-9 00:21:54

猜的很累啊!

为什么有些人就是不能写的清楚明了一些呢?
二维码

扫码加我 拉你入群

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

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

2013-7-9 00:36:03
h3327156 发表于 2013-7-8 23:27
**原本想说以下的话
楼主您应当对于您的服从N(450,150)进一步的说明,特别是后面的那个数字,
因为在 R ...
tw (hist x1 if y==0, bin(50) fcolor(green))  ///
   (hist x1 if y==1, bin(50) fcolor(none) lcolor(blue))   ///
   , legend(label(1 "N(450,150)") label(2 "N(300,150)"))

在你的基础上,添加一句,这样看容易些。
幸亏有你生成的y这个标示变量,否则混在一起,看不出背后的数据生成过程了
Graph.png
二维码

扫码加我 拉你入群

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

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

2013-7-9 01:30:07
h3327156 发表于 2013-7-8 23:27
**原本想说以下的话
楼主您应当对于您的服从N(450,150)进一步的说明,特别是后面的那个数字,
因为在 R ...
恩,谢谢谢谢!
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群