Центральная предельная теорема в [R] - PullRequest
Новая
       50

Центральная предельная теорема в [R]

0 голосов
/ 30 января 2012

Я работаю с языком [R], чтобы сгенерировать выборку из M = 32000 средних значений, каждое из которых рассчитывается путем усреднения 36 независимых значений случайного переменного непрерывного равномерного распределения (0, 1), генерируемого следующим образом: 

sampleA<-1:32000

for ( i in 1:32000){
    MuestraAUnif<- runif(36)
    sampleA[i]<-mean(MuestraAUnif)
}

Для сгенерированной выборки попросите меня вычислить относительную частоту наблюдаемых средних, превышающих L = 0,32 +4 * 1/100, и сравните ее с вероятностью (аппроксимируемой « Центральная предельная теорема »)что средние значения N больше, чем L. следующим образом: 

    L<- 0.32+4*1/100
    sigma<- sqrt(1/12) #(b-a)/12 
    miu = 0.5 #(a+b)/2
    greaterA <-sum(sampleA > L) #values of the sample greater than L are 23693
    xBar<- greaterA/length(sampleA) 
    X <- sum(sampleA) 
    n<-32000
    Zn<- (X - n*miu)/(sigma*sqrt(n))

    cat("P(xBar >",L,") = P(Z>", Zn, ")=","1 - P (Z < ", Zn,") =",1-pnorm(Zn),"\n") #print the theoretical prob Xbar greater than L
    cat("sum (sampleA >",L,")/","M=", n," para N =", 36,":",xBar, "\n") #print the sampling probability print when is greater than L

Вывод: 

P(xBar > 0.36 ) = P(Z> -3.961838 )= 1 - P (Z <  -3.961838 ) = 0.9999628 
sum (sampleA > 0.36 )/ M= 32000  para N = 36 : 0.7377187

Мой вопрос таков: почему до сих пор значения?они должны быть намного ближе (0.9999628 далеко от 0.7377187).Что-то не так с моей реализацией?Извините за мой английский.

1 Ответ

2 голосов
/ 30 января 2012

Melkhiah66. Ты все правильно сделал только поменяй MuestraAUnif<- runif(2) для MuestraAUnif<- runif(32) и это должно работать

...