topSegments {sesame} | R Documentation |
This is a utility function to show top differential methylated segments. The function takes a coefficient table as input and output the same table ordered by the sigificance of the segments.
topSegments(cf1)
cf1 |
coefficient table of one factor from DMR |
coefficient table ordered by adjusted p-value of segments
sesameDataCache("HM450") # in case not done yet data <- sesameDataGet('HM450.76.TCGA.matched') cf_list = summaryExtractCfList(DML(data$betas, ~type, meta=data$sampleInfo)) cf_list = DMR(data$betas, cf_list$typeTumour) topSegments(cf_list)