Как найти R-приближение для CLT

Question

Распределение по количеству неисправных лампочек является биномиальным распределением. то есть количество неисправных лампочек:

Используйте R для создания графического представления, которое сравнивает истинные значения pmf для числа найденных неисправных лампочек и их приближение с использованием центральной предельной теоремы.

Итак, я создал график для pmf в r, используя приведенный ниже код. Однако у меня возникла проблема с CLT-частью вопроса. Мы еще не рассмотрели эту часть учебного плана.

n <- 100
p <- 0.015
# pmf
x <- (0:n)
px <- dbinom(x,n,p) 
# Plot pmf
plot(x,px,type="h",ylab="p(x)" ,main=paste("Binomial(",n,",",p,")",sep=""))

Answer 1

Я не уверен, что это именно то, что вы хотите для графика CLT:

xclt <- seq(0,n,by = 0.001)
clt <- dnorm(xclt,mean = n*p,sd = sqrt(n*p*(1-p)))
lines(xclt,clt,col="red",lty= 2)
legend(80,0.3,legend = c("bin","clt"),col = c("black","red"),lty = c(1,2))

такой, что