Hola,

En la misma línea de la solución que ya te han proporcionado.
Incluyo el código añadido en rojo además de un breve comentario con las
instrucciones básicas a replicar en casos similares.

####################################
Mi <- 134
sigma <- 37.841
n <- 30
sims <- 100
cont <- rep(0, sims)
alfa <- 0.05

#Inicializar y dimensionar el data.frame donde se almacenarán los
resultados.
#El dimensionar no es estrictamente necesario, pero si se quieren modificar
data.frames grandes ayuda.
resul.df <- data.frame(lim.inf=rep(0,sims), lim.sup=rep(0,sims))

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))))

#Incluyo en cada iteración el valor de Lim.inf y Lim.sup en una nueva fila
del data.frame.
  resul.df$lim.inf[i] <- Lim.inf
  resul.df$lim.sup[i] <- Lim.sup

  if(Lim.inf < Mi && Mi < Lim.sup)
  {cont[i] <- 1}
}
sum(cont) / sims

####################################

Saludos,
Carlos Ortega
www.qualityexcellence.es

El 9 de septiembre de 2012 04:12, Lex B <bernalm.alexander@gmail.com>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
>
>
> --
>
> ALEXANDER BERNAL MUÑOZ
> U SANTO TOMAS - ESTADÍSTICA - BOGOTÁ
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

