This function computes the negative log-likelihood for a beta-binomial regression model where both the alpha and beta parameters are modeled as functions of predictors (mode 2).
log_likelihood2(params, X, Z, y, n, weights = NULL, lch = NULL)The negative log-likelihood of the model (large finite penalty if non-finite).
A numeric vector containing all model parameters. The first n_alpha elements are coefficients for the alpha model, and the remaining elements are coefficients for the beta model.
A matrix of predictors for the alpha model.
A matrix of predictors for the beta model.
A numeric vector of response values.
The maximum score (number of trials).
A numeric vector of weights for each observation (NULL = equal weights).
Optional precomputed lchoose(n, y). Since this term does
not depend on the parameters, passing it once avoids recomputation in
every optimizer iteration.
Uses a numerically stable implementation of the beta-binomial
log-probability via lbeta. The linear predictors of
log(alpha) and log(beta) are clamped to [-20, 20].