R语言求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

题目是这样的：从正态总体N(555,6.3^2)中随机抽取一个容量为36的样本，求样本的均值落在（53.8，56.8）的概率。直接计算可以得到概率为0.83
可用R模拟计算概率却差的很多，不知道这个模拟过程有什么问题还是代码有问题，代码：
p <- function(n){
     M <- rnorm(10000,55,6.3^2)
     x <- matrix(0,n,36)
     mx <- c(1:n)
     tj <- c(1:n)
     for(i in 1:n){
         x[i,] <- sample(M,36)
         mx <- mean(x[i,]) }
     tj[mx>53.8 & mx <56.8] <- 1
     tj[mx<53.8 | mx >56.8] <- 0
     p <- sum(tj)/n
     return (p)     }

结果：
> p(1000)
[1] 0.186

求高手看看有什么问题或者换种方法解决这道题，谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R语言 不知道 function matrix 样本

我去，被楼主绕晕了，rnorm(n, mean = 0, sd = 1)，rnorm中第三个参数是标准差不是方差，所以应该是rnorm(10000,55,6.3)

大概0.8多一点，没有问题

千5320 发表于 2016-6-20 18:52
我去，被楼主绕晕了，rnorm(n, mean = 0, sd = 1)，rnorm中第三个参数是标准差不是方差，所以应该是rnorm(1 ...

原来如此，困扰了一下午，谢谢！