Extract a table of the top-ranked CpGs from a linear model fit after a differential variability analysis.

1 |

`fit` |
list containing a linear model fit produced by |

`coef` |
column number or column name specifying which coefficient of the linear model fit is of interest. It should be the same coefficient that the differential variability testing was performed on. Default is last column of fit object. |

`number` |
maximum number of genes to list. Default is 10. |

`sort` |
logical, default is TRUE. Sorts output according the P-value. FALSE will return results in same order as fit object. |

This function summarises the results of a differential variability analysis performed with `varFit`

. The p-values from the comparison of interest are adjusted using Benjamini and Hochberg's false discovery rate with the function `p.adjust`

. The top ranked CpGs are selected by first ranking the adjusted p-values, then ranking the raw p-values. At this time no other sorting option is catered for.

Produces a dataframe with rows corresponding to the top CpGs and the following columns:

`genelist ` |
one or more columns of annotation for each CpG, if the gene information is available in |

`AvgVar ` |
average of the absolute or squared Levene residuals across all samples |

`DiffVar ` |
estimate of the difference in the Levene residuals corresponding to the comparison of interest |

`t ` |
moderated t-statistic |

`P.Value ` |
raw p-value |

`Adj.P.Value ` |
adjusted p-value |

Belinda Phipson

Phipson, B., and Oshlack, A. (2014). A method for detecting differential variability in methylation data shows CpG islands are highly variably methylated in cancers. *Genome Biology*, **15**:465.

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. *Journal of the Royal Statistical Society Series*, B, **57**, 289-300.

`varFit`

, `p.adjust`

1 2 3 4 5 6 7 8 9 10 11 12 13 | ```
# Randomly generate data for a 2 group problem with 100 CpG sites and 5 arrays in each group.
y<-matrix(rnorm(1000),ncol=10)
group<-factor(rep(c(1,2),each=5))
design<-model.matrix(~group)
# Fit linear model for differential variability
vfit<-varFit(y,design)
# Look at top table of results
topVar(vfit,coef=2)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.