getSigSegments: Retrieve the significantly gained and lost regions including the corresponding, original probes
Retrieve the significantly gained and lost regions including the corresponding, original probes. A significance level must be selected by the user.
getSigSegments(spm, sigLevels, chromosomes=NULL)
The sample point matrix to be plotted
The significance thresholds to be used
Takes a vector of chromosomes to be plotted. Defaults to all chromosomes.
Returns a sigSegments object containing the chromosome, start position, end position, average KC score and the mode of the KC score in that region of all segments passing the thresholds as set in 'sigLevels'.
Additionally, returns the IDs and indices of the probes and the positions in the sample point matrix within the significant regions.
The results are stored in two separate slots: 'gains' for gains and 'losses' for losses.
Use 'write.table' to save the results to file.
'sigLevels' should contain the significance thresholds in a list with the positive (gains) threshold in the 'pos' element and the negative (losses) threshold in the 'neg' element.
This is the format as returned by 'findSigLevelTrad' and 'findSigLevelFdr'.
References ~put references to the literature/web site here ~
spm1mb <- calcSpm(hsSampleData, hsMirrorLocs)
siglevel1mb <- findSigLevelTrad(hsSampleData, spm1mb, n=3)
sigSegments1mb <- getSigSegments(spm1mb, siglevel1mb)