www.cmm.bristol.ac.uk/forum

Thanks for the investigation and reply. Would it be ok to use:

binomialMCMC@residual$lev_2_resi_est_Intercept<-binomialMCMC@residual$lev_2_resi_est_Intercept/2

to replace the level 2 residuals? I would like to be able to produce some plots using caterpillarR.

Hi Chris.

I ran a slightly modified version of one of the exercises from the R2MLwiN JSS article (Zhang et al., 2016, p.26-27) as follows:

library(R2MLwiN)
data("bang1")
F6 = logit(use, cons) ~ 1 + age + lc + urban + (1|district)
binomialMCMC <- runMLwiN(Formula = F6, D = "Binomial", data = bang1 ...

Thanks for your reply. I have run the model in debugmode and managed to compare the residuals and also open the relevant temp .dta file. It seems as though the R2MLwiN output is multiplying the level 2 residuals by 2.

Is there something in particular that I could attach or provide to help ...

Hello.

I have been fitting a cross-classified logistic regression model using MCMC for a binary response variable and would like to extract the estimated level 2 residuals and match these to the level 2 group names. However, the residuals I get from using the resi.store=TRUE option in R2MLwiN are ...