[R-gregmisc-users] SF.net SVN: r-gregmisc: [1230] trunk/gmodels/R/estimable.R

[<wa...@us...>]

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.