Nothing
## uses untransform()
getD <- function (designD, beta, mask, parindx, link, fixed,
grouplevels, sessionlevels, parameter = 'D') {
if ((is.null(designD) | nrow(designD)==0) & (is.null(fixed[[parameter]]))) return(NULL)
if (ms(mask))
nmask <- max(sapply(mask, nrow))
else
nmask <- nrow(mask)
ngroup <- length(grouplevels)
nsession <- length(sessionlevels)
D <- array(dim = c(nmask, ngroup, nsession))
dimnames(D) <- list(1:nrow(D), grouplevels, sessionlevels)
if (!is.null(fixed[[parameter]])) {
D[,,] <- fixed[[parameter]]
}
else {
D[,,] <- designD %*% beta[parindx[[parameter]]] # linear predictor
D[,,] <- untransform (D, link[[parameter]])
# silently truncate D at zero
if (parameter %in% c('D')) D[D<0] <- 0
}
D
}
###############################################################################
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.