Nothing
#' Maximum likelihood function for model party
#' @keywords internal
#' @noRd
#' @importFrom utils flush.console
ll_party <- function(params){
name.params <- names(params)
beta <- params[ - which(name.params %in% c("phi_within", "phi_between",
"sigma"))]
beta <- unlist(beta)
sigma <- unlist(params[which(name.params %in% "sigma")])
phi_within <- unlist(params[which(name.params %in% "phi_within")])
phi_between <- unlist(params[which(name.params %in% "phi_between")])
cat("pars: ", phi_within, phi_between, sigma, beta, '\n')
flush.console()
n <- nrow(G_within)
I <- diag(n)
X <- as.matrix(X)
sigma_2 <- sigma
reg_par <- solve_block(I - (phi_within * G_within) - (phi_between * G_between))
reg_par_2 <- solve_block(I - (phi_within * t(G_within)) -
(phi_between * t(G_between)))
omega <- sigma_2 * reg_par %*% reg_par_2
omega_inv = (I - (phi_within * t(G_within)) - (phi_between * t(G_between))) *
(I - (phi_within * G_within) - (phi_between * G_between)) / sigma_2
det_reg_par <- as.numeric(determinant(reg_par %*% reg_par_2,
logarithm = T)$modulus)
other_params <- Y - reg_par %*% (X %*% beta)
res <- - (n / 2) * log(sigma_2) - (1 / 2) * det_reg_par - (1/2) *
t(other_params) %*% omega_inv %*% other_params
print( - res)
}
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.