DMRScan_package: DMRScan: An R-package for identification of Differentially...

In christpa/DMRScan: Detection of Differentially Methylated Regions

DMRScan: An R-package for identification of Differentially Metylated Regions

Description

DMRScan: An R-package for identification of Differentially Metylated Regions

Arguments

observations	An object of type GRangesList from makeCpGregions
windowSize	A sequence of windowSizes for the slidingWindow, must be an integer
windowThreshold	Optional argument with corresponding cut-off for each window. Will be estimated if not supplied.
...	Optional arguments to be pased to estimateThreshold, if no grid is specified.

Value

An object of type GRanges with signficantly differentially

Author(s)

Christian Page, [e-mail]

References

Not Published yet (Under revision)

Examples

                                                             
## nProbeoad methylation data from chromosome 22                             
data(DMRScan.methylationData)                                                
## nProbeoad phenotype (end-point for methylation data)                      
data(DMRScan.phenotypes)                                                     
                                                                             
## Test for an association between phenotype and Methylation                 
test.statistics <- apply(DMRScan.methylationData,1,function(x,y)             
  summary(glm(y ~ x, family = binomial(link = "logit")))$coefficients[2,3],  
                                                   y = DMRScan.phenotypes)   
## Set chromosomal position to each test-statistic                           
positions <- data.frame(matrix(as.integer(unlist(strsplit(names(test.statistics), split="chr|[.]"))), ncol = 3, byrow = TRUE))[,-1]
## Set clustering features                                                   
min.cpg <- 4  ## Minimum number of CpGs in a tested cluster                  
## Maxium distance (in base-pairs) within a cluster                          
## before it is broken up into two seperate cluster                          
max.gap <- 750                                                               
                                                                             
## Identify all clusters, and generate a list for each cluster               
regions <- makeCpGregions(observations = test.statistics,                    
                          chr = positions[,1], pos = positions[,2],          
                          maxGap = max.gap, minCpG = min.cpg)                
## Number of CpGs in the slidingWindows, can be either a single number       
## or a sequence of windowSizes                                              
windowSizes <- 3:7                                                           
nCpG        <- sum(sapply(regions, length)) ## Number of CpGs to be tested                  
                                                                             
# Estimate the windowThreshold, based on the number of CpGs and windowSizes  
windowThresholds <- estimateWindowThreshold(nProbe = nCpG,                   
               windowSize = windowSizes, method = "sampling", mcmc = 10000)  
## Run the slidingWindow                                                     
DMRScanResults   <- dmrscan(observations = regions,                          
                            windowSize = windowSizes,                        
                            windowThreshold = windowThresholds)              
## Print the result                                                          
print(DMRScanResults)                                                        

christpa/DMRScan documentation built on Nov. 16, 2024, 9:36 a.m.