Table 17 R-code for overlapAny matching with extension 2: split at breakpoints
# match |
|---|
> matchedIDs < − matchAnn2Ann(fData(CNdata)[,1], fData(CNdata)[,2], |
fData(CNdata)[,3], fData(GEdata)[,1], |
method = “overlap”) |
> # expand > matchedFeatures < − splitMatchingAtBreakpoints(matchedFeatures, CNdata) |
> # extract ids for object subsetting |
> matchedIDsGE < − lapply(matchedIDs, function(Z){ return(Z[, -2, drop=FALSE]) }) |
> matchedIDsCN < − lapply(matchedIDs, function(Z){ return(Z[, -1, drop=FALSE]) }) |
> # generate matched objects |
> GEdata < − ExpressionSet2weightedSubset(GEdata, matchedIDsGE, 1, 2, 3) |
> CNdata < − cghCall2weightedSubset(CNdata, matchedIDsCN, 1, 2, 3) |