Patricio Fuenmayor Viteri
2011-Feb-17 13:42 UTC
[R-es] Resumen de R-help-es, Vol 24, Envío 17
Hola Paco. Esto de las series temporales es un tema muy extenso y aveces apasionante.En R hay varias librarías que tiene la finalidad de tratarlas.Es correcto que trates de realizar regresiones en tu serie, pero debes especificar cuales son tus regresores, esa en la diferencia fundamental de este tipo de datos.Para modelarlas se asume que la serie puede explicarse por si misma, es decir no necesita otras variables, esto en el caso univariante... si es multivariante ya es otra cosa.Te recomiendo que revises estos documentos y estas librerías, y vas a tener mas clara la película y tendrás una idea de las ventajas que R te brinda. http://cran.r-project.org/web/views/TimeSeries.html http://cran.r-project.org/doc/contrib/Farnsworth-EconometricsInR.pdf http://cran.r-project.org/doc/contrib/Ricci-refcard-ts.pdf en base a esto y por experiencia te recomiendo que primero realices gráficas descriptivas, para ver que características puedes ver en tu serie.Aplica cualquiera de las técnicas que se tiene en la libraría forecast, y analiza los resultados y realiza pronósticos.. que en resumen es una de las finalidades de esto, y en base a los resultados, selecciona el modelo.También puedes intentar las técnicas de descomposición, que da una idea mas general, si tu serie es estacional.Yo uso la función auto.arima de la librería forecast, ayuda mucho para la selección automática del mejor modelo, pero ojo, no necesariamente el mejor modelo estadísticamente hablando, es el mas adecuado para explicar tu fenómeno. Saludos.Patricio> From: r-help-es-request@r-project.org > Subject: Resumen de R-help-es, Vol 24, Envío 17 > To: r-help-es@r-project.org > Date: Thu, 17 Feb 2011 12:00:06 +0100 > > Envíe los mensajes para la lista R-help-es a > r-help-es@r-project.org > > Para subscribirse o anular su subscripción a través de la WEB > https://stat.ethz.ch/mailman/listinfo/r-help-es > > O por correo electrónico, enviando un mensaje con el texto "help" en > el asunto (subject) o en el cuerpo a: > r-help-es-request@r-project.org > > Puede contactar con el responsable de la lista escribiendo a: > r-help-es-owner@r-project.org > > Si responde a algún contenido de este mensaje, por favor, edite la > linea del asunto (subject) para que el texto sea mas especifico que: > "Re: Contents of R-help-es digest...". Además, por favor, incluya en > la respuesta sólo aquellas partes del mensaje a las que está > respondiendo. > > > Asuntos del día: > > 1. Re: Regresión lineal para una serie temporal (Jorge Virto) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 17 Feb 2011 11:56:43 +0100 > From: Jorge Virto <jorge.virto@ehu.es> > To: r-help-es@r-project.org > Subject: Re: [R-es] Regresión lineal para una serie temporal > Message-ID: <201102171156.43956.jorge.virto@ehu.es> > Content-Type: text/plain; charset="utf-8" > > Hola Paco, > > creo que el problema está en que no le dices a R que tu primera variable es > una fecha. Prueba lo siguiente: > > datos <- read.csv("dat.csv",header=T) > > # Das formato fecha a los datos > # Fecha del tipo: Año-Mes-Día > fecha <- as.Date(datos[,1],format="%Y-%m-%d") > > reg <- lm(datos[,2]~fecha) > summary(reg) > > Un saludo, > > Jorge > > El Thursday 17 February 2011 11:38:40 Paco Pastor escribió: > > Hola > > > > Tienes razón, me falta un ejemplo. Mis datos son de la forma: > > > > > > "date" "temp" > > 01/01/2000 21.6 > > 02/01/2000 22.4 > > ..... .... > > 31/12/2010 22.3 > > > > No tengo dos variables x-y para hacer la regresión sino solamente una > > variable "temp". Leo los datos, creo un objeto zoo para representar la > > serie temporal pero no sé como seguir para hacer la regresión lineal. > > > > Adjunto los datos de ejemplo guardados con write.csv > > > > Gracias > > > > El 17/02/2011 11:18, Oscar Perpiñan Lamigueiro escribió: > > > El Thu, 17 Feb 2011 10:52:42 +0100 > > > > > > Paco Pastor<paco@ceam.es> escribió: > > >> Hola a todos > > >> > > >> Estoy intentando hacer un análisis rápido de una serie temporal de > > >> datos diarios pero me encuentro con algunos problemas. Me gustaría en > > >> primera instancia hacer una regresión lineal pero no encuentro la > > >> forma. > > >> > > >> Tras leer los datos diarios creo un objeto de la clase zoo y sobre > > >> éste no puedo utilizar lm(). He leído algo sobre dynlm pero no > > >> encuentro la forma. > > >> > > >> Se agradece cualquier sugerencia, > > >> > > >> Paco > > > > > > A falta de un ejemplo, me invento uno con lo que he entendido de tu > > > mensaje: > > > > > > library(zoo) > > > library(lattice) > > > > > > ##Dos variables relacionadas entre sí > > > x=rnorm(100) > > > y=2*x+1 > > > > > > ##Indice temporal > > > idx=seq(as.POSIXct(''2011-01-01''), length=100, by=''day'') > > > > > > ##Creo el objeto zoo y lo pinto > > > z<- zoo(data.frame(x, y), idx) > > > xyplot(z, superpose=TRUE) > > > > > > ##Hago un ajuste entre las variables > > > lmZ<- lm(y~x, z) > > > summary(lmZ) > > > > > > Saludos. > > > > > > Oscar. > > > > > > ------------------- > > > Oscar Perpiñán Lamigueiro > > > Dpto. Ingeniería Eléctrica > > > EUITI-UPM > > > > ------------------------------ > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > > > Fin de Resumen de R-help-es, Vol 24, Envío 17 > *********************************************[[alternative HTML version deleted]]