[R-sig-ME] current r-forge version fails R CMD check ... ?

[Reinhold Kliegl]

Not sure this helps, but I was told once that the following should hold:

model at ST  # Check on proper convergence: no zeros on diagonal

The dataframe D clearly violates this expectation. I seem to recall
that I have seen different solutions for degenerated data before;
context was a comparison between SAS and lmer.

I do NOT get  different results across runs for various models that
are well behaved with respect to the above criterion.

Reinhold

On Mon, Aug 3, 2009 at 4:36 PM, Fabian
Scheipl<Fabian.Scheipl at stat.uni-muenchen.de> wrote:
> Comparison of m2 and m3 works for me.
>
> I use R-2.9.0 with lme4_0.999375-31,  Matrix_0.999375-30 on Windoze XP
> Professional SP2.
>
> Code is below.
>
> Best,
> Fabian
>
>> sessionInfo()
> R version 2.9.0 (2009-04-17)
> i386-pc-mingw32
>
> locale:
> LC_COLLATE=German_Germany.1252;LC_CTYPE=German_Germany.1252;LC_MONETARY=German_Germany.1252;LC_NUMERIC=C;LC_TIME=German_Germany.1252
>
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods   base
>
> other attached packages:
> [1] lme4_0.999375-31   Matrix_0.999375-30 lattice_0.17-25
>
> loaded via a namespace (and not attached):
> [1] grid_2.9.0  tools_2.4.1
>
>> all(replicate(50, {
> + set.seed(1)
> + D <-  data.frame(y= rnorm(20,10), ff = gl(4,5),
> +         x1=rnorm(20,3), x2=rnorm(20,7),
> +         x3=rnorm(20,1))
> + m2 <- lmer(y ~ x1 + (x2|ff), data = D)
> + m3 <- lmer(y ~ (x2|ff) + x1, data = D)
> +
> + identical(ranef(m2), ranef(m3))
> + }))
> [1] TRUE
>
>> set.seed(1)
>> all(replicate(50, {
> + D <-  data.frame(y= rnorm(20,10), ff = gl(4,5),
> +         x1=rnorm(20,3), x2=rnorm(20,7),
> +         x3=rnorm(20,1))
> + m2 <- lmer(y ~ x1 + (x2|ff), data = D)
> + m3 <- lmer(y ~ (x2|ff) + x1, data = D)
> +
> + identical(ranef(m2), ranef(m3))
> + }))
> [1] TRUE
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>