When selecting a linear regression model it is sometimes useful to drop independent variables from the model and to examine the resulting sub-model.
g02dfc updates the
decomposition used in the computation of the linear regression model. The
decomposition may come from
g02dac,
g02dcc,
g02dec or a previous call to
g02dfc.
For the general linear regression model with
independent variables fitted,
g02dac or
g02dec computes a
decomposition of the (weighted) independent variables and forms an upper triangular matrix
and a vector
. To remove an independent variable
and
have to be updated. The column of
corresponding to the variable to be dropped is removed and the matrix is then restored to upper triangular form by applying a series of Givens rotations. The rotations are then applied to
. Note that only the first
elements of
are affected.
The method used means that while the updated values of
and
are computed an updated value of
from the
decomposition is not available so a call to
g02dec cannot be made after a call to
g02dfc.
g02ddc can be used to calculate the parameter estimates,
, from the information provided by
g02dfc.
Hammarling S (1985) The singular value decomposition in multivariate statistics SIGNUM Newsl. 20(3) 2–25
There will inevitably be some loss in accuracy in fitting a model by dropping terms from a more complex model rather than fitting it afresh using
g02dac.
Background information to multithreading can be found in the
Multithreading documentation.
None.
A dataset consisting of 12 observations on four independent variables and one dependent variable is read in. The full model, including a mean term, is fitted using
g02dac. The value of
indx is read in and that variable dropped from the regression. The parameter estimates are calculated by
g02ddc and printed. This process is repeated until
indx is 0.