Optimized Mixed Model Solver
mixed_solve_optimized.RdAn optimized implementation of mixed model estimation that precomputes expensive matrix operations and can be reused across multiple voxels for significant performance improvements.
Usage
mixed_solve_optimized(
X,
Z,
Y,
K = NULL,
workspace = NULL,
compute_se = FALSE,
n_threads = 0
)Arguments
- X
Fixed effects design matrix (n × p)
- Z
Random effects design matrix (n × q)
- Y
Response data - can be a vector (single voxel) or matrix (n × V for multiple voxels)
- K
Kinship/covariance matrix for random effects (q × q). Defaults to identity.
- workspace
Precomputed workspace (optional, will compute if NULL)
- compute_se
Whether to compute standard errors (default: FALSE)
- n_threads
Number of OpenMP threads for multi-voxel (0 = auto)