This function updates the number of site level clusters for one patient level cluster via RJMCMC, when there is no patient in this cluster. Only used for updateRJ
function.
RJi_empty(w, K, Gamma, R, C,theta, tau, m, q, hyper_delta)
RJi_empty(w, K, Gamma, R, C,theta, tau, m, q, 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 current site level clustering membership
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 site level covariates
The hyper-parameter with default value being 1
Yuliang Li