This function updates the number of site level clusters for one patient level cluster via RJMCMC. Only used for updateRJ function.
RJi(w, K, Gamma, Beta, X, Y, Z, R, delta, mu, mu_star, c, sigma_square, C,
theta, tau, m, n,q,T0, hyper_delta)RJi(w, K, Gamma, Beta, X, Y, Z, R, delta, mu, mu_star, c, sigma_square, C, theta, tau, m, n,q,T0, hyper_delta) returns a list with following variables:
The numbers of site level clusters
The weights for site level clusters
Linear coefficients for site level covariates
The site level clusteirng membership
The DPP related kernel matrice
The weights for site level clusters. Should be a vector.
The number of site level cluster for the patient level cluster.
The linear coefficients for site level covariates. Should be a matrix.
The linear coefficients for patient level covariates. Should be a vector.
The design matrix for patient level clusters.
The observed CAL value matrix
The design matrix for site level clusters.
The current site level clustering membership
The missing indicator matrix
The current CAL mean vector
The latent value vector for missingness model
The linear coefficients for missingness model
The current noise variance
The DPP related kernel matrices. Should be an array of 3 dimensions
The DPP hyper-parameter for site level
A fixed DPP hyper-parameter, which we suggest high value, say 10^5
The number of sites
The number of patients
The number of site level covariates
The number of teeth
The hyper-parameter with default value being 1
Yuliang Li