I have some functions that uses an old version of the lme4 package where
mcmcsamp was implemented for binomial models. WHen I try to run it on the
latest lmer implementation I get messages saying "Update not yet written".

Is there any chance it may be written some time soon? The old version used
to work fine for getting reasonable confidence intervals for variance
components. If not are there alternatives that would not take too much
rewriting. Would MCMCsamp have similar syntax?

Here is some sample models from my functions.

#A univariate model:
    fitres1<-lmer(cbind(pos,neg)~1+(1|STUDY),
    family=binomial(link = "logit"),data=TPR )
    spl1<-mcmcsamp(fitres1,Nmcmc+5000)[-1:-5000,]

# a bivariate model
fitres<-lmer(cbind(pos,neg)~(GROUP)*(logit_-1)+(logit_-1|STUDY),
        family=binomial(link = "logit"),data=dd)
         spl<-mcmcsamp(fitres,Nmcmc+5000)[-1:-5000,]



With many thanks to all you R folk for making things available.

-- 
Gillian M Raab
10 Ainslie Place EH3 6AS
tel 0131 226 6234
mobile 07748 678 551

	[[alternative HTML version deleted]]


