Revision: 1230
http://r-gregmisc.svn.sourceforge.net/r-gregmisc/?rev=1230&view=rev
Author: warnes
Date: 2007-12-07 14:09:13 -0800 (Fri, 07 Dec 2007)
Log Message:
-----------
Add support for lme models to estimable()
Modified Paths:
--------------
trunk/gmodels/R/estimable.R
Modified: trunk/gmodels/R/estimable.R
===================================================================
--- trunk/gmodels/R/estimable.R 2007-12-07 22:07:26 UTC (rev 1229)
+++ trunk/gmodels/R/estimable.R 2007-12-07 22:09:13 UTC (rev 1230)
@@ -184,14 +184,21 @@
else if ("lm" %in% class(obj))
{
cf <- summary.lm(obj)$coefficients[, 1]
- vcv <- summary.lm(obj)$cov.unscaled * summary.lm(obj)$sigma^2
if ("glm" %in% class(obj))
- {
- vcv <- summary(obj)$cov.scaled
- }
+ vcv <- summary(obj)$cov.scaled
+ else
+ vcv <- summary.lm(obj)$cov.unscaled * summary.lm(obj)$sigma^2
}
+ else if ("lme" %in% class(obj))
+ {
+ s.o <- summary(obj)
+ cf <- s.o$tTable[,1]
+ se <- s.o$tTable[, 2]
+ rho <- s.o$cor
+ vcv <- rho * outer(se, se)
+ }
else
- stop("obj must be of class 'lm', 'glm', 'aov', 'gee' or 'geese'")
+ stop("obj must be of class 'lm', 'glm', 'aov', 'gee', 'geese', or 'lme'.")
u <- (cm %*% cf)-beta0
vcv.u <- cm %*% vcv %*% t(cm)
W <- t(u) %*% solve(vcv.u) %*% u
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
