Giving 2 sub_models, build a correlated bivaraite operator based on K = D(theta, eta) $$D(\theta, \rho) = \begin{pmatrix} \cos(\theta) + \rho \sin(\theta) & -\sin(\theta) \sqrt{1+\rho^2} \\ \sin(\theta) - \rho \cos(\theta) & \cos(\theta) \sqrt{1+\rho^2} \end{pmatrix}$$ Exact same implementation as in paper. Fix noise sd=1.
bv_normal(
mesh,
sub_models,
rho = 0,
c1 = 1,
c2 = 1,
group = NULL,
share_param = FALSE,
fix_bv_theta = FALSE,
...
)
mesh for build the model
a list of sub_models (total 2 sub_models)
the parameter related to correlation
the noise sd for 1st sub_model
the noise sd for 2nd sub_model
group vector, can be inherited from ngme() function
TRUE if share the same parameter for 2 sub_models (of same type)
ignore
a list of specification of model