Nothing
#' mle_compet_lim
#' @keywords internal
#' @noRd
#' @importFrom bbmle parnames mle2
mle_het <- function(y, X, G, z, side, starting.values, boundL, boundU) {
theta <- starting.values
theta_L <- boundL
theta_U <- boundU
if(side == "het_l"){
bbmle::parnames(ll_het_left) <- names(theta)
fit <- mle3(ll_het_left, optimizer = "constrOptim",
method = "Nelder-Mead", start = theta, parnames = names(theta),
ui = rbind(c(rep(0, length(theta) - 3), - 1, - 1,0),
c(rep(0, length(theta) - 3), 1, 1, 0),
c(rep(0, length(theta) - 1), 1)),
ci = c( - 1 / max(colSums(G)), - 1 / max(colSums(G)), 0),
vecpar = TRUE, skip.hessian = TRUE,
data = list(Y = y, X = X, G = G, z = z))
} else if (side == "het_r"){
bbmle::parnames(ll_het_right) <- names(theta)
fit <- mle3(ll_het_right, optimizer = "constrOptim",
method = "Nelder-Mead", start = theta, parnames = names(theta),
ui = rbind(c(rep(0, length(theta) - 3), - 1, - 1, 0),
c(rep(0, length(theta) - 3), 1, 1, 0),
c(rep(0, length(theta) - 1), 1)),
ci = c( - 1 / max(rowSums(G)), - 1 / max(rowSums(G)), 0),
vecpar = TRUE, skip.hessian = TRUE,
data = list(Y = y, X = X, G = G, z = z))
}
return(fit)
}
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.