[Lme4-commits] r1718 - pkg/lme4/tests

[noreply at r-forge.r-project.org]

Author: mmaechler
Date: 2012-05-07 18:26:48 +0200 (Mon, 07 May 2012)
New Revision: 1718

Modified:
   pkg/lme4/tests/PIRLSfail.R
Log:
- the deviance function at least is available
- cosmetic


Modified: pkg/lme4/tests/PIRLSfail.R
===================================================================
--- pkg/lme4/tests/PIRLSfail.R	2012-05-04 15:42:36 UTC (rev 1717)
+++ pkg/lme4/tests/PIRLSfail.R	2012-05-07 16:26:48 UTC (rev 1718)
@@ -1,12 +1,20 @@
 ## example from HSAUR2 package
-## library("multcomp")
-## trees513 <- subset(trees513, !species %in% c("fir", "ash/maple/elm/lime", "softwood (other)"))
-## trees513$species <- trees513$species[,drop = TRUE]
-## levels(trees513$species)[nlevels(trees513$species)] <- "hardwood"
-## save("trees513",file="trees513.RData")
+if(FALSE) {
+    library("multcomp")
+    trees513 <- subset(trees513, !species %in%
+		       c("fir", "ash/maple/elm/lime", "softwood (other)"))
+    trees513$species <- trees513$species[,drop = TRUE]
+    levels(trees513$species)[nlevels(trees513$species)] <- "hardwood"
+    save("trees513", file="trees513.RData")
+}
+
 library("lme4")
 load("trees513.RData")
+dfun <- glmer(damage ~ species - 1 + (1 | lattice / plot),
+              data = trees513, family = binomial, devFunOnly = TRUE)
+ls.str(environment(dfun))# and you can investigate ...
+
 ## FIXME: try() to pass checks
-mmod <- try(lmer(damage ~ species - 1 + (1 | lattice / plot),
+mmod <- try(glmer(damage ~ species - 1 + (1 | lattice / plot),
               data = trees513, family = binomial()))