R/summary.spermanova.R

In laura20vg/permanova: MANOVA and PERMANOVA

summary.spermanova <- function(perm, latex=FALSE){
  cat(" ###### PERMANOVA : MANOVA based on distances #######\n")
  cat("________________________________________________\n")
  print(perm$call)
  cat("Number of permutations:", perm$nperm, "\n\n")
  SC<- c(round(perm$Inicial$BSS,3), round(perm$Inicial$WSS,2), round(perm$Inicial$TSS,3))
  grados <- c(perm$Inicial$glb, perm$Inicial$glw, perm$Inicial$glt)
  Fexp <- c(round(perm$Inicial$Fexp,3), "", "")
  pvalor <- c(round(perm$pval,10), "", "")
  res <- data.frame(SC, grados, Fexp, pvalor, row.names = c("Between Groups", "Within Groups", "TOTAL"))
  print(res)

  if (latex){
    print(xtable(res, caption="PERMANOVA results"))
  }
}