Learn R Programming

BAREB (version 0.1.2)

RJi_empty: Function to update the number of site level clusters for one patient level cluster via RJMCMC, when there is no patient in this cluster

Description

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.

Usage

RJi_empty(w, K, Gamma,  R,  C,theta, tau, m,  q, hyper_delta)

Value

RJi_empty(w, K, Gamma, R, C,theta, tau, m, q, hyper_delta) returns a list with following variables:

K

The numbers of site level clusters

w

The weights for site level clusters

Gamma

Linear coefficients for site level covariates

R

The site level clusteirng membership

C

The DPP related kernel matrice

Arguments

w

The weights for site level clusters. Should be a vector.

K

The number of site level cluster for the patient level cluster.

Gamma

The linear coefficients for site level covariates. Should be a matrix.

R

The current site level clustering membership

C

The DPP related kernel matrices. Should be an array of 3 dimensions

theta

The DPP hyper-parameter for site level

tau

A fixed DPP hyper-parameter, which we suggest high value, say 10^5

m

The number of sites

q

The number of site level covariates

hyper_delta

The hyper-parameter with default value being 1

Author

Yuliang Li