[R-es] Regresión lineal para una serie temporal
Oscar Perpiñan Lamigueiro
oscar.perpinan en upm.es
Jue Feb 17 12:10:28 CET 2011
El Thu, 17 Feb 2011 11:38:40 +0100
Paco Pastor <paco en ceam.es> 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
>
OK.
Entonces el objeto zoo lo creas así:
dat <- read.zoo(file='dat.csv', sep=',', dec='.',
format='%Y-%m-%d', header=TRUE)
Y cuando hablas de regresión lineal con una sola variable, ¿te refieres
a ajustar frente al tiempo?
lmZ <- lm(coredata(dat)~index(dat))
summary(lmZ)
¿o necesitas algo más elaborado?:
dat.fit <- StructTS(dat, type='level')
plot(dat)
lines(fitted(dat.fit), col='green')
Saludos.
Oscar.
-------------------
Oscar Perpiñán Lamigueiro
Dpto. Ingeniería Eléctrica
EUITI-UPM
Más información sobre la lista de distribución R-help-es