README.md

Harman

The removal of batch effects from datasets using a PCA and constrained optimisation based technique.

Harman is a PCA and constrained optimisation based technique that maximises the removal of batch effects from datasets, with the constraint that the probability of overcorrection (i.e. removing genuine biological signal along with batch noise) is kept to a fraction which is set by the end-user (Oytam et al, 2016; DOI:10.1186/s12859-016-1212-5).

The Harman approach can be generalised to any high dimensional dataset where an experimental factor of interest to be kept is declared and another (usually technical) factor is declared to be removed, with the constraint that the factor to be kept and removed are not completely confounded.

To compile the package on a command line: Windows: Rcmd.exe INSTALL --preclean --no-multiarch --with-keep.source Harman Linux: R CMD INSTALL --preclean --build Harman -l ~/R/library

To set up help files: roxygen2::roxygenize(package.dir = "/path/to/harman")

N.B. * The roxygenize() function should create the NAMESPACE file also.

To use harman:

install.packages("Harman")
library("Harman")
data(OLF)
olf.harman <- harman(olf.data, expt=olf.info$Treatment, batch=olf.info$Batch)
olf.data.corrected <- reconstructData(olf.harman)

Required packages for R CMD Check to pass (Packages used in Vignettes)

if (!requireNamespace("BiocManager", quietly=TRUE))
    install.packages("BiocManager")
BiocManager::install()
BiocManager::install(c("BiocGenerics", "BiocStyle", "knitr", "rmarkdown", "RUnit",
"missMethyl", "RColorBrewer", "bladderbatch", "limma", "minfi", "lumi",
"msmsEDA", "affydata", "minfiData"))

Also required to build vignettes is a version of Texlive for the pdflatex program and a version of the NetCDF library.

Further documentation and the Matlab version is available here: http://www.bioinformatics.csiro.au/harman/



Try the Harman package in your browser

Any scripts or data that you put into this service are public.

Harman documentation built on Nov. 8, 2020, 7:50 p.m.