R/mle_compet_lim.R

In econet: Estimation of Parameter-Dependent Network Centrality Measures

#' mle_compet_lim
#' @keywords internal
#' @noRd
#' @importFrom  bbmle parnames mle2
mle_compet_lim <- function(y, X, G, z, starting.values, boundL, boundU, ind_fixed_effects) {

  theta <- starting.values
  theta_L <- boundL
  theta_U <- boundU

  if (is.null(ind_fixed_effects)) {
    X = X[, -1]
  }

  bbmle::parnames(ll_compet_lim) <- names(theta)
  fit <- mle2(ll_compet_lim, method = "L-BFGS-B", start = theta,
              parnames = names(theta), vecpar = TRUE, lower = theta_L,
              upper = theta_U, control = list(maxit = 1000),
              data = list(Y = y, X = X, G = G, z = z))

  return(fit)

}