模拟分布

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

问题：生成一个数据集，该数据集中大约30%是服从N（450,150），其余部分服从N（300,150），请问如何操作？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：如何操作 数据集 如何

**原本想说以下的话
楼主您应当对于您的服从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)

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

=_=

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

所以我演练的回答，也就不废话了！ 换提问的人去猜吧！

猜的很累啊！

为什么有些人就是不能写的清楚明了一些呢？

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这个标示变量，否则混在一起，看不出背后的数据生成过程了

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

恩，谢谢谢谢！

蓝色 发表于 2013-7-9 00:36
tw (hist x1 if y==0, bin(50) fcolor(green))  ///
   (hist x1 if y==1, bin(50) fcolor(none) lcolor ...

非常感谢！提问确是一门艺术啊！

学习了