Buenas noches Andres,

Una manera es usando lapply():

# simulacion

mis.phis <- seq(0.15,0.90,0.05)

simula <- lapply(mis.phis, function(PHI)

      arima.sim(list(order = c(1,0,0), ar = PHI), n = 1000, sd = 0.5)

      )

names(simula) <- mis.phis


 # media para cada phi

lapply(simula, mean)


Saludos,

Jorge Ivan Velez



2010/7/20 Andres rodriguez trujillo <>

>
> Cordial Saludo
>
> Ando trabajando con el comando arima.sim de la libreria stats que por medio
> de ella me realiza simulaciones de procesos ARIMA, ARMA, AR o MA, y pues
> necesito simular cierto numero de procesos AR(1) en realidad 1000 variando
> el valor del phi del modelo AR(1), simular un proceso AR(1) es de la sgte
> manera,
>
> arima.sim(list(order = c(1,0,0), ar = 0.15), n = 1000,sd=0.5)  donde el
> argumento "ar" es el phi, esto es una simulacion de un proceso Ar(1) con
> phi=0.15
>
> y pues lo que en realidad necesito es hacer el mismo numero de simulaciones
> pero con diferentes phis  0.15, 0.20, 0.25, 0.30, 0.35,..., 0.90
> y el algoritmo que he realizado no me arroja lo que necesito, pues creo y
> asumo que me esta tomando el primer valor de la simulacion variando el valor
> de phi y esta mal, he aqui el algoritmo
>
> l<-1
> d<-array(0,1000)
> for (i in seq(0.15,0.90,0.05)){
> d[l]<-arima.sim(list(order = c(1,0,0), ar = i), n = 1000,sd=0.5)
> l<-l+1
> }
> d
>
> Hubo 16 avisos (use warnings() para verlos)
>
> luego vienen las 16 primeras simulaciones que son el numero total de los
> phis y despues Ceros
>
> por utlimo es posible calcular la media de cada una de las simulaciones
> para cada phi
>
> Agradesco la ayuda y muchas Gracias.
>
> Andres.
>
>
> _________________________________________________________________
>
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

	[[alternative HTML version deleted]]

