An embedded and charset-unspecified text was scrubbed... Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080818/93600a91/attachment.pl>
Hi Fr?nzi, mod.sum <- summary(model) mod.sum @ sigma Hank On Aug 18, 2008, at 9:34 AM, Fr?nzi Korner wrote:> Hi Oliver > > > > you once asked this question at the R help list. > > > > Hi all, > > > > I try to fit a glmm model with binomial distribution and I would to > > verify that the scale parameter is close to 1... > > > > the lmer function gives the following result : > > Estimated scale (compare to 1 ) 0.766783 > > > > But I would like to know how this estimation (0.766783) is performed, > > and I would > > like to know if it is possible to find this estimation with the > > different results given by the function lmer. > > > > > > Have you got an answer? I am now working with R 2.7.1.. This R- > version does > not give the scale parameter in the summary of an lmer-object. > > Therefore, I would like to calculate the scale parameter by hand, > but I did > not find out how. > > > > Thanks and best > > Fr?nzi > > > > > > ********************************************************************** > ** > > Dr. Fr?nzi Korner-Nievergelt > > oikostat - Statistische Analysen und Beratung > > Ausserdorf 43 > > CH - 6218 Ettiswil > > Tel.: +41 (0) 41 980 49 22 > > www.oikostat.ch > > > > ********************************************************************** > *** > > > > Schweizerische Vogelwarte > > CH - 6204 Sempach > > > > www.vogelwarte.ch > > ********************************************************************** > *** > > > > > [[alternative HTML version deleted]] > > <ATT00001.txt>
Fr?nzi Korner <fraenzi.korner <at> oikostat.ch> writes:> Have you got an answer? I am now working with R 2.7.1.. This R-version does > not give the scale parameter in the summary of an lmer-object. > > Therefore, I would like to calculate the scale parameter by hand, but I did > not find out how. >One quick answer is that you can retrieve the scale parameter by using the sigma function, which is not exported from the lme4 namespace (but probably? should be): lme4:::sigma(obj) The other answer is that you might want to post this on the r-sig-mixed list.
On Mon, Aug 18, 2008 at 9:27 AM, Ben Bolker <bolker at ufl.edu> wrote:> Fr?nzi Korner <fraenzi.korner <at> oikostat.ch> writes:>> Have you got an answer? I am now working with R 2.7.1.. This R-version does >> not give the scale parameter in the summary of an lmer-object.>> Therefore, I would like to calculate the scale parameter by hand, but I did >> not find out how.> One quick answer is that you can retrieve the scale parameter > by using the sigma function, which is not exported from the lme4 > namespace (but probably? should be):> lme4:::sigma(obj)I would export the function if I was more confident that this is the appropriate value for comparison (and also considered the name more carefully).> The other answer is that you might want to post this > on the r-sig-mixed list.I think that should be the R-SIG-Mixed-Models list.