Adjust for positional and batch effects using ComBat function
In this package, we will carefully compare the effects of positions and batches on your datasets. Then a appropriate correction order will be provided.
You should provide the sample Sentrix information and the methylation levels in matrix file without NA file. The example files shown in the Data file.
library(devtools)
install_github("ChuanJ/JCpackage")
library(posibatch)
dat <-as.matrix(read.csv("beta.csv",header=T,sep=",",row.names=1))
Sentrix <- read.csv("Sentrix.csv",header=T,sep=",")
results <- posibatches(dat, Sentrix, batch=TRUE, par.prior=TRUE, prior.plots=FALSE, mean.only.posi=FALSE, mean.only.batch=FALSE)
write.csv(results,"dataAfterBatch&position.csv")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.