[Vegan-commits] r2766 - pkg/vegan/R

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

Author: jarioksa
Date: 2013-11-28 15:09:19 +0100 (Thu, 28 Nov 2013)
New Revision: 2766

Modified:
   pkg/vegan/R/anovacca.byterm.R
Log:
you must use constrained inertia with model='full' (it does not matter with other 'model's)

Modified: pkg/vegan/R/anovacca.byterm.R
===================================================================
--- pkg/vegan/R/anovacca.byterm.R	2013-11-27 18:17:25 UTC (rev 2765)
+++ pkg/vegan/R/anovacca.byterm.R	2013-11-28 14:09:19 UTC (rev 2766)
@@ -69,11 +69,11 @@
     Chisq <- sapply(mods, function(x) x$chi[2]) - chibig
     Fstat <- (Chisq/Df)/(chibig/dfbig)
     ## Simulated F-values
-    Fval <- sapply(mods, function(x) x$den)
+    Fval <- sapply(mods, function(x) x$num)
     ## Had we an empty model we need to clone the denominator
     if (length(Fval) == 1)
         Fval <- matrix(Fval, nrow=nperm)
-    Fval <- sweep(Fval, 1, big$den, "-")
+    Fval <- sweep(-Fval, 1, big$num, "+")
     Fval <- sweep(Fval, 2, Df, "/")
     Fval <- sweep(Fval, 1, scale, "/")
     ## Simulated P-values