[BioC] Help on cghMCR units of gapAllowed

Nathalie Conte nac at sanger.ac.uk
Thu Dec 8 16:51:29 CET 2011


HI,
 I just want a precision on cghMCR package for those familiar with it.
The function cghMCR will have an argument about gap allowed which will 
fuse segment of DNA below this size. What is the scale is it bp or Kb?In 
the manual, it seems quite clear that it is pb but in one example, it 
seems to be kb...
|segments| 	|segments| is a data frame extracted from the "output" 
element of the object returned by segment of the package DNAcopy or 
|getSegments 
<http://rss.acs.unt.edu/Rdoc/library/cghMCR/html/cghMCR-class.html>|
|gapAllowed| 	|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| 	|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| 	|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| 	|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/losses with ocurrence rate grater than the 
threshold values are declared as MCRs

example given by the manual, specify 500kb and this is written like 
this... should be 500 000 if the unit is bp

  require("cghMCR")
  data("sampleData")
  segments <- getSegments(sampleData)
  cghmcr <- cghMCR(segments, gapAllowed = 500, alteredLow = 0.20,
                   alteredHigh = 0.80, recurrence = 50)

I am a little confused , could you please help
thanks






-- 
 The Wellcome Trust Sanger Institute is operated by Genome Research 
 Limited, a charity registered in England with number 1021457 and a 
 company registered in England with number 2742969, whose registered 
 office is 215 Euston Road, London, NW1 2BE.



More information about the Bioconductor mailing list