[BioC] [devteam-bioc] getGmt error

Martin Morgan mtmorgan at fhcrc.org
Fri Sep 6 15:21:31 CEST 2013


On 09/06/2013 06:17 AM, Martin Morgan wrote:
>
>    txt[4694] = sub("NM_009369\t", "", txt[4694])
>    writeLines(txt, "c2.all.v4.0.orig_MODIFIED_.gmt")

maybe a better way is to avoid writing to the intermediate file

   txt = readLines("c2.all.v4.0.orig_MODIFIED_.gmt")
   ## remove duplicate entry in gene set
   txt[4694] = sub("NM_009369\t", "", txt[4694])
   getGmt(textConnection(txt))

and in this way more tightly coupling the transformation with your analysis.

Martin
-- 
Computational Biology / Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N.
PO Box 19024 Seattle, WA 98109

Location: Arnold Building M1 B861
Phone: (206) 667-2793



More information about the Bioconductor mailing list