pmfanova: pmfanova

View source: R/pmf-anova.R

pmfanovaR Documentation



perform principle compponent analysis on ramclustR object dataset, export plots


  ramclustObj = RC, = NULL, = "SpecAbund", = "cmpd",
  subset = c(""), = NULL,
  posthoc = NULL,
  effectsplots = TRUE,
  delim = "-",
  plots = TRUE,
  pcut = 0.05,
  p.adj = "BH",
  which.quan = NULL,
  filter = TRUE,
  summary.statistics = TRUE,
  output.summary = TRUE



ramclustR object to perform ANOVA on

character; optionally assign a name to the anova. valuable if you wish to perform more than one ANOVA on a sample set. If NULL, name will be autogenerated from the and subset information.

character; which dataset (SpecAbund or SpecAbundAve) to perform PCA on.

how should metabolites columns be labelled? one of 'ann' or 'cmpd', typically.


character or integer vector; if character, must be even length. If you wish to perform ANOVA only when your factor called 'treatment' is a 'trt' sample and when 'time' is '3', then you would use i.e. c("treatment", "trt", "time", "3"). vector length must always be even, and with 'factor' followed by 'level'. If an integer vector is provided, only row numbers matching those integers are retained.

character; the model you wish to run. i.e. "Trt * Time" or "Trt * Time + (1|Block)"


character; if you wish to perform post hoc testing, define the posthoc test to use. i.e. if is "Trt * Time" you can use "Trt|Time" to specific 'treatment within time' comparisons.


logical; if TRUE, plot effect plots in pdf format. If plotting fails, can set to 'FALSE' to get tabular results only.


character; what is the character delimiting factors in the sample names. generally '-'


logical; if plotting is any issue set this to FALSE.


numeric; what is the pvalue cutoff for determining significance? Used only for plotting.


character; what p.adjust method should be used. generally 'fdr' or 'bh' (equivalent). see ?p.adjust


logical, TRUE by default. when $cmpd.use slot is present (from function), only cmpds that passed cv filtering are used. If you wish to change that behavior, rerun the function with a really high CV threshold.


logical, TRUE by default. If true, calculate mean and standard deviation for each level as defined by the model.


logical, TRUE by default. If true, saves full model results for all ANOVA as .Rdata file, and exports all summaries to a .txt file.


chracter vector. which factors should be interpreted as quantitative? i.e. c("time", "dose")


This function performs standard or mixed model anova using a ramclusR object as input


returns a ramclustR object. new R object in $pca slot. Optionally, new R object in $AuerGervini slot if npc = "auto".


Corey Broeckling

cbroeckl/ documentation built on Sept. 30, 2022, 5:40 p.m.