Nothing
obj.cml <- function(para, map, data, ref, inv.V, bet0, sample.info, outcome){
para[map$all.bet] <- bet0
names(para)[map$all.bet] <- names(bet0)
ncase <- sample.info$ncase
nctrl <- sample.info$nctrl
nmodel <- length(map$bet)
lam <- para[map$lam]
the <- para[map$the]
fx <- as.matrix(data[, names(the), drop = FALSE])
y <- data[, outcome]
tilt <- tilt.cc(para, map, ref)
Delta <- tilt$Delta
delta <- tilt$delta
g <- gfunction.cc(para, map, ref, Delta, delta, ncase, nctrl)
pr <- as.vector(1/(1+g %*% lam))
plogL <- sum(y * log(Delta)) + sum(log(pr))
names(plogL) <- NULL
plogL
}
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.