Title: | Find chromosome regions showing common gains/losses |
---|---|
Description: | This package provides functions to identify genomic regions of interests based on segmented copy number data from multiple samples. |
Authors: | J. Zhang and B. Feng |
Maintainer: | J. Zhang <[email protected]> |
License: | LGPL |
Version: | 1.65.0 |
Built: | 2024-11-29 04:25:03 UTC |
Source: | https://github.com/bioc/cghMCR |
Instantiates a cghMCR object using the parameters passed.
cghMCR(segments, gapAllowed = 500, alteredLow = 0.03, alteredHigh = 0.97, spanLimit = 2e+07, recurrence = 75, thresholdType = c("quantile", "value"))
cghMCR(segments, gapAllowed = 500, alteredLow = 0.03, alteredHigh = 0.97, spanLimit = 2e+07, recurrence = 75, thresholdType = c("quantile", "value"))
segments |
|
gapAllowed |
|
alteredLow |
|
alteredHigh |
|
recurrence |
|
spanLimit |
|
thresholdType |
|
The function is just a constructor of the cghMCR
class
for the instantiation of a cghMCR object
An object of the cghMCR class
The function is a contribution of The Center for Applied Cancer Science of Dana-Farber Cancer Institute
Jianhua Zhang
References on S4 class
data("segData") cghmcr <- cghMCR(segData, gapAllowed = 500, alteredLow = 0.20, alteredHigh = 0.80, recurrence = 50)
data("segData") cghmcr <- cghMCR(segData, gapAllowed = 500, alteredLow = 0.20, alteredHigh = 0.80, recurrence = 50)
Objects of this class provides the functionalities to detecting chromosome regions that show gains or losses across differnet samples
Objects can be created by calls of the form new("cghMCR", ...)
.
A constructor cghMCR
may be used to instantiate object
of this class
DNASeg
:Object of class "data.frame"
containing
segmentation data derived from segmentation analysis using segment
DNAData
:Object of class "data.frame"
containing
raw data derived used for the segmentation analysisfrom segmentation analysis
altered
:Object of class "data.frame"
containing
data for the altered regions
gapAllowed
:Object of class gapAllowed
is an
integer specifying low threshold of base pair number to separate two
adjacent segments, belower which the two segments will be joined as
an altered span
alteredLow
:Object of class alteredLow
is a
positive number between 0 and 1 specifying the lower reshold
percential value. Only segments with values falling below this
threshold are considered as altered span
alteredHigh
:Object of class alteredHigh
is a
positive number between 0 and 1 specifying the upper reshold
percential value. Only segments with values falling over this
threshold are considered as altered span
recurrence
:Object of class recurrence
is an
integer between 1 and 100 that specifies the rate of occurrence for
a gain or loss that are observed across sample. Only gains or losses
with ocurrence rate grater than the threshold values are declared as MCRs
spanLimit
:Object of class spanLimit
is an
integer that defines the leangh of altered spans that can be
considered as locus. It is not of any use at this time
thresholdType
:A character string that can be either "quantile", "value" to indicate the type of the value for recurrence
signature(object = "cghMCR")
: identifies minimum
common regions of gains/losses across samples
The function is a contribution of The Center for Applied Cancer Science of Dana-Farber Cancer Institute
Jianhua Zhang
require("CNTools") data("sampleData") cghmcr <- cghMCR(sampleData[sampleData[, "ID"] %in% sample(unique(sampleData[, "ID"]), 20), ], gapAllowed = 500, alteredLow = 0.20, alteredHigh = 0.80, recurrence = 50)
require("CNTools") data("sampleData") cghmcr <- cghMCR(sampleData[sampleData[, "ID"] %in% sample(unique(sampleData[, "ID"]), 20), ], gapAllowed = 500, alteredLow = 0.20, alteredHigh = 0.80, recurrence = 50)
Methods for function colnames
or rownames
in Package ‘cghMCR’ to extract row or column names from a SGOL object
Extract colnum/row names from SGOL object
This function takes a data frame generated by MCR
and
then append probe ids corresponding to each MCR as a column to the
data frame.
mergeMCRProbes(mcr, rawData)
mergeMCRProbes(mcr, rawData)
mcr |
|
rawData |
|
The mcr
data frame passed must have the first column for
chromosome numbers, the 7th column for the starting positions of
the MCRs, and the 8th column for the ending positions of the MCRs.
A data frame with MCRs and the corresponding probe ids and other data.
The function is a contribution of The Center for Applied Cancer Science of Dana-Farber Cancer Institute
Jianhua Zhang
data("segData") cghmcr <- cghMCR(segData, gapAllowed = 500, alteredLow = 0.20, alteredHigh = 0.80, recurrence = 50) mcrs <- MCR(cghmcr) mcrs <- mergeMCRProbes(mcrs, segData[["data"]])
data("segData") cghmcr <- cghMCR(segData, gapAllowed = 500, alteredLow = 0.20, alteredHigh = 0.80, recurrence = 50) mcrs <- MCR(cghmcr) mcrs <- mergeMCRProbes(mcrs, segData[["data"]])
This function takes a DNAcopy object generated by
running getSegments
and then plots the original data
along with the segments identified.
## S3 method for class 'DNAcopy' plot(x, ..., save = FALSE, layout)
## S3 method for class 'DNAcopy' plot(x, ..., save = FALSE, layout)
x |
|
... |
|
save |
|
layout |
|
This function only works with the DNAcopy returned by
getSegments
. It will work with the DNAcopy returned by
segment of the DNAcopy package when the a column containing the probe
ids is added as the first column of the "data" element of the object.
This function returns invisible() or the name of the plot depending on whether save is set to FALSE of TRUE.
The function is a contribution of The Center for Applied Cancer Science of Dana-Farber Cancer Institute
Jianhua Zhang
The DNAcopy package
if(interactive()){ require("cghMCR") data("segData") plot(segData) }
if(interactive()){ require("cghMCR") data("segData") plot(segData) }
Segments of gains or losses along chromosomes can be calculated based on segmentation data derived from the segment function of the DNAcopy package
Objects can be created by calls of the form new("SGOL", ...)
or
using a constructor SGOL.
gol
:Object of class "matrix"
holding the gain
or loss data for chromosomal segments
threshold
:Object of class "vector"
of length 2
indicating the lower and upper thresholds below/over which data
points will be included in the calculation of SGOL score using the
method defined by method
method
:Object of class "function"
giving the
method used to calculate SGOL scores. Common methods include sum,
median, and mean
signature(object = "SGOL")
: extracts SGOL scores
signature(object = "SGOL")
: gets the name of
the function used to calculate the SGOL scores
signature(x = "SGOL", y = "ANY",...)
: plots the data
signature(object = "SGOL")
: gets the
thredhold used for the calculation
Jianhua Zhang
The SGOL score is a modified version of the GISTIC score published in PNAS 104: 20007-20012
showClass("SGOL") require(CNTools)
showClass("SGOL") require(CNTools)