Sample subject specific random effects
b_samp(Omega, beta, R, theta, X, q, Z, ZtZ, N)
current covariance matrix
current fixed effect coefficients vector
component of current outcome vector (Y.I=cos(theta)*R, Y.II=sin(theta)*R)
design matrix model parameters (differs per person)
dimension Z(number of random effects)
design matrix for random effects (differs per person)
transpose(Z)*Z
sample size at second level