R语言与概率论

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

当p=0.15，n=5000时，怎么画出各个位置的小球个数的理论上的密度分布曲线？

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：R语言 概率论 概率论

N <- c(0:20)
P <- data.frame( Prob = choose(20, N)* 0.15^(20-N) *0.85^N )
P$N <- round(2000*P$Prob)

#画图：
barplot(P$N,col = "red")

plot(P$Prob, type = "b")

什么意思啊，里面的参数不是很懂啊

死蠢猫 发表于 2016-6-24 17:19
什么意思啊，里面的参数不是很懂啊

加了注释你再看看吧，不知道你到底要哪样的图。

1. N <- c(0:20) #21个格子，0到20，代表会向右边掉多少次
   
2. P <- data.frame( Prob = choose(20, N)* 0.15^(20-N) *0.85^N )
   
3.             #计算掉到每个格子的概率，排列问题
   
4. P$Num <- round(5000*P$Prob)
   
5.             #四舍五入，5000个球掉入每个格子球数的数学期望
   
6. 
   
7. #画图：
   
8. barplot(P$Num,col = "red") #条形图
   
9. plot(P$Prob, type = "b") #点线图

复制代码

谢谢 爱迪生_陈！借花献佛，稍微改造了一下程序。
m = 20  # 钉板层数，也是最底层钉子数
p = 0.15  # 落在钉子左侧的概率
n = 5000  # 玻璃球数
space <- c(0:m)  # 最下面21个格子
# 计算掉到每个格子的概率，排列问题
result <- data.frame(Prob = choose(m, space)* p^(m-space) *(1-p)^ space)
# 四舍五入，5000个球掉入每个格子球数的数学期望
result$Num <- round(n*result$Prob)
            
#画图：
barplot(result$Num, col = "red")  # 条形图
plot(result$Prob, type = "b")     # 点线图