| NEWS | R Documentation |
It is now possible to use the MQL estimation technique as an alternative to PQL.
As an alternative to extending a logit model with random effects, it is now possible to add an over-dispersion parameter to logit model.
In addition to approximate the ML estimator, MQL and PQL have a variant that approximates the REML estimator.
There is now a simulate() method for objects returned
by mblogit() or mclogit() (but only for those
without random effects).
The algorithm for fitting random-effects models tended to stop prematurely returning the starting values obtained using a methods of moments. It has been completely refactored and proceeds similar to the PQL algorithm in Professor Brian Ripley's MASS package:
Now an inner step, in which a linear mixed model is fitted to a working dependent variable is nested into outer step iterations, in which the the working dependent variable is updated.
Also, the PQL algorithm no longer builds on starting values from a no-random-effects model, because surprisingly this makes the algorithm more stable and not less.
As a consequence, the algorithm does a much better job at avoiding divergence or running into numerical difficulties.
The PQL estimator for random-effects model uses a (inverse) Cholesky factor parametrisation, which makes sure that random-effects (Co-)Variance matrices are always positive (semi-)definite.
mclogit now complains about (i.e. throws an error
exception) when the random effects structure cannot be estimated,
e.g. because random effects are constant within choice sets
and therefore drop out by the formation of conditional logits.
mblogit now handles responses with only two columns.
mblogit now can deal with matrix responses that have
rows that sum to zero.
mclogit and mblogit now return a component
named "df.residual" instead of "residual.df".
mclogit now supports conditional logit models with
random slopes.
mblogit now supports multinomial baseline models with
random intercept and random slopes.
predict methods of objects created by mclogit
and mblogit are better in handling missing data.
New nobs and extractAIC methods for mclogit objects, so that drop1.default should work with these.
New function mblogit to fit multinomial baseline logit models.
mclogit mclogit.fit now allow user-provided starting values.
getSummary methods now return "contrasts" and "xlevels" components.
Fixed prediction method for mclogit results.
Corrected handling of weights and standard errors of prediction.
Matrices returned by the mclogit method of vcov() have
row and column names.
The number of observations is now displayed where it was not before.
nobs is defined in mclogit.fit.rePQL.
mclogit.fit and mclogit.fit.rePQL are exported to enable their use by other packages.