En particular lo resolví de la siguiente manera:

Mi <- 134
sigma <- 37.841
n <- 30
sims <- 100
cont <- rep(0, sims)
alfa <- 0.05
*res=matrix(0,sims,2)*
for(i in 1:sims) {
   x <- rnorm(n, Mi, sigma)
   Lim.inf <- (mean(x) - (qt(1-alfa  / 2, n-1) * (sd(x) / sqrt(n))))
   Lim.sup <- (mean(x) + (qt(1-alfa  / 2, n-1) * (sd(x) / sqrt(n))))
   if(Lim.inf < Mi && Mi < Lim.sup)
   {cont[i] <- 1}
*res[i,1]=Lim.inf
   res[i,2]=Lim.sup*
   }
sum(cont) / sims

Le agregue una matriz que recogerá los resultados y por cada iteración 
ingresa un resultado en la línea correspondiente. De esa forma salva 
todos los valores, si desearas incluir solamente los que cumplan la 
condición que colocas deberías incluirlos en el if, creo.

Puede haber otra opción más eficiente, de repente otros con mayor 
experiencia lo pueden comentar mejor.

Saludos

Fernando Macedo


El 08/09/12 23:12, Lex B escribió:
> Buenas noches, las sentencias de repeticion me han generado muchos
> problemas, en este caso, la idea es almacenar como  filas de un data.frame
> los resultados de las funciones Lim.inf  y Lim.sup, les agradezco su ayuda,
> el codigo que quiero correr es el siguiente, y lo hace bastante bien, pero
> no puedo ni visualizar y mucho menos trabajar sobre los resultados de estas
> dos funciones.
>
>
> Mi <- 134
> sigma <- 37.841
> n <- 30
> sims <- 100
> cont <- rep(0, sims)
> alfa <- 0.05
> for(i in 1:sims) {
> x <- rnorm(n, Mi, sigma)
> Lim.inf <- (mean(x) - (qt(1-alfa  / 2, n-1) * (sd(x) / sqrt(n))))
> Lim.sup <- (mean(x) + (qt(1-alfa  / 2, n-1) * (sd(x) / sqrt(n))))
> if(Lim.inf < Mi && Mi < Lim.sup)
> {cont[i] <- 1}
> }
> sum(cont) / sims
>
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


	[[alternative HTML version deleted]]

