Search Archives

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20140408/9dc48f3e/attachment-0001.pl>

Hello, is it possible to install the flexlambda and master branch of lme4 at the same time: library(devtools) install_github("lme4", "lme4") install_github("lme4", "flexlambda") such that I can load them as different packages. (Obviously not both at...

Hello, I am currently testing whether I should include certain random effects in my lmer model or not. I use the anova function for that. My procedure so far is to fit the model with a function call to lmer...

Hello, For testing/research purposes I compiled R from source with blas as a shared library. I then went on to compile openblas from source, tuned it to Sandybridge and linked R against. The crucial step being: cd /usr/local...

Dear Ben, dear Tobias, the devtools solution is actually pretty nice: library(devtools) dev_mode(, path = "~/your/path/to/new/library/here") will switch your R repl from > to d> you can then simply install install_github("lme4", "lme4", ref...

Hello, fitting linear mixed models it is often suggested that testing for random effects is not the best idea; mainly because the value of the random effects parameters lie at the boundary of the parameter space. Hence, it is preferred...

Hello all, warning ahead: This will be quite a long post. But have faith, it is necessary. I run parametric bootstraps in order to compare two models. I bootstrap the likelihood ratio between the models. For anybody acquainted with this...

Hello, I am performing a priori power simulations for mixed-effect models based on previous experiments. This works out quite nicely. I extract parts of my parameters from a previous model I fitted: prev_mod <- lmer(Y ~ A + (B | Context...

Hi Malcolm, Interesting. I just read the reference manual for the ordinal package (v. July 2, 2014) and indeed at the beginning it states "random slopes (random coefficient models) are not yet implemented" (p. 3). If this i indeed true...

Dear Vince, Ben and Steve, thank you! I'm reading the paper and the suggested sections. It's really good and really helpful! Just to make sure I understand you correctly: The individual variance-covariance matrices for the random effects...

Hello, I've posted yesterday investigating about random effects in proportional odds cummulative mixed effects model; specifically about clmm() from the ordinal package. I was doing some more reading about categorical data analysis (e.g. Agresti (2010): Analysis of Ordinal...

Hi Jarrod, Hi Malcolm, Thank you Jarrod! In the meantime I wrote Rune and asked him whether clmm() is currently able to fit random slope models and if he could update his reference manual should this be the case. He...

Hello all, I fitted a proportional ordered logit model and I have some questions regarding my interpretation of the R output and calculating certain values. I will illustrate this with the wine dataset which can be found in the ordinal...

Thanks Ben! Figured it out myself pretty quickly. Seems you also have to transpose: t(chol(VarCorr(fm1)[[1]]/sigma(fm1)) if you want the representation to be identical to "getME(fm1, "Tlist")[[1]] Thanks again Vince, Steve and Ben...

I ran some more tests and I can safely say that the culprit in Method 2 and Method 3 is the call to: dattemp <- mod at frame and subsequently using it as the data argument in the foreach() loop calls...

> * is the difference between p=0.004 and p=0.006 _practically_ > important? No, I don?t think that it is important at all. Mostly because I don?t rely on p values. (But a lot of people still like...