##NOTE: Comment out in exmaple see vignette for more detailed usage ##
## Load genome bed file ##
data(hg19_genes_bed)
## Load curated DNA methylation bed file ##
#data(methylationdata)
#methylation <- methylation[-which(is.na(methylation[,5])),]
#methylation[, 5] <- replace(methylation[,5],methylation[,5] == 0,
#min(subset(methylation[, 5], methylation[, 5] !=0 ), na.rm=TRUE))
## Create a PvalueAnnotation with defaults for promoter size##
test_annotation<-makePvalueAnnotation(data=hg19_genes, gene_name_col=5)
## Load DNA methylation into PvalueAnnotation ##
#test_annotation <- annotateModification(pvalue_annotation=test_annotation,
#methylation, weight_by=c(promoter="distance", body="distance"), verbose=TRUE,
#mod_corr=FALSE, mod_type="methylation")
## Extract GRanges with modification data ##
#extractModification(pvalue_annotation=test_annotation)
Run the code above in your browser using DataLab