This function estimates the theta4 parameter using the Millar method.
Usage
estimate_millar_theta4(
class_predictions_mixed_sample,
PHI_matrix,
np,
use_accelerated_em = TRUE,
...
)
Arguments
- class_predictions_mixed_sample
Numeric vector of predicted classes for the mixed sample
- PHI_matrix
Confusion matrix used for correction
- np
Integer number of populations
- use_accelerated_em
Logical; whether to use the accelerated EM version
- ...
Additional arguments passed to internal functions
Value
Numeric value of the estimated theta4 parameter
Examples
if (FALSE) { # \dontrun{
pred_classes <- c(1, 2, 1, 2)
phi_mat <- matrix(c(0.8, 0.2, 0.2, 0.8), 2, 2)
theta4 <- estimate_millar_theta4(pred_classes, phi_mat, np = 2)
} # }