Revision: 1521
http://r-gregmisc.svn.sourceforge.net/r-gregmisc/?rev=1521&view=rev
Author: warnes
Date: 2011-12-14 18:17:29 +0000 (Wed, 14 Dec 2011)
Log Message:
-----------
Improve formatting of ci.mer().
Modified Paths:
--------------
trunk/gmodels/R/ci.R
Modified: trunk/gmodels/R/ci.R
===================================================================
--- trunk/gmodels/R/ci.R 2011-12-14 18:14:03 UTC (rev 1520)
+++ trunk/gmodels/R/ci.R 2011-12-14 18:17:29 UTC (rev 1521)
@@ -68,26 +68,26 @@
retval
}
-ci.mer <- function(x,
- confidence=0.95,
- alpha=1-confidence,
- sim.mer=TRUE,
- n.sim=1000,
- ...
- )
+ci.mer <- function (x,
+ confidence = 0.95,
+ alpha = 1 - confidence,
+ sim.mer = TRUE,
+ n.sim = 1e4,
+ ...)
{
-## if(!(require(coda, quietly=TRUE) & require(Matrix, quietly=TRUE)))
-## stop("coda and Matrix packages required for ci.mer")
-
- x.effects <- fixef(x)
- n <- length(x.effects)
-
- retval <- est.mer(obj = x, cm = diag(n), beta0 = rep(0, n),
- conf.int = confidence, show.beta0 = FALSE,
- n.sim = n.sim)[,c("Estimate", "Lower.CI", "Upper.CI",
- "Std. Error", "p value")]
-
- colnames(retval)[c(2:3, 5)] <- c("CI lower", "CI upper", "p-value")
- rownames(retval) <- names(x.effects)
- retval
+ x.effects <- x@fixef
+ n <- length(x.effects)
+
+ retval <- gmodels:::est.mer(obj = x,
+ cm = diag(n),
+ beta0 = rep(0, n),
+ conf.int = confidence,
+ show.beta0 = FALSE,
+ n.sim = n.sim)
+
+ retval <- retval[, c("Estimate", "Lower.CI", "Upper.CI", "Std. Error", "p value")]
+ colnames(retval)[c(2:3, 5)] <- c("CI lower", "CI upper", "p-value")
+ rownames(retval) <- names(x.effects)
+
+ retval
}
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
