Hello all,
Does anyone know how to constrain/force specific coefficients when running
lm()?
I need to run recresid() {strucchange package} on the residuals of
forecast.lm, but forecast.lm's coefficients must be determined by
parameter.estimation.lm
I could estimate forecast.lm without lm() and use some other kind of
optimisation, but recresid() requires an object with class lm. recresid()
allows you to specify a formula, rather than an lm object, but it looks like
coefficients are estimated this way too and can't be forced.
Here is a bit of code to compensate for my poor explanation:.
# Estimate the coefficients of model
parameter.estimation.lm = lm(formula = y ~ x1 + x2, data =
estimation.dataset)
# How do I force the coefficients in forecast.lm to the coeff estimation
from parameter.estimation.lm??
forecast.lm = lm(formula = y ~ x1 + x2, data = forecast.dataset)
# Because I need recursive residuals from the application of the
coefficients from parameter.estimation.lm to a different dataset
recresid(forecast.lm)
Thanks in advance guys,
R.
[[alternative HTML version deleted]]
Rick Ram
2005-Jul-28 18:23 UTC
[R] Fwd: Forcing coefficents in lm(), recursive residuals, etc.
Resending cos I think this didn't get through for some reason... apologies
if it arrives twice!
---------- Forwarded message ----------
From: Rick Ram <r.ramyar@gmail.com>
Date: 28-Jul-2005 18:03
Subject: Forcing coefficents in lm(), recursive residuals, etc.
To: R-help <r-help@stat.math.ethz.ch>
Hello all,
Does anyone know how to constrain/force specific coefficients when running
lm()?
I need to run recresid() {strucchange package} on the residuals of
forecast.lm, but forecast.lm's coefficients must be determined by
parameter.estimation.lm
I could estimate forecast.lm without lm() and use some other kind of
optimisation, but recresid() requires an object with class lm. recresid()
allows you to specify a formula, rather than an lm object, but it looks like
coefficients are estimated this way too and can't be forced.
Here is a bit of code to compensate for my poor explanation:.
# Estimate the coefficients of model
parameter.estimation.lm = lm(formula = y ~ x1 + x2, data =
estimation.dataset)
# How do I force the coefficients in forecast.lm to the coeff estimation
from parameter.estimation.lm??
forecast.lm = lm(formula = y ~ x1 + x2, data = forecast.dataset)
# Because I need recursive residuals from the application of the
coefficients from parameter.estimation.lm to a different dataset
recresid(forecast.lm)
Thanks in advance guys,
R.
[[alternative HTML version deleted]]
Possibly Parallel Threads
- adding observations to lm for fast recursive residuals?
- R strucchange question: recursive-based CUSUM
- lm() and factors appending
- forecasting linear regression from lagged variable
- 'breackpoints' (package 'strucchange'): 2 blocking error messages when using for multiple regression model testing