模拟R中的二项分布

Question

我试图在p中模拟R中的二项式分布T~B（10，p） p~U（0,1）。

我创建了以下代码：

 Bin_sim<-function(n) {
Prob<- sample(0:10, size = n, replace=T)
Prob2<-rbinom(n=n, size=10, prob=(Prob/10))
Prob3<-Prob[Prob2==8]

table(Prob3) / length(Prob3)
}

以下代码生成此结果：

 Bin_sim(10000)
Prob3
          3           4           5           6           7           8 
0.001322751 0.009259259 0.039682540 0.141534392 0.268518519 0.320105820 
          9 
0.219576720 

如果我想模拟T = 8的所有成功概率的分布，我是否正朝着正确的方向前进？ 提前感谢，帮助会很有用。