ANOVA method is applying a specified log normal model to all the features.

## ANOVA dialog

To setup ANOVA model or the alternative Welch's ANOVA (which is used on normally distributed data that violates the assumption of homogeneity of variance), select factors from sample attribute. The factors can be categorical or numeric attribute. Click on a check button to select and click Add factors button to add it to the model (Figure 1).

When more than one factor is selected, click **Add interaction** button to add interaction term of the selected factors.

Once a factor is added to the model, you can specify whether the factor is a random effect (check **Random** check box) or not.

Most factors in an analysis of variance are fixed factors, i.e. the levels of that factor represent all the levels of interest. Examples of fixed factors include gender, treatment, genotype, etc. However, in experiments that are more complex, a factor can be a random effect, meaning the levels of the factor only represent a random sample of all of the levels of interest. Examples of random effects include subject and batch. Consider the example where one factor is type (with levels normal and diseased), and another factor is subject (the subjects selected for the experiment). In this example, “Treatment” is a fixed factor since the levels treated and control represent all conditions of interest. “Subject”, on the other hand, is a random effect since the subjects are only a random sample of all the levels of that factor. When model has both fixed and random effect, it is called a mixed model.

When more than one factor is added to the model, click on the *Cross tabulation* link at the bottom to view the relationship between the factors in a different browser tab (Figure 2).

Once the model is set, click on **Next** button to setup comparisons (contrasts) (Figure 3).

Start by choosing a factor or interaction from the *Factor* drop-down list. The subgroups of the factor or interaction will be displayed in the left panel; click to select one or more level(s) or subgroup name(s) and move them to one of the boxes on the right. The ratio/fold change calculation on the comparison will use the group in the top box as numerator, and the group in the bottom box as the denominator. When multiple levels (groups) are in either numerator or denominator box(es), in *Combine* mode, click on** Add comparison** button to combine all numerator levels and combine all denominator levels in a single comparison in the *Comparison* table below; in Pairwise, click on **Add comparison** button will split all numerator levels and denominator levels into a factorial set of comparisons – in other words, it will add every numerator level paired with every denominator level comparisons to the* Comparison* table . Multiple comparisons from different factors can be added from the specified model.

## ANOVA advanced options

Click on the *Configure* to customize *Advanced options* (Figure 4)

**Low-expression feature** and **Multiple test correction** sections are the same as the matching GSA advanced option, see above GSA advanced options.

**Report option**

Use only reliable estimation results: There are situations when a model estimation procedure does not fail outright, but still encounters some difficulties. In this case, it can even generate p-value and fold change on the comparisons, but they are not reliable, i.e. they can be misleading. Therefore, the default of Use only reliable estimation results is set

**Yes**.Report p-value for effects: If set to

**No**, only the p-value of comparison will be displayed on the report, the p-value of the factors and interaction terms are not shown in the report table. When you choose**Yes**in addition to the comparison’s p-value, type III p-values are displayed for all the non-random terms in the model.Shrinkage to error term variance: by default,

**None**is select, which is lognormal model. Limma-trend and Limma-voom options are lognormal with shrinkage. (Limma-trend is the same as the GSA default option–lognormal with shrinkage). Shrinkage options are recommended for small sample size design, no random effects can be included when performing shrinkage.Report partial correlations: If the model has a numeric factor(s), when choosing

**Yes**, partial correlation coefficient(s) of the numeric factor(s) will be displayed in the result table. When choosing**No**, partial correlation coefficients are not shown.Data has been log transformed with base: showing the current scale of the input data on this task.

## ANOVA report

Since there is only one model for all features, so there is no pie charts design models and response distribution information. The Gene list table format is the same as the GSA report.

## References

- Benjamini, Y., Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing, JRSS, B, 57, 289-300.
- Storey JD. (2003) The positive false discovery rate: A Bayesian interpretation and the q-value. Annals of Statistics, 31: 2013-2035.
- Auer, 2011, A two-stage Poisson model for testing RNA-Seq
- Burnham, Anderson, 2010, Model selection and multimodel inference
- Law C, Voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biology, 2014 15:R29.
- http://cole-trapnell-lab.github.io/cufflinks/cuffdiff/index.html#cuffdiff-output-files
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biology, 2010

## Additional Assistance

If you need additional assistance, please visit our support page to submit a help ticket or find phone numbers for regional support.

Your Rating: | Results: | 41 | rates |

## 1 Comment

## Melissa del Rosario

author: wxw