Nothing
residuals.mmglmlong1 <- function (object, ...){
subjects <- unique(object$longitude)
N <- length(subjects)
k <- length(object$y)
n <- k/N
err <- rep(NA, k)
for (i in 1:N){
a <- (object$longitude==subjects[i])
if (object$glmfamily$family=="binomial") size <- object$size[a]
else size <- NA
w <- mmglm1(object$y[a], object$Pi, object$delta,
object$glmfamily, object$beta, object$Xdesign[a,],
sigma=object$sigma, nonstat=object$nonstat,
size=size, msg=FALSE)
err[((i-1)*n+1):(i*n)] <- residuals(w)
}
return(err)
}
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.