求帮忙解决随机数的生成问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

已知\[ e ~ 0.5*N(-1,2.5^2)+0.5*N(1,0.5^2)\]，在R语言中如何产生e 的随机数

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：求帮忙 随机数 LaTeX atex Late 如何

e ~ 0.5N(-1,2.5^2)+0.5N(1,0.5^2)

1. random_num <- function(n,mean1,sigma1,mean2,sigma2){
   
2.   num <- numeric(n)
   
3.   for(i in seq_along(num)){
   
4.     num[i] <- 0.5*rnorm(1,mean1,sigma1) + 0.5*rnorm(1,mean2,sigma2)
   
5.   }
   
6.   num
   
7. }

复制代码

然后你的那个分布随机数可以用

1. random_num(100,-1,2.5,1,0.5)  #这里生成了100个随机数

复制代码

houyunhuang 发表于 2016-10-10 16:35
然后你的那个分布随机数可以用

我看的论文e服从上面那个分布，论文上面讲了E(e)=0，Median(e)=0.67，Mode(e)=1,其中Mode表示众数。
可是按你的这种方式产生的随机数中位数和众数都不能近似的等于它们。

shizhangbo 发表于 2016-10-10 17:14
我看的论文e服从上面那个分布，论文上面讲了E(e)=0，Median(e)=0.67，Mode(e)=1,其中Mode表示众数。
可是 ...

我尝试了生成10万个随机数，均值为0.00220842，中位数是0.003366811，密度最大处为0.05444197。可能这样生成随机数有问题

正态分布的线性组合不是该服从正态分布？

shizhangbo 发表于 2016-10-10 17:14
我看的论文e服从上面那个分布，论文上面讲了E(e)=0，Median(e)=0.67，Mode(e)=1,其中Mode表示众数。
可是 ...

好象，论文是错的！！

你好,我通过问题看到你在阅读modal模型的论文,可以请教你点问题并交流一下吗?

你好,我通过问题看到你在阅读modal模型的论文,可以请教你点问题并交流一下吗?

你好,我通过问题看到你在阅读modal模型的论文,可以请教你点问题并交流一下吗?