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]]
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]]
Hola el otro dia intente hacer una tecnica multivariante (kmeans) con una matriz de datos que es un spatial.grid.data.frame. Pero R me dice que no existe tal método para un objeto de clase S4. He estado leyendo sobre el asunto de los objetos S3 yS4 pero por desgracia, no entiendo mucho lo que he encontrado. ¿Hay alguien por aqui que pueda decirme por qué no puedo hacer este analisis con este tipo de objetos, si he de transformarlo a S3 o me "fastidiará" alguna propiedad de la matriz de datos, o si hay algún análisis semejante para este tipo de objetos?? Gracias! saludos. Arantzazu Blanco Bernardeau Dpto de Química Agrícola, Geología y Edafología Universidad de Murcia-Campus de Espinardo> From: jorgeivanvelez@gmail.com > Date: Tue, 20 Jul 2010 21:49:49 -0400 > To: andreswk25@hotmail.com > CC: r-help-es@r-project.org > Subject: Re: [R-es] simulaciones Ar(1) > > 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]] > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es_________________________________________________________________ Accede a tu Hotmail en un solo clic ¡Descárgate Internet Explorer 8 y e[[elided Hotmail spam]] [[alternative HTML version deleted]]
Hola, ¿qué tal? De la página de ayuda de kmeans: x numeric matrix of data, or an object that can be coerced to such a matrix (such as a numeric vector or a data frame with all numeric columns) A kmeans le tienes que pasar un objeto como el que describe la ayuda. Un spatial.grid.data.frame no lo es. Pero uno de sus "slots" es el data.frame data. Sólo tienes que extraerlo. En http://www.stat.auckland.ac.nz/S-Workshop/Gentleman/S4Objects.pdf puedes ver cómo extraer datos de los "slots" de un objeto S4 (básicamente, usando @ en lugar de $). Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 27 de julio de 2010 15:28, Arantzazu Blanco Bernardeau <aramucia en hotmail.com> escribió:> > > Hola > el otro dia intente hacer una tecnica multivariante (kmeans) con una matriz de datos que es un spatial.grid.data.frame. Pero R me dice que no existe tal método para un objeto de clase S4. He estado leyendo sobre el asunto de los objetos S3 yS4 pero por desgracia, no entiendo mucho lo que he encontrado. > ¿Hay alguien por aqui que pueda decirme por qué no puedo hacer este analisis con este tipo de objetos, si he de transformarlo a S3 o me "fastidiará" alguna propiedad de la matriz de datos, o si hay algún análisis semejante para este tipo de objetos?? > Gracias! saludos. > Arantzazu Blanco Bernardeau > Dpto de Química Agrícola, Geología y Edafología > Universidad de Murcia-Campus de Espinardo
Hola, Normalmente un objeto S4 tiene unos métodos asociados para extraer información en formatos "convencionales". Además del documento que te recomienda Carlos, te puede resultar útil este otro algo más extenso: http://cran.r-project.org/doc/contrib/Genolini-S4tutorialV0-5en.pdf. En este "tutorial" recomienda emplear estos métodos frente al acceso directo al slot con @. En el caso concreto del SpatialGridDataFrame, hay un método para obtener información en formato "matrix", y se llama as.matrix: http://www.dpi.inpe.br/gilberto/tutorials/software/R-contrib/sp/html/SpatialGridDataFrame-class.html. Con lo que ahí obtengas podrás emplear kmeans. Los autores del paquete sp (el que contiene a la clase en cuestión) escribieron un libro y tienen una página web con información relacionada. Quizás te pueda interesar: http://r-spatial.sourceforge.net/. Y aprovecho el correo para indicarte un paquete llamado cluster (http://cran.r-project.org/web/packages/cluster/) que contiene métodos más robustos para hacer tareas de clasificación (por ejemplo, la función pam). Saludos. Oscar. -- Oscar Perpiñán Lamigueiro Profesor Asociado Dpto. de Ingeniería Eléctrica EUITI-UPM On Tue, 27 Jul 2010 17:27:18 +0200 "Carlos J. Gil Bellosta " <cgb en datanalytics.com> wrote:> Hola, ¿qué tal? > > De la página de ayuda de kmeans: > > x numeric matrix of data, or an object that can be coerced to such a > matrix (such as a numeric vector or a data frame with all numeric > columns) > > A kmeans le tienes que pasar un objeto como el que describe la ayuda. > Un spatial.grid.data.frame no lo es. Pero uno de sus "slots" es el > data.frame data. Sólo tienes que extraerlo. > > En http://www.stat.auckland.ac.nz/S-Workshop/Gentleman/S4Objects.pdf > puedes ver cómo extraer datos de los "slots" de un objeto S4 > (básicamente, usando @ en lugar de $). > > Un saludo, > > Carlos J. Gil Bellosta > http://www.datanalytics.com > > > > El día 27 de julio de 2010 15:28, Arantzazu Blanco Bernardeau > <aramucia en hotmail.com> escribió: > > > > > > Hola > > el otro dia intente hacer una tecnica multivariante (kmeans) con una matriz de datos que es un spatial.grid.data.frame. Pero R me dice que no existe tal método para un objeto de clase S4. He estado leyendo sobre el asunto de los objetos S3 yS4 pero por desgracia, no entiendo mucho lo que he encontrado. > > ¿Hay alguien por aqui que pueda decirme por qué no puedo hacer este analisis con este tipo de objetos, si he de transformarlo a S3 o me "fastidiará" alguna propiedad de la matriz de datos, o si hay algún análisis semejante para este tipo de objetos?? > > Gracias! saludos. > > Arantzazu Blanco Bernardeau > > Dpto de Química Agrícola, Geología y Edafología > > Universidad de Murcia-Campus de Espinardo > > _______________________________________________ > R-help-es mailing list > R-help-es en r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es