[Lme4-commits] r1405 - in pkg/lme4Eigen: R man

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Fri Sep 23 11:14:15 CEST 2011


Author: mmaechler
Date: 2011-09-23 11:14:15 +0200 (Fri, 23 Sep 2011)
New Revision: 1405

Modified:
   pkg/lme4Eigen/R/lmer.R
   pkg/lme4Eigen/man/ranef.Rd
Log:
*must* use lme4Eigen::: inside devFun() [evaluated from C++]

Modified: pkg/lme4Eigen/R/lmer.R
===================================================================
--- pkg/lme4Eigen/R/lmer.R	2011-09-23 08:20:06 UTC (rev 1404)
+++ pkg/lme4Eigen/R/lmer.R	2011-09-23 09:14:15 UTC (rev 1405)
@@ -174,14 +174,14 @@
         parent.env(rho) <- parent.frame()
         devfun <- if (compDev) {
             function(theta)
-                .Call(glmerLaplace, pp$ptr, resp$ptr,
+                .Call(lme4Eigen:::glmerLaplace, pp$ptr, resp$ptr,
                       theta, u0, beta0, verbose, FALSE, tolPwrss)
         } else {
             function(theta) {
                 pp$u0 <- u0
                 pp$beta0 <- beta0
                 pp$theta <- theta
-                pwrssUpdate(pp, resp, verbose, tol=tolPwrss)
+                lme4Eigen:::pwrssUpdate(pp, resp, verbose, tol=tolPwrss)
                 resp$Laplace(pp$ldL2(), pp$ldRX2(), pp$sqrL(0))
             }
         }
@@ -197,14 +197,14 @@
             rho$control <- control
             devfunb <- if (compDev) {
                 function(pars)
-                    .Call(glmerLaplace, pp$ptr, resp$ptr, pars[dpars],
+                    .Call(lme4Eigen:::glmerLaplace, pp$ptr, resp$ptr, pars[dpars],
                           u0, pars[-dpars], verbose, TRUE, tolPwrss)
             } else {
                 function(pars) {
                     pp$u0 <- u0
                     pp$theta <- pars[dpars]
                     pp$beta0 <- pars[-dpars]
-                    pwrssUpdate(pp, resp, verbose, uOnly=TRUE, tol=tolPwrss)
+                    lme4Eigen:::pwrssUpdate(pp, resp, verbose, uOnly=TRUE, tol=tolPwrss)
                     resp$Laplace(pp$ldL2(), pp$ldRX2(), pp$sqrL(0))
                 }
             }
@@ -757,6 +757,7 @@
 	ans <- ans[whchL]
 
 	if (postVar) {
+            .NotYetUsed("postVar=TRUE")## FIXME
 	    vv <- .Call(reTrmsCondVar, re, sigma(object))
 	    for (i in seq_along(ans))
 		attr(ans[[i]], "postVar") <- vv[[i]]

Modified: pkg/lme4Eigen/man/ranef.Rd
===================================================================
--- pkg/lme4Eigen/man/ranef.Rd	2011-09-23 08:20:06 UTC (rev 1404)
+++ pkg/lme4Eigen/man/ranef.Rd	2011-09-23 09:14:15 UTC (rev 1405)
@@ -68,9 +68,11 @@
 fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), sleepstudy)
 fm3 <- lmer(diameter ~ (1|plate) + (1|sample), Penicillin)
 ranef(fm1)
+if(FALSE) { ##-- postVar=TRUE is not yet implemented -- FIXME
 str(rr1 <- ranef(fm1, postVar = TRUE))
 dotplot(rr1,scales = list(x = list(relation = 'free')))[["Subject"]]
 str(ranef(fm2, postVar = TRUE))
+}
 op <- options(digits = 4)
 ranef(fm3, drop = TRUE)
 options(op)



More information about the Lme4-commits mailing list