[R] Mac/PC differences in lmer results
bgunter@4567 @end|ng |rom gm@||@com
Thu May 30 16:34:57 CEST 2019
The BLAS in use on each?
"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Thu, May 30, 2019 at 5:27 AM Nicolas Schuck <nico.schuck using gmail.com>
> Dear fellow R coders,
> I am observing differences in results obtained using glmer when using a
> Mac or Linux computer versus a PC. Specifically, I am talking about a
> relatively complex glmer model with a nested random effects structure. The
> model is set up in the following way:
> gcctrl = glmerControl(optimizer=c('nloptwrap'), optCtrl = list(maxfun =
> 500000), calc.derivs = FALSE)
> glmer_pre_instr1 = glmer(
> formula = cbind(FREQ, NSAMP-FREQ) ~ FDIST_minz + poly(RFREQ,2) + ROI
> + (1 + FDIST_minz + RFREQ + ROI|ID/COL),
> data = cdf_pre_instr,
> family = binomial,
> control = gcctrl)
> Code and data of an example for which I find reproducible, non-negligible
> differences between Mac/Win can be found here:
> https://gitlab.com/nschuck/glmer_sandbox/tree/master <
> The differences between the fitted models seem to be most pronounced
> regarding the estimated correlation structure of the random effects terms.
> Mac and Linux yield very similar results, but Windows deviates quite a bit
> in some cases. This has a large impact on p values obtained when performing
> model comparisons. I have tried this on Mac OS 10.14, Windows 10 and Ubuntu
> and Debian. All systems I have tried are using lme 1.1.21 and R 3.5+.
> Does anyone have an idea what the underlying cause might be?
> [[alternative HTML version deleted]]
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
[[alternative HTML version deleted]]
More information about the R-help