Search Archives

It's one of several possibilities ... Phillip On 19/12/18 18:26, Souheyla GHEBGHOUB wrote: > Is MCMCglmm the correct r function for multinomial multilevel mixed > logistic regression ? > > [[alternative HTML version deleted]] > > _______________________________________________ > R-sig-mixed-models at r-project.org...

I have no idea -- you may have to dig through the nlme source code to find out (https://github.com/cran/nlme/tree/master/R). nlme is largely in maintenance? mode at this point -- it's been 6 years since...

On 26/04/2020 11:06, Rolf Turner wrote: > Can the wise denizens of this list confirm to me the problem is *only* > the former? > > Be that as it may, is not still true that in general log likelihood is...

Can you show the output of sessionInfo() when run in vanilla R and RStudio? Phillip On 20/10/20 6:07 pm, Bingsong Zhang wrote: > Dear all, > > Did anyone encounter installation problems in R? I can successfully install > glmmTMB in...

On 16/5/21 11:28 pm, Ben Bolker wrote: > ?You can look at your results to figure out which coefficients are > missing, or you can use attr(getME(m,"X"), "col.dropped") to see which > columns of the original...

Yes. There will of course be greater uncertainty for conditions/groups with less data. The best way to see how this impacts your inference is to simulate balanced and unbalanced data for your hypothesis and look at the difference in...

This was discussed recently on the list in the context lme, but the same solution works for lme4: https://stat.ethz.ch/pipermail/r-sig-mixed-models/2021q1/029243.html Phillip On 21/01/2021 02:10, Sijia Huang wrote...

Dear all, We are happy to announce that registration is now open for the 2nd ed of Practical Mixed Effect Regression Modelling for Psychology and Language Science at the Radboud University summer school, taught by me and Laurel Brehm. The...

Without taking the time to look that up, you could also just use the lsmeans package, which has several different options for multiple-comparisons corrections and also uses the least-square means approach. Best, Phillip > On 17 Oct 2016, at...

See also https://doingbayesiandataanalysis.blogspot.com/2018/09/analyzing-ordinal-data-with-metric.html On 3/19/23 06:40, Karl Ove Hufthammer wrote: > Francesco Romano skreiv 18.03.2023 15:44: >> Thus, I would be very grateful to anyone...

Have you looked at the PyMC3 documentation? e.g. https://docs.pymc.io/notebooks/hierarchical_partial_pooling.html https://docs.pymc.io/notebooks/multilevel_modeling.html On 14/09/2019 00:33, John Poe wrote: > Hello everyone, > > I'm working...

I haven't tried Bambi just yet, but the PyMC3 package upon which it is built has been really quite nice in my experience. The only reason why I wasn't using PyMC3 more extensively is that it was a...

Yes. Well, to be more precise 0 + CAT and 1 + CAT estimate different contrasts and so the correlation you're estimating correspond to whatever contrast comes out. Same deal for 1 + CAT with sum vs treatment vs. Helmert coding. On...

For completeness, you can also get use getME(model, "beta") to get the fixed effects estimates as a simple, unnamed vector. :) On 17/12/19 4:40 am, Rolf Turner wrote: > > > On 17/12/19 1:49 am, Voeten, C...

You can split the specification of your grouping to achieve this, at least in part: lmer(y ~ A * B * C + (A * C | group) + (B|group) , data = data) Note that life gets tricky with the interaction terms. Phillip On 03/10...

Can you be more specific about the structure of your data and which object you're calling car::Anova() on? A small reproducible example is best, but details about the structure of the data and how you calculated the mixed...

This is a well-documented issue: || doesn't zero correlations between a categorical variable's levels. As far as I know, there are software-development/technical reasons for this, not statistical ones. The afex package has an implementation that zeroes...

See this section of the GLMM FAQ for a discussion of effect sizes and variable importance: https://bbolker.github.io/mixedmodels-misc/glmmFAQ.html#model-summaries-goodness-of-fit-decomposition-of-variance-etc. Phillip On 14/4/22 3:56...

This is generally not advisable because the sampling distribution of the variances is highly skewed. This means that the standard error isn't a particularly informative statistic.? This has been discussed in various amounts of detail in various places, but...

If I understand correctly, "fixed effects" in econometrics are simply categorical variables, especially ones with a large number of levels. There are "fixed" in the sense that they are observed at fixed (discrete) levels instead of as continuously. I don...