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.
log_likelihood2(params, X, Z, y, n, weights = NULL)
The negative log-likelihood of the model.
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 in the beta-binomial distribution).
A numeric vector of weights for each observation. If NULL, equal weights are used.
This function uses a numerically stable implementation of the beta-binomial log-probability. It allows for weighted observations, which can be useful for various modeling scenarios.