[Lme4-commits] r1526 - branches/roxygen/tests

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Thu Jan 26 22:56:06 CET 2012


Author: dmbates
Date: 2012-01-26 22:56:06 +0100 (Thu, 26 Jan 2012)
New Revision: 1526

Modified:
   branches/roxygen/tests/vcov-etc.R
Log:
Use Matrix::diag explicitly.


Modified: branches/roxygen/tests/vcov-etc.R
===================================================================
--- branches/roxygen/tests/vcov-etc.R	2012-01-26 21:55:34 UTC (rev 1525)
+++ branches/roxygen/tests/vcov-etc.R	2012-01-26 21:56:06 UTC (rev 1526)
@@ -23,16 +23,17 @@
 	  all.equal(V.d <- vcov(fm1),
 		    V.s <- vcov(fm1.s)#, tol = 1e-9
 	  ),
-	  all.equal(diag(V.d), unname(se1.d)^2, tol= 1e-12)
+	  all.equal(Matrix::diag(V.d), unname(se1.d)^2, tol= 1e-12)
 	  ,
 	  all.equal(unname(se1.d),
-                    c(0.576012259837672, rep.int(0.518683836068417,3)),
+                    c(0.576011960125099, rep.int(0.518683987372292,3L)),
                     tol = 1.e-8)
           )
 
 ### -------------------------- a "large" example -------------------------
 str(InstEval)
 
+if (FALSE) {                            # sparse X is not currently implemented
 ## this works
 system.time(
 fm7 <- lmer(y ~ d + service + studage + lectage + (1|s),
@@ -65,4 +66,5 @@
 qqmath(rr,strip=FALSE)$d
 dotplot(rr,strip=FALSE)$`dept:service`
 }
+}
 cat('Time elapsed: ', proc.time(),'\n') # for ``statistical reasons''



More information about the Lme4-commits mailing list