Added functionality to the gene set testing functions to allow
the user to limit the input CpGs to genomic regions of interest with the
genomic.features
parameter.
Added sig.genes
parameter to the gene set testing functions to
allow the user to add the significant genes overlapping each gene set to the
output.
Added a fix to gsameth to avoid p-values of exactly zero being reported.
Modified the default for coef
in the varFit
function so the
default is all columns of the design matrix used to compute Levene Residuals.
Updated vignette to include examples of how to perform gene set testing following a region level analysis
Bug fix: Fixed bug that resulted in negative value for odds. Bug only occurred in rare cases when CpG universe was smaller than all CpGs on array AND if some large collections had more genes than the universe.
Added fract.counts
logical parameter to enable fractional counting
of CpGs to be turned off in gometh, gsameth, goregion and gsaregion. This
switches off the adjustment for CpGs that map to multiple genes. This is
primarily intended for comparison with earlier versions and we strongly
encourage it to be set to TRUE.
Updated gsameth to take into account that some CpGs map to multiple genes. For a small number of gene families, this previously caused their associated GO categories/gene sets to be erroneously overrepresented and thus highly significant.
Updated gometh to call gsameth to perform statistical testing for GO categories and KEGG pathways; goana and kegga are no longer used.
Replaced use of the AnnotationDbi 'toTable' function with 'select' to speed up execution.
The 'topGSA' function is now used to look at the top ranked results from both gometh and gsameth. The vignette has been modified accordingly.
Updated getMappedEntrezIDs, gometh and gsameth to to speed up execution by taking the array annotation in as an optional argument.
missMethyl now uses the latest IlluminaHumanMethylationEPICanno.ilm10b2.hg19 annotation by default for EPIC arrays.
Added getAdjusted function for extracting RUVm adjusted data for visualisation purposes
Updated vignette to demonstrate use of getAdjusted function
Vignette now includes an example of how to handle cases with RUVm where number of samples is greater than number of Illumina negative controls
modified gene set testing functions gometh and gsameth to accommodate EPIC arrays
New gene set testing function gsameth() added, kegg testing functionality
added to gometh
.
This version of the package contains functions to perform SWAN normalisation, RUV normalisation and differential methylation analysis, differential variability analysis and gene ontology testing for the 450K array.
Added the gometh function to perform gene ontology analysis.
Added functions to perform RUV normalization and differential methylation analysis.
This version of the package has all calculations in matrix formulation and handles nuisance parameters appropriately.
Added new function, contrasts.varFit.
Modified getLeveneResiduals
Modified the varFit function to accept DGEList objects for differential variability testing for RNA-Seq data.
First version of package contains functions to perform SWAN normalisation and differential variability analysis for DNA methylation data from Illumina's Infinium HumanMethylation450 beadchip.