【问题标题】:Simulate birth with Rstudio - draw pmf and cdf使用 Rstudio 模拟出生 - 绘制 pmf 和 cdf
【发布时间】:2015-01-07 11:17:35
【问题描述】:

我有一个小练习要使用 Rstudio 来解决我的统计考试。 我试着把它翻译成英文,所以如果有不清楚的地方请向我解释。

“模拟 100,000 次出生并使用以下概率:男性 51.3%,女性 48.7%,使用 sample 函数。

  • 检查获得的男性和女性数量与理论百分比的差异。

  • 画出本实验的概率函数的 PMF 和 CDF(以 50 次出生为样本)。

  • 计算分布的均值和方差。"

我得到了51356个男性和48644个女性,相差56。

但是现在,我怎样才能画出概率函数的PMF和CDF呢?

这里我放了用来模拟分娩的代码:

mysample <- data.frame(sample(c("M","F"),100000,replace=T,prob=c(0.513,0.487)))
names(mysample)<-c("Gender")
males <- subset(mysample, Gender=="M")
females <- subset(mysample,Gender=="F")

theoricM <- 100000*0.513
theoricF <- 100000*0.487
realM <- as.integer(nrow(maschi))
realF <- as.integer(nrow(femmine))

#create a data frame to show differences
result <-data.frame(realM,theoricM,realF,theoricF)
names(result)<- c("Males","Theoric Males","Females","Theoric Females")

结果:

希望有人可以帮助我,我知道对于有 R 经验的人来说这是一个非常简单的问题,但我刚开始使用这种语言。

感谢所有愿意回复的人。

编辑:

我试过这段代码:

x <- 1:50
plot(x,dbinom(x ,size = 50,prob = 0.513),type="l", ylab="PMF", main="Binomial Distribution PMF")

结果是:

我想我的理解是,由于概率非常接近 1/2,在一组 50 名新生儿中,男性人数将非常接近 25。这是什么情节?而且,这是正确的做法吗?

【问题讨论】:

  • 提示:?table, ?cumsum
  • 对不起,我听不懂,你能解释得更好吗?
  • 我不愿意在考试中给你太多帮助。实际上,在重新阅读问题时,我认为您应该查看 ?dbinom(您被要求绘制 理论 采样分布,而不是从大量模拟中观察到的分布)。
  • 这不是考试(当然!),这是为考试做准备的练习;所以我必须学习如何解决它才能通过考试。但是,我已经直觉它是二项分布,我无法理解如何创建理论上的 pmf;例如在这段代码中:如何设置可能的选择只有两个(M 和 F)x &lt;- 1:50 plot(x,dbinom(x ,size = 50,prob = c(0.513,0.487)),type="l", ylab="PMF", main="Binomial Distribution PMF") 我尝试了类似 data &lt;- sample(c("M","F"),50,replace=T,prob=c(0.513,0.487)) 而不是 x,但它不是数字
  • 您应该将这些内容添加到您的问题中,以便更清楚您尝试了什么。 (你已经接近了。)提示#2:对于二项式,您不需要指定两个概率 - 只需指定您算作“成功”的任何结果的概率。 (对于具有 >2 个类别的多项分布,问题会更加困难......)

标签: r cdf


【解决方案1】:

您的代码(和结论)在我看来是正确的。

使用type="h" 绘制“高密度”图可能在图形上更好;这更清楚地表明x 的非整数值的概率为零。

x <- 1:50
par(las=1,bty="l") ## cosmetic
plot(x,dbinom(x ,size = 50,prob = 0.513),type="h", ylab="PMF", 
     main="Binomial Distribution PMF")

(绘制 CDF/CMF 时,您可能需要使用 type="s"type="S";请参阅 ?plot

【讨论】:

  • 好的,非常感谢。我会记住你的建议。
猜你喜欢
  • 1970-01-01
  • 1970-01-01
  • 2011-06-01
  • 2014-10-06
  • 1970-01-01
  • 2013-09-09
  • 2013-10-23
  • 1970-01-01
  • 2010-12-25
相关资源
最近更新 更多