Nothing
residuals.DirichletRegModel <- function(object, type=c("standardized",
# "score",
"composite",
"raw"),
...){
Y <- object$Y
fitted.vals <- fitted(object, mu=TRUE, alpha=TRUE, phi=TRUE)
M <- fitted.vals[["mu"]]
A <- fitted.vals[["alpha"]]
f <- fitted.vals[["phi"]]
raw.res <- Y - M
type <- match.arg(type)
wghts <- object$weights
V <- apply(M, 2, function(m) (m*(1-m))/(1+f) )
switch(type,
"standardized" = {
res <- raw.res/sqrt(V)
},
# "score" = {
# res <- digamma(colSums(A)) - digamma(A) + log(Y)
# },
"composite" = {
res <- rowSums((raw.res/sqrt(V))^2)
},
"raw" = {
res <- raw.res
}
)
return(res)
}
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.