Hola:
 
Tengo un modelo de regresión lineal en el cual las componentes son cosenoidales, y lo construyo del siguiente modo:
 
 
modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]"
x1[, 1] = cos(2 * pi * t / periods[1])
x2[, 1] = sin(2 * pi * t / periods[1])
for (i in 2:nComp) {
  x1[, i] = cos(2 * pi * t / periods[i])
  x2[, i] = sin(2 * pi * t / periods[i])
  modelo = paste(modelo, " + x1[,", i, "] + x2[,", i, "]", sep="")
}
res = lm (as.formula(modelo))
 
Este modelo deriva en resultados correctos de su análisis correspondiente (summary(res)). Sin embargo, si añado regresión polinómica, el summary no me devuelve los estimadores correctos para las componentes lineal, cuadrática y cúbica.
 
La línea del modelo la modifico del siguiente modo:
 
modelo = "y ~ I(t) + I(t^2) + I(t^3) + x1[, 1] + x2[, 1]"
 
He probado también así:
 
modelo = "y ~  t +  t^2 + t^3 + x1[, 1] + x2[, 1]"
 
¡Pero no hay manera!
 
¿Me podéis orientar?.
 
Gracias
 
Eva
	[[alternative HTML version deleted]]

