powered by
Calculate covariance matrix in the MVN-ICAR
.ICAR_calc(s, J, clam)
Sigma_s = (I - W)^(-1) * Q, W, Q
split points, J + 2
number of split points
controls neighbor interactions, in range (0, 1)