[R-es] Cosinor

[Ana Pérez V.]

Hola, ¿qué tal?:

Tengo la siguiente duda:

Para el caso del Cosinor Simple:

Y = M +  A * cos(2*pi*t/T + phi) + error

Cuando T conocido, podemos asumir:

 Y = M + beta * X1 + gamma * X2 + error, 

siendo X1 = cos(2*pi*t/T) y X2 = sin(2*pi*t/T)

En este caso (Cosinor simple) el modelado lo llevo a cabo del siguiente modo: 

lm( y ~ X1 + X2)

Para el caso de componentes múltiples:

Y = M + ∑Aj * cos(2*pi*t/Tj + phij) + error

Si Tj conocidos: Y = M + ∑(betaj*X1j + gammaj * X2j) + error

Modelado:  y ~ x1tot+ x2tot, donde los xtot (x totales) resultan de la suma de los parciales (donde cada parcial corresponde a un Tj).

Sin embargo, en el caso del modelo generalizado (componente cuadrática, cúbica, etc.):

Y = M + alfa1*t + alfa2*t2 +  (∑Aj * cos(2*pi*t/Tj + phij)) + error

¿Qué facilidades ofrece R para el estudio de este modelo?. Se hace manualmente?.

Gracias.

Un saludo.