Evaluates the loglikelihood of a Poisson-Tweedie generalized linear mixed model with random intercept, using the adaptive Gauss-Hermite quadrature rule.
loglik.pt.1re(beta, D, a, Sigma, y, X, Z, id, offset = NULL, GHk = 10,
tol = 9.88131291682493e-324, GHs = NULL)
Vector of regression coefficients
Dispersion parameter (must be > 1)
Power parameter (must be < 1)
A matrix with the variance of the random intercept
Response vector (discrete)
Design matrix for the fixed effects
Design matrix for the random effects
Id indicator (it should be numeric)
Offset term to be added to the linear predictor
Number of quadrature points (default is 10)
Tolerance value for the evaluation of the probability mass function of the Poisson-Tweedie distribution
Quadrature points at which to evaluate the loglikelihood. If NULL
(default value),
the GH quadrature points are computed internally
The loglikelihood value obtained using a Gauss-Hermite quadrature
approximation with GHk
quadrature points.
ptmixed
and the examples therein