Nothing
RaoScoreStatistic <- function(fm, X, dispersion = 1) {
fam <- fm$family
mus <- fm$fitted
tots <- fm$prior.weights
variances <- tots * fam$variance(mus)
dmu.deta <- tots * fam$mu.eta(fm$linear.predictor)
Info <- crossprod((D <- X * dmu.deta)/sqrt(variances))/dispersion
qScores <- t(D/variances) %*% ((fm$y - mus) * tots)/dispersion
t(qScores) %*% chol2inv(chol(Info)) %*% qScores
}
ordinaryDeviance <- function(fm, dispersion = 1) {
LP <- fm$linear.predictor
y <- fm$y
fam <- fm$family
mu <- fam$linkinv(LP)
wt <- fm$prior.weights
sum(fam$dev.resid(y, mu, wt))/dispersion
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.