GR2 <-
function(object)
{
if(class(object)[[1]] != "glmssn") return("Not a glmssn object")
W <- object$sampinfo$X
Vi <- object$estimates$Vi
z <- object$sampinfo$z
betahat <- object$estimates$betahat
muhat <- sum(Vi %*% z)/sum(Vi)
ones <- matrix(1, ncol = 1, nrow = length(z))
1 - t(z - W %*% betahat) %*% Vi %*% (z - W %*% betahat)/
t(z - ones %*% muhat) %*% Vi %*% (z - ones %*% muhat)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.