Get methylation difference gene
methyDiff(
cpgData,
sampleGroup,
combineMethod = RobustRankAggreg::rhoScores,
missing_value = "knn",
region = "Body",
model = "cpg",
adjust.method = "BH"
)
data.frame of cpg beta value
vector of sample group
method to combine the cpg pvalues
Method to impute missing expression data, one of "zero" and "knn".
region of genes, one of "Body", "TSS1500", "TSS200", "3'UTR", "1stExon", "5'UTR", and "IGR".
if "cpg", step1: calculate difference cpgs; step2: calculate difference genes. if "gene", step1: calculate the methylation level of genes; step2: calculate difference genes.
character string specifying the method used to adjust p-values for multiple testing. See p.adjust for possible values.