r - 伯努利试验的模拟测试，用于增加 R 中的样本量

Question

在 R 软件中，我正在尝试为不同样本大小的伯努利试验生成模拟图，如下所示

我已经为给定的样本量生成了比例，并重复了这个过程 1000 次。但我不明白，如何将图形绘制为给定的图形。

m=1000 ###### replications
n=10  ###### sample size
my=function(){
n=10 ######Trials
u=rbinom(n,1,0.5)
}
Simulated_data=t((replicate(m, my())))
Proportions<-apply(Simulated_data,1,sum)/n
Proportions

Answer 1

我想你想要的是这个。

R <- 50
n <- 1000

set.seed(42)
res <- replicate(R, rbinom(n, 1, .5))
p <- apply(res, 2, function(x) cumsum(x)/seq_along(x))
plot(1:1000, ylim=c(0, 1), type="n")
apply(p, 2, function(x) lines(1:1000, x, col="grey50"))
abline(h=.5)