[R-es] Regresión lineal para una serie temporal

[Oscar Perpiñan Lamigueiro]

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