[R] leverage

Thomas Lumley tlumley at u.washington.edu
Wed Aug 6 16:44:20 CEST 2003


On Wed, 6 Aug 2003, jane murray wrote:

> Hi
> Can anyone help with the technique of obtaining leverages from a
> conditional logistic regression model?  The code lm.influence does not seem
> to work for this data.

If you use method="approximate" you can get delta-betas (`influence'
rather than `leverage') with resid(model, "dfbeta")

For the exact conditional likelihood there isn't really a shortcut, but
if N is the number of observations:
  sapply(1:N, function(i) coef(update(model,subset=-i)))
is perfectly feasible.

	-thomas




More information about the R-help mailing list