Nothing
pearson2 <- function(object, ...) {
dispersion <-
with(object, {
y.hat <- predict(y, coefficients[1:p], X[,1:p], offset)
a.hat <- predict_s(y, coefficients[-(1:p)], X[,-(1:p)])
pearson <- sum(((y - y.hat)^2) /
variance(y, y.hat, a.hat))
pearson / df.residual})
return(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.