Nothing
gfunction.alp.lo <- function(para, map, data){
data$'(Intercept)' <- 1
nmodel <- length(map$bet)
g.alp <- list()
n <- nrow(data)
nlam <- max(map$lam)
offset <- max(map$the)
k <- 0
for(i in 1:nmodel){
id.a <- alp.index.lo(map, i)
alp.exist <- !is.null(id.a)
if(alp.exist){
alp <- para[id.a]
}else{
alp <- NULL
}
id.b <- map$bet[[i]]
bet <- para[id.b]
gam <- c(alp, bet)
rx <- as.matrix(data[, names(gam), drop = FALSE])
tmp <- as.vector(exp(rx %*% gam))
yi <- tmp/(1+tmp)
if(!alp.exist){
rm(id.a, id.b, alp.exist)
next
}
for(j in id.a){
rx0 <- rx[, names(para)[j]]
ga <- matrix(0, nrow = n, ncol = nlam)
ga[, id.a - offset] <- -rx[, names(alp), drop = FALSE] * (rx0 * yi * (1-yi))
ga[, id.b - offset] <- -rx[, names(bet), drop = FALSE] * (rx0 * yi * (1-yi))
k <- k + 1
g.alp[[k]] <- ga
rm(ga)
}
rm(id.a, id.b, alp.exist)
}
if(length(g.alp) == 0){
g.alp <- NULL
}
g.alp
}
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.