微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何在没有for循环的情况下重复生成随机样本1000次并绘制alpha估计量的图?

如何解决如何在没有for循环的情况下重复生成随机样本1000次并绘制alpha估计量的图?

我必须根据beta分布(alpha = 5,beta = 1)生成一个n = 20的随机样本,并计算alpha的MLE。我使用功能rbeta和ebeta做到了这一点:

alfa = 5
beta=1
n = 20
sample<- rbeta(n,alfa,beta)
mle_alfa <- ebeta(proba,method = "mle")$parameters[1]

但是我必须重复指向它1000次(即生成1000个样本)并绘制直方图,箱形图和q-qplots

enter image description here

我试图像这样在循环中重复它,但是我不知道它是否正确以及如何使用它来绘制曲线。你能帮我么?谢谢!

sample2<- sapply(1:1000,function(x) rbeta(n,beta))
sample2<- as.numeric(sample2)

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。