Computes the sum of the vector deviance and other intermediate results
glmdev(y, ni, ci, wa, vtheta, offset = 0, icase = .dFvGet()$ics)A list with the following components:
2*sum_i abs(Li-Ti)
The estimates of theta_i
The values of Li
The alues of Ti
The vector of observations
The number of trial at xi in the binomial case (ics=2). Otherwise ni=1 for each xi.
The constants ci
The vector of ai=b/|Axi|
The vector of xi^T
Optional offset added to the linear predictor.
Set ics=1 for Bernoulli case, ics=2 for Binomial case and ics=3 for Poisson case
Kuensch, H.R., Stefanski L.A., Carroll R.J. (1989). Conditionally unbiased bounded-influence estimation in general regression models, with application to generalized linear models. Journal of the American Statistical Association, 84, 460-466.
Marazzi, A. (1993). Algorithms, Routines, and S-functions for robust Statistics. Chapman and Hall, New York.
Marazzi A. (1997). Object oriented S-plus functions for robust discrete generalized linear models available in the doc folder of this package.