[Lme4-commits] r1471 - testing/lme4Eigen

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

Author: dmbates
Date: 2011-12-06 18:42:13 +0100 (Tue, 06 Dec 2011)
New Revision: 1471

Added:
   testing/lme4Eigen/optimxcbpp.R
Log:
Check settings of tolPwrss and various optimx algorithms on the cbpp example.


Added: testing/lme4Eigen/optimxcbpp.R
===================================================================
--- testing/lme4Eigen/optimxcbpp.R	                        (rev 0)
+++ testing/lme4Eigen/optimxcbpp.R	2011-12-06 17:42:13 UTC (rev 1471)
@@ -0,0 +1,17 @@
+library(optimx)
+library(lme4Eigen)
+form <- cbind(incidence, size - incidence) ~ period + (1 | herd)
+
+ff <- glmer(form, cbpp, binomial, devFunOnly=2L)
+environment(ff)$tolPwrss                # current default (but should be changed)
+optimx(with(environment(ff), c(pp$theta, pp$beta0)), ff,
+       lower=c(0, rep.int(-Inf, 4L)), control=list(all.methods=TRUE))
+
+ff <- glmer(form, cbpp, binomial, devFunOnly=2L, tolPwrss=1e-9)
+optimx(with(environment(ff), c(pp$theta, pp$beta0)), ff,
+       lower=c(0, rep.int(-Inf, 4L)), control=list(all.methods=TRUE))
+
+ff <- glmer(form, cbpp, binomial, devFunOnly=2L, tolPwrss=1e-10)
+optimx(with(environment(ff), c(pp$theta, pp$beta0)), ff,
+       lower=c(0, rep.int(-Inf, 4L)), control=list(all.methods=TRUE))
+