Sample ICAR random effects with component-wise sum-to-zero constraints
sample_icar(
A,
tau = 1,
kappa = 1e+06,
isolate = c("independent", "drop"),
isolate_prec = NULL
)Numeric vector of length nrow(A).
Sparse adjacency matrix (dgCMatrix preferred).
Positive precision multiplier.
Large penalty enforcing sum-to-zero (default 1e6).
How to handle isolated nodes: "independent" or "drop".
Precision for isolated nodes if kept independent.