[Vegan-commits] r491 - in pkg: R inst

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Mon Sep 1 16:45:07 CEST 2008


Author: jarioksa
Date: 2008-09-01 16:45:07 +0200 (Mon, 01 Sep 2008)
New Revision: 491

Modified:
   pkg/R/permutest.cca.R
   pkg/R/print.permutest.cca.R
   pkg/inst/ChangeLog
Log:
permutest.cca tells if it analysed all eigenvalues or only the first

Modified: pkg/R/permutest.cca.R
===================================================================
--- pkg/R/permutest.cca.R	2008-08-31 17:04:27 UTC (rev 490)
+++ pkg/R/permutest.cca.R	2008-09-01 14:45:07 UTC (rev 491)
@@ -88,7 +88,7 @@
     }
     sol <- list(call = x$call, model = model, F.0 = F.0, F.perm = F.perm, 
                 chi = c(Chi.z, Chi.xz), num = num, den = den, df = c(q, 
-                                                              r), nperm = permutations, method = x$method, Random.seed = seed)
+                                                              r), nperm = permutations, method = x$method, first = first,  Random.seed = seed)
     if (!missing(strata)) {
         sol$strata <- deparse(substitute(strata))
         sol$stratum.values <- strata

Modified: pkg/R/print.permutest.cca.R
===================================================================
--- pkg/R/print.permutest.cca.R	2008-08-31 17:04:27 UTC (rev 490)
+++ pkg/R/print.permutest.cca.R	2008-09-01 14:45:07 UTC (rev 491)
@@ -4,7 +4,11 @@
     cat("\nPermutation test for", x$method, "\nCall:\n")
     cat(deparse(x$call), "\n\n")
     Pval <- sum(x$F.perm >= x$F.0)/x$nperm
-    cat("Test for significance of all constrained eigenvalues\n")
+    cat("Permutation test for ")
+    if (x$first)
+        cat("first constrained eigenvalue\n")
+    else
+        cat("all constrained eigenvalues\n")
     cat("Pseudo-F:\t", x$F.0, "\n")
     cat("Significance:\t", format.pval(Pval, eps = 1/x$nperm), 
         "\n")

Modified: pkg/inst/ChangeLog
===================================================================
--- pkg/inst/ChangeLog	2008-08-31 17:04:27 UTC (rev 490)
+++ pkg/inst/ChangeLog	2008-09-01 14:45:07 UTC (rev 491)
@@ -29,10 +29,12 @@
 	needs evaluate the community within calling function. (The
 	environment() drives me absolutely crazy, says J.O.)
 
-	* specaccum: empy (all zeros) species removed from input
-	data. These confused sd estimation in method = "exact" which gave
-	all NA.
+	* specaccum: empty (all zeros) species removed from input
+	data. These confused method = "exact" which gave sd = NA.
 
+	* permutest.cca: tells whether the test was for all eigenvalues or
+	for the first (this is implicit in anova.cca).
+
 	* vegan-internal.Rd: brief documentation of internal functions for
 	developers with some guidelines of good practices.
 	



More information about the Vegan-commits mailing list