[R-es] Regresión lineal para una serie temporal
Paco, Ya que tienes una serie temporal (temperaturas que evolucionan a lo largo del tiempo), y parece que tu objetivo es estudiar esta evolución (regresión respecto al tiempo). Resulta que es precisamente el objeto básico de la teorÃa de las series temporales. Es decir que puedes utilizar los modelos clásicos de este área: modelos ArMa, Arima, .... y puedes también incluir una tendencia y un componente estacional. Echa, por ejemplo, un vistazo al paquete "tseries". Un saludo. Olivier -- ____________________________________ Olivier G. Nuñez Email: onunez en iberstat.es Tel : +34 663 03 69 09 Web: http://www.iberstat.es ____________________________________ El 17/02/2011, a las 12:42, Paco Pastor escribió:
Hola Jorge y Óscar
Efectivamente, quiero hacer la regresión de la variable respecto al
tiempo.
He probado vuestras dos sugerencias y obtengo la regresión pero
resulta que dan valores ligeramente diferentes. ¿Qué método debo
escoger? supongo que dependerá del tipo de datos con los que se
trabaja. Además, ésta es solamente una primera aproximación a la
tendencia de la serie y luego tendrá que venir algo más complejo.
library(zoo)
temperatura=read.csv("temp_maxima_C.dat",header=F,sep=",")
tdat=temperatura[2:8]
mitja=rowMeans(tdat, na.rm=TRUE)
temp2=cbind(temperatura,mitja)
tfecha=as.Date(temp2[,1],"%d/%m/%Y")
tdata=zoo(temp2$mitja,tfecha)
lm1=lm(coredata(tdata)~index(tdata))
lm1
Call:
lm(formula = coredata(tdata) ~ index(tdata))
Coefficients:
(Intercept) index(tdata)
2.136e+01 1.119e-04
lm2=lm(tdata~tfecha)
lm2
Call:
lm(formula = tdata ~ tfecha)
Coefficients:
(Intercept) tfecha
2.155e+01 8.312e-05
Gracias
El 17/02/2011 11:56, Jorge Virto escribió:
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
-- ----------- Francisco Pastor Meteorology department, Instituto Universitario CEAM-UMH http://www.ceam.es ----------- mail: paco en ceam.es skype: paco.pastor.guzman Researcher ID: http://www.researcherid.com/rid/B-8331-2008 Cosis profile: http://www.cosis.net/profile/francisco.pastor ----------- Parque Tecnologico, C/ Charles R. Darwin, 14 46980 PATERNA (Valencia), Spain Tlf. 96 131 82 27 - Fax. 96 131 81 90
_______________________________________________ R-help-es mailing list R-help-es en r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es