Имитационный тест для испытаний Бернулли для увеличения размера выборки в R - PullRequest
2 голосов
/ 29 мая 2020

В программном обеспечении R я пытаюсь создать график моделирования для испытания Бернулли для различных размеров выборки, как показано ниже Simulation of Bernoulli Trials for increasing sample size

Я создал пропорции для данного размера выборки и повторил процесс 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

1 Ответ

2 голосов
/ 29 мая 2020

Я думаю, что вам нужно это.

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)

enter image description here

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...