NSAPH-Projects/BNP-For-Unmeasured-Confounding

Questions

Closed this issue a year ago · 0 comments

katehu commented a year ago

BNP-For-Unmeasured-Confounding/pm25-cardiovasculardisease-example/2_model_bnp_2.R

Line 198 in 3984375

for (l in which(check_q == 1 | check_q == 2)) {

Should we use

which(check_q <3)

what if check_q ==0 ?

BNP-For-Unmeasured-Confounding/pm25-cardiovasculardisease-example/2_model_bnp_2.R

Line 200 in 3984375

m_eta = cbind(m_eta, rep(eta_prior[1] / eta_prior[2], 2))

Why we repeat (eta_prior[1] / eta_prior[2]) twice for each l?