EdgeR: how to save results for list of genes ?
2
1
Entering edit mode
Guest User ★ 13k
@guest-user-4897
Last seen 9.6 years ago
Hello, how can I save my list of differentially expressed genes after having my top DE gene visualized with topTags(lrt)? (I followed step by step the example of carcinoma vs matched normal tissue) thanks, anna -- output of sessionInfo(): R version 2.15.1 (2012-06-22) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=French_France.1252 LC_CTYPE=French_France.1252 [3] LC_MONETARY=French_France.1252 LC_NUMERIC=C [5] LC_TIME=French_France.1252 attached base packages: [1] splines stats graphics grDevices utils datasets methods [8] base other attached packages: [1] edgeR_3.0.0 limma_3.14.1 loaded via a namespace (and not attached): [1] tools_2.15.1 -- Sent via the guest posting facility at bioconductor.org.
• 6.8k views
ADD COMMENT
1
Entering edit mode
@gordon-smyth
Last seen 9 hours ago
WEHI, Melbourne, Australia
tab <- topTags(lrt) write.table(tab, file="mygenelist.txt") For more details, type ?write.table Gordon > Date: Wed, 24 Oct 2012 02:32:48 -0700 (PDT) > From: "anna [guest]" <guest at="" bioconductor.org=""> > To: bioconductor at r-project.org, m.nadira at yahoo.fr > Subject: [BioC] EdgeR: how to save results for list of genes ? > > > Hello, > how can I save my list of differentially expressed genes after having my top DE gene visualized with topTags(lrt)? > > (I followed step by step the example of carcinoma vs matched normal tissue) > > thanks, > anna > > -- output of sessionInfo(): > > R version 2.15.1 (2012-06-22) > Platform: i386-pc-mingw32/i386 (32-bit) > > locale: > [1] LC_COLLATE=French_France.1252 LC_CTYPE=French_France.1252 > [3] LC_MONETARY=French_France.1252 LC_NUMERIC=C > [5] LC_TIME=French_France.1252 > > attached base packages: > [1] splines stats graphics grDevices utils datasets methods > [8] base > > other attached packages: > [1] edgeR_3.0.0 limma_3.14.1 > > loaded via a namespace (and not attached): > [1] tools_2.15.1 > ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}
ADD COMMENT
0
Entering edit mode
@nadira-mrabet-5573
Last seen 9.6 years ago
Hello, I mean I want to wite the whole list not only the top genes! I read the help, but I still don't do it correctly. what to do after lrt <- glmLRT(fit) to save data ? thanks, anna tab <- topTags(lrt)   write.table(tab, file="mygenelist.txt") For more details, type ?write.table Gordon >Date: Wed, 24 Oct 2012 02:32:48 -0700 (PDT) >From: "anna [guest]" <guest at="" bioconductor.org=""> >To: bioconductor at r-project.org, m.nadira at yahoo.fr >Subject: [BioC] EdgeR: how to save results for list of genes ? >>> Hello, >how can I save my list of differentially expressed genes after having my top DE gene visualized with topTags(lrt)? >> (I followed step by step the example of carcinoma vs matched normal tissue) >> thanks, >anna [[alternative HTML version deleted]]
ADD COMMENT
0
Entering edit mode
Hi, On Thu, Oct 25, 2012 at 9:55 AM, Nadira MRABET <m.nadira at="" yahoo.fr=""> wrote: > Hello, > I mean I want to wite the whole list not only the top genes! `topTags` takes a paremeter `n` which tells it how many rows of the table to return. You can set that number to `Inf` to get the whole thing. HTH, -steve -- Steve Lianoglou Graduate Student: Computational Systems Biology | Memorial Sloan-Kettering Cancer Center | Weill Medical College of Cornell University Contact Info: http://cbio.mskcc.org/~lianos/contact
ADD REPLY
0
Entering edit mode
Hello, yes now it works like this: > tab <- topTags(lrt, n=Inf) > write.table(tab, file="mygenelist.txt") thank you Steve, anna >________________________________ > De : Steve Lianoglou <mailinglist.honeypot@gmail.com> >À : Nadira MRABET <m.nadira@yahoo.fr> >Cc : "bioconductor@r-project.org" <bioconductor@r-project.org> >Envoyé le : Jeudi 25 octobre 2012 16h08 >Objet : Re: [BioC] EdgeR: how to save results for list of genes ? > >Hi, > >On Thu, Oct 25, 2012 at 9:55 AM, Nadira MRABET <m.nadira@yahoo.fr> wrote: >> Hello, >> I mean I want to wite the whole list not only the top genes! > >`topTags` takes a paremeter `n` which tells it how many rows of the >table to return. You can set that number to `Inf` to get the whole >thing. > >HTH, >-steve > >-- >Steve Lianoglou >Graduate Student: Computational Systems Biology >| Memorial Sloan-Kettering Cancer Center >| Weill Medical College of Cornell University >Contact Info: http://cbio.mskcc.org/~lianos/contact > > > [[alternative HTML version deleted]]
ADD REPLY

Login before adding your answer.

Traffic: 715 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6