[Vegan-commits] r686 - in pkg/vegan: R inst

[noreply at r-forge.r-project.org]

Author: jarioksa
Date: 2009-02-17 09:45:51 +0100 (Tue, 17 Feb 2009)
New Revision: 686

Modified:
   pkg/vegan/R/rad.null.R
   pkg/vegan/inst/ChangeLog
Log:
rad.null: failed with families gaussian and Gamma

Modified: pkg/vegan/R/rad.null.R
===================================================================
--- pkg/vegan/R/rad.null.R	2009-02-15 15:21:51 UTC (rev 685)
+++ pkg/vegan/R/rad.null.R	2009-02-17 08:45:51 UTC (rev 686)
@@ -9,14 +9,16 @@
     wt <- rep(1, nsp)
     if (nsp > 0) { 
         fit <- rev(cumsum(1/nsp:1)/nsp) * sum(x)
+        res <- dev.resids(x, fit, wt)
+        deviance <- sum(res)
         aic <- aicfun(x, nsp, fit, wt, deviance)
     }
     else {
         fit <- NA
         aic <- NA
+        res <- NA
+        deviance <- NA
     }
-    res <- dev.resids(x, fit, wt)
-    deviance <- sum(res)
     residuals <- x - fit
     rdf <- nsp
     p <- NA

Modified: pkg/vegan/inst/ChangeLog
===================================================================
--- pkg/vegan/inst/ChangeLog	2009-02-15 15:21:51 UTC (rev 685)
+++ pkg/vegan/inst/ChangeLog	2009-02-17 08:45:51 UTC (rev 686)
@@ -4,6 +4,8 @@
 
 Version 1.6-13 (opened Feb 15, 2009)
 
+	* rad.null: failed with family = gaussian and Gamma.
+
 	* simulate.rda: new method to simulate responses with random error
 	for an rda() result object. The function uses Normal error, and
 	estimates the sd separately for each species from their rda