findSegmentsFast: Find Segments in the provided CN data for set parameters

In mknoll/cnAnalysis450k: Copy Number Analysis for 450k Illumina Methylation Arrays

Description

Uses data from minfis getCN() function and normalizes probe-wise against control CN data. Segments are identified with changepoints cpr.var() function (BinSeg) Differes compared to findSegments in: - no plotting - only wilcoxon tests - differences are calculated - only limited set of return values

Usage

findSegmentsFast(data, ctrl, ctrlAll, arrayType = "auto")

Arguments

data	CN data to evaluate
ctrl	CN data of controls, levels to test to (1 mean / median over all ctrl samples)
ctrlAll	CN data of all control samples
arrayType	"auto","450k", "EPIC"; auto -> tries to automatically determine the array type (450k, EPIC)

Value

data containing chr, startCG, endCG, segmentmedian, p-value, samplename

Examples

norm <- minfi::getCN(minfi::preprocessRaw(minfiData::RGsetEx))
ctrlAll <- norm[,5,drop=FALSE]
ctrl <- norm[,4] #ctrl <- apply(ctrlAll, 1, "median")
samples <- norm[,1,drop=FALSE]
findSegmentsFast(samples,ctrl, ctrlAll)[1:4,]

mknoll/cnAnalysis450k documentation built on May 23, 2019, 2:01 a.m.