now returns p-values

MomX · Sep 28, 2017 · eb600d0 · eb600d0
1 parent 4b0099b
commit eb600d0
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 4 deletions.
diff --git a/R/mult-Manova.R b/R/mult-Manova.R
@@ -316,21 +316,23 @@ MANOVA_PW.PCA <- function(x,
   stars <- as.character(stars)
   nl <- nlevels(fac)
   # bloody ugly
-  stars.tab <- matrix(NA, nrow = nl - 1, ncol = nl,
+  pvalue.tab <- stars.tab <- matrix(NA, nrow = nl - 1, ncol = nl,
                       dimnames = list(levels(fac)[-nl], levels(fac)))
   stars.tab <- as.table(stars.tab)
   k <- 1
   for (i in 1:(nl - 1)) {
     for (j in (i + 1):nl) {
       stars.tab[i, j] <- stars[k]
+      pvalue.tab[i, j] <- res[, 6][k]
       k <- k + 1
     }
   }
   cat("$stars.tab\n")
   print(stars.tab)
   cat("\n$summary (see also $manovas)\n")
   print(res, digits=4)
-  invisible(list(manovas = manovas, summary = res, stars.tab = stars.tab))
+  invisible(list(manovas = manovas, summary = res,
+                 stars.tab = stars.tab, pvalue.tab = pvalue.tab))
 }
 
 # #' @rdname MANOVA_PW

diff --git a/man/MANOVA.Rd b/man/MANOVA.Rd
diff --git a/man/MANOVA_PW.Rd b/man/MANOVA_PW.Rd