Nothing
obj.lm <- function(para, map, data, ref, inv.V, bet0, outcome){
nmodel <- length(map$bet)
lam <- para[map$lam]
sigma <- para[map$the[1]]
the <- para[map$the[-1]]
fx <- as.matrix(data[, names(the), drop = FALSE])
y <- data[, outcome]
res <- y - as.vector(fx %*% the)
g <- gfunction.lm(para, map, ref)
pr <- as.vector(1/(1+g %*% lam))
n <- nrow(data)
bet <- para[map$all.bet]
plogL <- -n/2 * log(sigma) - 1/2/sigma * sum(res^2) + sum(log(pr)) -1/2 * as.vector(t(bet0 - bet) %*% inv.V %*% (bet0 - bet))
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.