Skip to content

Latest commit

 

History

History
31 lines (27 loc) · 967 Bytes

README.md

File metadata and controls

31 lines (27 loc) · 967 Bytes

JCPackage

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.

Install from GitHub:

library(devtools)
install_github("ChuanJ/JCpackage")

library(posibatch)

Usage:

input the methylation levels

dat <-as.matrix(read.csv("beta.csv",header=T,sep=",",row.names=1))

input the Sentrix information

Sentrix <- read.csv("Sentrix.csv",header=T,sep=",")

Correct the positional effect and batch effect

results <- posibatches(dat, Sentrix, batch=TRUE, par.prior=TRUE, prior.plots=FALSE, mean.only.posi=FALSE, mean.only.batch=FALSE)

output the file

write.csv(results,"dataAfterBatch&position.csv")