[R-sig-ME] lmer, poisson family and mcmcsamp

Amelie LESCROEL lescroel at cebc.cnrs.fr
Tue Feb 17 10:55:51 CET 2009


Thanks Chris, I saw this one but I was wondering whether the situation 
changed or was about to change. What about Pr(>|z|)? Are they reliable?
Amelie

Christopher David Desjardins a écrit :
> Regarding mcmcsamp() see 
> https://mailman.stat.ethz.ch/pipermail/r-sig-mixed-models/2008q4/001660.html specifically 
>
>
> ...  In particular mcmcsamp is not currently working in that branch.  
> My goal is to get mcmcsamp stable and to get all the other parts 
> working in this branch then release it for testing and finally release 
> lme4-1.0.  The subtext in that message is that only the models I 
> describe above will be fit by the lme4-1.0 package.
>
>
> On Feb 16, 2009, at 9:59 AM, Amelie LESCROEL wrote:
>
>> Dear all,
>>
>> Could someone explain to me why summary(lmer.object) is giving z 
>> values and Pr(>|z|) when the lmer.object was fitted with 
>> family=poisson? I used to fit linear mixed models with 
>> family=gaussian and link="identity" or link="log", then used 
>> mcmcsamp() and HPDinterval() for computing "confidence" intervals and 
>> considered that a fixed effect was significant when the 95% HPD 
>> interval did not include zero. Now, with family=poisson, when I use 
>> mcmcsamp(), I get the following:
>>
>> > set.seed(101); mcmcfm1 <- mcmcsamp(fm1, n=10000)
>> Error in .local(object, n, verbose, ...) : Update not yet written
>>
>> but I get z values and Pr(>|z|). Why is it so? Can we "trust" this 
>> probability for judging of the fixed effects significance? Would the 
>> mcmcsamp method be available soon for families other than gaussian?
>>
>> Thanks in advance,
>>
>> Amélie
>>
>> PS: When will lme4_1.0 be released? Many thanks to Douglas Bates for 
>> developing lme4 in interaction with users, and for sharing it. Please 
>> don't give up on making it available in R. lme4 is a wonderful tool 
>> which is helping me understanding penguin's foraging behaviour (and 
>> advanced stats!) as it is helping others in plant ecology or human 
>> psychology.
>>
>> > sessionInfo()
>> R version 2.7.2 (2008-08-25)
>> i386-pc-mingw32
>>
>> locale:
>> LC_COLLATE=French_France.1252;LC_CTYPE=French_France.1252;LC_MONETARY=French_France.1252;LC_NUMERIC=C;LC_TIME=French_France.1252 
>>
>>
>> attached base packages:
>> [1] stats     graphics  grDevices utils     datasets  methods   base
>> other attached packages:
>> [1] lme4_0.999375-28   Matrix_0.999375-17 lattice_0.17-17    gdata_2.4.2
>> loaded via a namespace (and not attached):
>> [1] grid_2.7.2   gtools_2.5.0 tools_2.7.2
>>
>>
>>
>> __________ Information from ESET Mail Security, version of virus 
>> signature database 3857 (20090216) __________
>>
>> The message was checked by ESET Mail Security.
>> http://www.eset.com
>>
>> _______________________________________________
>> R-sig-mixed-models at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
> -----------------
> Christopher David Desjardins
> Ph.D. Student
> Quantitative Methods in Education
> Department of Educational Psychology
> University of  Minnesota
> http://blog.lib.umn.edu/desja004/educationalpsychology/
>
>
>
>
>
>
>
> __________ Information from ESET Mail Security, version of virus 
> signature database 3857 (20090216) __________
>
> The message was checked by ESET Mail Security.
> http://www.eset.com
>
>
>
>
>



__________ Information from ESET Mail Security, version of virus signature database 3860 (20090217) __________

The message was checked by ESET Mail Security.
http://www.eset.com




More information about the R-sig-mixed-models mailing list