R/ll_compet_lim.R

In econet: Estimation of Parameter-Dependent Network Centrality Measures

#' Maximum likelihood function for model compet_het
#' @keywords internal
#' @noRd
#' @importFrom utils flush.console
ll_compet_lim <- function(params){

  name.params <- names(params)
  beta <- params[ - which(name.params %in% c("alpha", "phi", "sigma"))]
  beta <- unlist(beta)
  sigma <- unlist(params[which(name.params %in% "sigma")])
  phi <- unlist(params[which(name.params %in% "phi")])
  alpha <- unlist(params[which(name.params %in% "alpha")])

  cat("pars: ", alpha, phi, sigma, beta, '\n')
  flush.console()

  n <- nrow(G)
  I <- diag(n)
  X <- as.matrix(X)
  z <- diag(z)

  sigma_2 <- sigma
  reg_par <- solve_block(I - phi * G)

  other_params <- Y - alpha * (reg_par %*% rep(1, n)) - (X %*% beta)

  res <- - (n / 2) * log(sigma_2)  - (1 / 2) *
    t(other_params) %*% other_params / sigma_2

  print( - res)
}