Hola,

Si lo he entendido bien, te están pidiendo el intervalo de confianza para
cada "y" de cada iteración.
Y que para cada una de esas "y" veas si su intervalo de confianza calculado
captura el intervalo de confianza teórico.

Lo que has hecho es calcular el intervalo de confianza de las 1000 "y"
juntas.

Yo lo haría así:


test.inf<-0
test.sup<-0
 for(i in 1:1000)
 {
 y<-rnorm(10, 2.5, 0.16)
 test.inf[i]<-t.test(y, conf.level=0.95)$conf.int[1]
 test.sup[i]<-t.test(y, conf.level=0.95)$conf.int[2]
 }
test.df<-data.frame(inf=test.inf, sup=test.sup, dif.inf=test.inf-(2.5-0.16)
,dif.sup=test.sup-(2.5+0.16))
good.df<-test.df[test.df$dif.inf>0 & test.df$dif.sup>0, ]

por.val<-dim(good.df)[1]*100/1000
print(c("Porcentaje Correcto: ", por.val))

También comentarte que además de un mensaje desesperado en el "Subject"
ayuda el incluir un poco más de detalle del tipo de problema que quieres
resolver...

Carlos Ortega
www.qualityexcellence.es
http://qualityexcellence.wordpress.com/ (*Nuevo Blog*)




2011/4/16 Angela Cadena <angelacadenaramirez@gmail.com>

> hola!
>
> podrian ayudarme con esto:
>
>
> Genere 1000 muestras aleatorias de tamaño  n=10 de una población
> normal con media
>  =2,5 y desviación  0,16. Con cada una de las muestras construya un
> intervalo de
> confianza al 95% para la media y para la varianza. Cuente que
> porcentaje de los 1000
> intervalos atrapan la verdadera media y la verdadera varianza poblacional
> respectivamente. Comente los resultados encontrados.
>
> he hecho esto
>
> promedio=0
> varianza=0
>  for(i in 1:1000)
>  {
>  y<-rnorm(10, 2.5, 0.16)
>  promedio[i]=mean(y)
>  varianza[i]=var(y)
>  }
> promedio
> varianza
> #para generar 1000 prom y 1000 varianzas para la muestra de n=10,
> mu=2.5 y dev=0.16
>
> intmedia=c(promedio)            #intmedia =iintervalo de confianza para la
> media
> t.test(intmedia, conf.level=0.95)
>
>       One Sample t-test
>
> data:  intmedia
> t = 1575.343, df = 999, p-value < 2.2e-16
> alternative hypothesis: true mean is not equal to 0
> 95 percent confidence interval:
>  2.494522 2.500744
> sample estimates:
> mean of x
>  2.497633
>
> # el intervalo de conf para la media seria ( 2.494522 , 2.500744 )
>
>
>  ¿para la varianza seria igual?
>  ¿como hago para contar los intervalos que atrapan las verderas medias
> y las verdaderas varianzas?
>
> soy novata en esto, por favor ayudenme!
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

	[[alternative HTML version deleted]]

