Me estoy convirtiendo un alto dimensiones modelo a un menor modelo tridimensional mediante la instalación de una ventana deslizante de un lineal (paramétrico) modelo y examinando la evolución de los valores de los parámetros a lo largo del tiempo. Me voy de 6.3 millones de puntos alrededor de 2500 valores de 6 parámetros.
La física dice que la intersección debe ser un valor constante, pero cuando uso LM se mueve alrededor. Creo que el movimiento es debido al ruido, y que es la causa de otros valores del parámetro de no indicar correctamente. Me gustaría establecer un valor constante.
Cómo hacer un modelo lineal en R que se ha prescrito para interceptar (no cero).
Código actual:
for (i in 1:(n-k)){
fit <- lm(y ~ x1 + x2 + I(x3^2) + x4 + x5 + x6 , data=data[i:(i+k),])
#STORE PARAMETERS INTO VARIABLES
... #truncated for brevity
}
Código que no hace el trabajo:
fit <- lm(y ~ I(9.81) + x1 + x2 + I(x3^2) + x4 + x5 + x6 , data=data[i:(i+k),])
fit <- lm(y ~ 9.81 + x1 + x2 + I(x3^2) + x4 + x5 + x6 , data=data[i:(i+k),])
Pregunta:
- ¿Cómo puedo recetar el constante?
- He intentado buscar por esto tanto en google y la CV - hay un vocabulario que me estoy perdiendo?
- Puede usted comentar sobre cómo algo como AIC o R2 son afectados por este modelo? Yo prefiero usar AIC o BIC y creo que, como modelo de los criterios de selección se deben tener en cuenta los parámetros, pero el R2 cambios (creo) de una manera fundamental entre los dos.
- He intentado buscar en la CV para una respuesta a esta pregunta, pero no la encontró. Una solución alternativa que se propuso, pero su forma es sustancialmente diferente de lo que fue solicitado. Es sobre el masaje de las entradas, no sobre el formato de los comandos sin alterar fundamentalmente la de datos. La respuesta que me gustaba (y más útiles) es acerca de la forma de la fórmula introducida, no se trata de crear nuevas variables.
Como de costumbre, comentarios y sugerencias son solicitados.