Time course experiment....
2
0
Entering edit mode
viritha kaza ▴ 580
@viritha-kaza-4318
Last seen 9.6 years ago
Naomi Altman <naomi at="" ...=""> writes: > > I did not get the original posting, but doesn't Sohail just need > "TopTable" for this? > > --Naomi Hi Naomi, I had a similiar issue and used TopTable but I am getting only first 10 but how should I write the code to get selected probesets based on an adjusted p-value less than or equal to 0.05. code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) topTable(fit2,coef=2,adjust="BH") I have 3 control and 3 treated samples.The eset has the probenames and the expression matrix. So my question is how do I get the data of all the probesets which pass this criteria? Thanks, Viritha > > At 08:51 AM 1/2/2006, Gordon Smyth wrote: > >Dear Sohail, > > > >Well, there are lots of ways to generate such a table. Perhaps the simplest is > > > > fitsel <- fit2[sel.dif, ] > > as.data.frame( fitsel ) > > > >Best wishes > >Gordon > > > > >Date: Tue, 20 Dec 2005 14:03:43 -0500 > > >From: "Khan, Sohail" <khan at="" ...=""> > > >Subject: [BioC] Time course experiment.... > > >To: <bioconductor at="" ...=""> > > > > > >Dear List, > > > > > >I have performed a time course analysis using limma, as described in > > >"Bioinformatics and Computational Biology Solutions ........". > > >How can I get a list of differentially expressed genes? I've tried > > >the code below: > > >sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr") <0.05 > > >This produces a logical vector, right?. I would like a table of > > >differentially expressed genes with p vales etc. Sorry, if I missed > > >this in the limma user's guide. Thanks for any suggestions. > > > > > > > > >Sohail Khan > > >Scientific Programmer > > >COLD SPRING HARBOR LABORATORY > > >Genome Research Center > > >500 Sunnyside Boulevard > > >Woodbury, NY 11797 > > >(516)422-4076 > > > >_______________________________________________ > >Bioconductor mailing list > >Bioconductor at ... > >https://stat.ethz.ch/mailman/listinfo/bioconductor > > Naomi S. Altman 814-865-3791 (voice) > Associate Professor > Dept. of Statistics 814-863-7114 (fax) > Penn State University 814-865-1348 (Statistics) > University Park, PA 16802-2111 >
limma limma • 1.3k views
ADD COMMENT
0
Entering edit mode
@manca-marco-path-4295
Last seen 9.6 years ago
Hi Viritha, just do code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) top.list<-topTable(fit2,coef=2,adjust="BH") your.results <- top.list[top.list$adj.P.Val < 0.05, ] ...it should do the trick... All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX Maastricht E-mail: m.manca at maastrichtuniversity.nl Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confidential and solely for the use of the intended recipient. It may contain material protected by privacy or attorney-client privilege. If you are not the intended recipient or the person responsible for delivering to the intended recipient, be advised that you have received this email in error and that any use is STRICTLY PROHIBITED. If you have received this email in error please notify us by telephone on +31626441205 Dr Marco MANCA ********************************************************************** *********************************************** ________________________________________ Da: bioconductor-bounces at r-project.org [bioconductor-bounces at r-project.org] per conto di viritha [viritha.k at gmail.com] Inviato: marted? 8 febbraio 2011 18.14 A: bioconductor at stat.math.ethz.ch Oggetto: Re: [BioC] Time course experiment.... Naomi Altman <naomi at="" ...=""> writes: > > I did not get the original posting, but doesn't Sohail just need > "TopTable" for this? > > --Naomi Hi Naomi, I had a similiar issue and used TopTable but I am getting only first 10 but how should I write the code to get selected probesets based on an adjusted p-value less than or equal to 0.05. code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) topTable(fit2,coef=2,adjust="BH") I have 3 control and 3 treated samples.The eset has the probenames and the expression matrix. So my question is how do I get the data of all the probesets which pass this criteria? Thanks, Viritha > > At 08:51 AM 1/2/2006, Gordon Smyth wrote: > >Dear Sohail, > > > >Well, there are lots of ways to generate such a table. Perhaps the simplest is > > > > fitsel <- fit2[sel.dif, ] > > as.data.frame( fitsel ) > > > >Best wishes > >Gordon > > > > >Date: Tue, 20 Dec 2005 14:03:43 -0500 > > >From: "Khan, Sohail" <khan at="" ...=""> > > >Subject: [BioC] Time course experiment.... > > >To: <bioconductor at="" ...=""> > > > > > >Dear List, > > > > > >I have performed a time course analysis using limma, as described in > > >"Bioinformatics and Computational Biology Solutions ........". > > >How can I get a list of differentially expressed genes? I've tried > > >the code below: > > >sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr") <0.05 > > >This produces a logical vector, right?. I would like a table of > > >differentially expressed genes with p vales etc. Sorry, if I missed > > >this in the limma user's guide. Thanks for any suggestions. > > > > > > > > >Sohail Khan > > >Scientific Programmer > > >COLD SPRING HARBOR LABORATORY > > >Genome Research Center > > >500 Sunnyside Boulevard > > >Woodbury, NY 11797 > > >(516)422-4076 > > > >_______________________________________________ > >Bioconductor mailing list > >Bioconductor at ... > >https://stat.ethz.ch/mailman/listinfo/bioconductor > > Naomi S. Altman 814-865-3791 (voice) > Associate Professor > Dept. of Statistics 814-863-7114 (fax) > Penn State University 814-865-1348 (Statistics) > University Park, PA 16802-2111 > _______________________________________________ Bioconductor mailing list Bioconductor at r-project.org https://stat.ethz.ch/mailman/listinfo/bioconductor Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD COMMENT
0
Entering edit mode
Hi Marco, Thanks for your reply. But topTable by default only gives top 10.So how do I know the number in topTable fuction which satisfies the criteria of 0.05 that I need to give. waiting for your reply... Thanks, Viritha On Tue, Feb 8, 2011 at 12:35 PM, Manca Marco (PATH) < m.manca@maastrichtuniversity.nl> wrote: > Hi Viritha, > > just do > > code: > design <- model.matrix(~factor(rep(1:2, each=3))) > fit <- lmFit(eset, design) > fit2 <- eBayes(fit) > top.list<-topTable(fit2,coef=2,adjust="BH") > your.results <- top.list[top.list$adj.P.Val < 0.05, ] > > ...it should do the trick... > > All the best, Marco > > > -- > Marco Manca, MD > University of Maastricht > Faculty of Health, Medicine and Life Sciences (FHML) > Cardiovascular Research Institute (CARIM) > > Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) > Visiting address: Experimental Vascular Pathology group, Dept of Pathology > - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX > Maastricht > > E-mail: m.manca@maastrichtuniversity.nl > Office telephone: +31(0)433874633 > Personal mobile: +31(0)626441205 > Twitter: @markomanka > > > > ******************************************************************** ************************************************* > > This email and any files transmitted with it are confidential and solely > for the use of the intended recipient. > > It may contain material protected by privacy or attorney-client privilege. > If you are not the intended recipient or the person responsible for > > delivering to the intended recipient, be advised that you have received > this email in error and that any use is STRICTLY PROHIBITED. > > If you have received this email in error please notify us by telephone on > +31626441205 Dr Marco MANCA > > > ******************************************************************** ************************************************* > ________________________________________ > Da: bioconductor-bounces@r-project.org [bioconductor- bounces@r-project.org] > per conto di viritha [viritha.k@gmail.com] > Inviato: martedì 8 febbraio 2011 18.14 > A: bioconductor@stat.math.ethz.ch > Oggetto: Re: [BioC] Time course experiment.... > > Naomi Altman <naomi@...> writes: > > > > > I did not get the original posting, but doesn't Sohail just need > > "TopTable" for this? > > > > --Naomi > > Hi Naomi, > I had a similiar issue and used TopTable but I am getting only first 10 but > how > should I write the code to get selected probesets based on an adjusted > p-value > less than or equal to 0.05. > > code: > design <- model.matrix(~factor(rep(1:2, each=3))) > fit <- lmFit(eset, design) > fit2 <- eBayes(fit) > topTable(fit2,coef=2,adjust="BH") > > I have 3 control and 3 treated samples.The eset has the probenames and the > expression matrix. > So my question is how do I get the data of all the probesets which pass > this > criteria? > Thanks, > Viritha > > > > > At 08:51 AM 1/2/2006, Gordon Smyth wrote: > > >Dear Sohail, > > > > > >Well, there are lots of ways to generate such a table. Perhaps the > simplest > is > > > > > > fitsel <- fit2[sel.dif, ] > > > as.data.frame( fitsel ) > > > > > >Best wishes > > >Gordon > > > > > > >Date: Tue, 20 Dec 2005 14:03:43 -0500 > > > >From: "Khan, Sohail" <khan@...> > > > >Subject: [BioC] Time course experiment.... > > > >To: <bioconductor@...> > > > > > > > >Dear List, > > > > > > > >I have performed a time course analysis using limma, as described in > > > >"Bioinformatics and Computational Biology Solutions ........". > > > >How can I get a list of differentially expressed genes? I've tried > > > >the code below: > > > >sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr") <0.05 > > > >This produces a logical vector, right?. I would like a table of > > > >differentially expressed genes with p vales etc. Sorry, if I missed > > > >this in the limma user's guide. Thanks for any suggestions. > > > > > > > > > > > >Sohail Khan > > > >Scientific Programmer > > > >COLD SPRING HARBOR LABORATORY > > > >Genome Research Center > > > >500 Sunnyside Boulevard > > > >Woodbury, NY 11797 > > > >(516)422-4076 > > > > > >_______________________________________________ > > >Bioconductor mailing list > > >Bioconductor@... > > >https://stat.ethz.ch/mailman/listinfo/bioconductor > > > > Naomi S. Altman 814-865-3791 (voice) > > Associate Professor > > Dept. of Statistics 814-863-7114 (fax) > > Penn State University 814-865-1348 (Statistics) > > University Park, PA 16802-2111 > > > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > [[alternative HTML version deleted]]
ADD REPLY
0
Entering edit mode
Hi Viritha, On 02/08/2011 07:01 PM, viritha kaza wrote: > Hi Marco, > Thanks for your reply. > But topTable by default only gives top 10.So how do I know the number in > topTable fuction which satisfies the criteria of 0.05 that I need to give. > waiting for your reply... > Thanks, > Viritha the help page of topTable says that you can specify the number of returned genes by using the 'number' parameter. You could set this to the maximum (ie the number of genes you tested) and then filter. HTH. James. > > On Tue, Feb 8, 2011 at 12:35 PM, Manca Marco (PATH)< > m.manca at maastrichtuniversity.nl> wrote: > >> Hi Viritha, >> >> just do >> >> code: >> design<- model.matrix(~factor(rep(1:2, each=3))) >> fit<- lmFit(eset, design) >> fit2<- eBayes(fit) >> top.list<-topTable(fit2,coef=2,adjust="BH") >> your.results<- top.list[top.list$adj.P.Val< 0.05, ] >> >> ...it should do the trick... >> >> All the best, Marco >> >> >> -- >> Marco Manca, MD >> University of Maastricht >> Faculty of Health, Medicine and Life Sciences (FHML) >> Cardiovascular Research Institute (CARIM) >> >> Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) >> Visiting address: Experimental Vascular Pathology group, Dept of Pathology >> - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX >> Maastricht >> >> E-mail: m.manca at maastrichtuniversity.nl >> Office telephone: +31(0)433874633 >> Personal mobile: +31(0)626441205 >> Twitter: @markomanka >> >> >> >> ******************************************************************* ************************************************** >> >> This email and any files transmitted with it are confidential and solely >> for the use of the intended recipient. >> >> It may contain material protected by privacy or attorney-client privilege. >> If you are not the intended recipient or the person responsible for >> >> delivering to the intended recipient, be advised that you have received >> this email in error and that any use is STRICTLY PROHIBITED. >> >> If you have received this email in error please notify us by telephone on >> +31626441205 Dr Marco MANCA >> >> >> ******************************************************************* ************************************************** >> ________________________________________ >> Da: bioconductor-bounces at r-project.org [bioconductor-bounces at r-project.org] >> per conto di viritha [viritha.k at gmail.com] >> Inviato: marted? 8 febbraio 2011 18.14 >> A: bioconductor at stat.math.ethz.ch >> Oggetto: Re: [BioC] Time course experiment.... >> >> Naomi Altman<naomi at="" ...=""> writes: >> >>> >>> I did not get the original posting, but doesn't Sohail just need >>> "TopTable" for this? >>> >>> --Naomi >> >> Hi Naomi, >> I had a similiar issue and used TopTable but I am getting only first 10 but >> how >> should I write the code to get selected probesets based on an adjusted >> p-value >> less than or equal to 0.05. >> >> code: >> design<- model.matrix(~factor(rep(1:2, each=3))) >> fit<- lmFit(eset, design) >> fit2<- eBayes(fit) >> topTable(fit2,coef=2,adjust="BH") >> >> I have 3 control and 3 treated samples.The eset has the probenames and the >> expression matrix. >> So my question is how do I get the data of all the probesets which pass >> this >> criteria? >> Thanks, >> Viritha >> >>> >>> At 08:51 AM 1/2/2006, Gordon Smyth wrote: >>>> Dear Sohail, >>>> >>>> Well, there are lots of ways to generate such a table. Perhaps the >> simplest >> is >>>> >>>> fitsel<- fit2[sel.dif, ] >>>> as.data.frame( fitsel ) >>>> >>>> Best wishes >>>> Gordon >>>> >>>>> Date: Tue, 20 Dec 2005 14:03:43 -0500 >>>>> From: "Khan, Sohail"<khan at="" ...=""> >>>>> Subject: [BioC] Time course experiment.... >>>>> To:<bioconductor at="" ...=""> >>>>> >>>>> Dear List, >>>>> >>>>> I have performed a time course analysis using limma, as described in >>>>> "Bioinformatics and Computational Biology Solutions ........". >>>>> How can I get a list of differentially expressed genes? I've tried >>>>> the code below: >>>>> sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr")<0.05 >>>>> This produces a logical vector, right?. I would like a table of >>>>> differentially expressed genes with p vales etc. Sorry, if I missed >>>>> this in the limma user's guide. Thanks for any suggestions. >>>>> >>>>> >>>>> Sohail Khan >>>>> Scientific Programmer >>>>> COLD SPRING HARBOR LABORATORY >>>>> Genome Research Center >>>>> 500 Sunnyside Boulevard >>>>> Woodbury, NY 11797 >>>>> (516)422-4076 >>>> >>>> _______________________________________________ >>>> Bioconductor mailing list >>>> Bioconductor at ... >>>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>> >>> Naomi S. Altman 814-865-3791 (voice) >>> Associate Professor >>> Dept. of Statistics 814-863-7114 (fax) >>> Penn State University 814-865-1348 (Statistics) >>> University Park, PA 16802-2111 >>> >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > > [[alternative HTML version deleted]] > > > > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD REPLY
0
Entering edit mode
Hi James, Thanks !!!!!!! Thats a good idea. -Viritha On Tue, Feb 8, 2011 at 2:30 PM, James F. Reid <james.reid@ifom-ieo- campus.it=""> wrote: > Hi Viritha, > > > On 02/08/2011 07:01 PM, viritha kaza wrote: > >> Hi Marco, >> Thanks for your reply. >> But topTable by default only gives top 10.So how do I know the number in >> topTable fuction which satisfies the criteria of 0.05 that I need to give. >> waiting for your reply... >> Thanks, >> Viritha >> > > the help page of topTable says that you can specify the number of returned > genes by using the 'number' parameter. > You could set this to the maximum (ie the number of genes you tested) and > then filter. > > HTH. > James. > > >> On Tue, Feb 8, 2011 at 12:35 PM, Manca Marco (PATH)< >> m.manca@maastrichtuniversity.nl> wrote: >> >> Hi Viritha, >>> >>> just do >>> >>> code: >>> design<- model.matrix(~factor(rep(1:2, each=3))) >>> fit<- lmFit(eset, design) >>> fit2<- eBayes(fit) >>> top.list<-topTable(fit2,coef=2,adjust="BH") >>> your.results<- top.list[top.list$adj.P.Val< 0.05, ] >>> >>> ...it should do the trick... >>> >>> All the best, Marco >>> >>> >>> -- >>> Marco Manca, MD >>> University of Maastricht >>> Faculty of Health, Medicine and Life Sciences (FHML) >>> Cardiovascular Research Institute (CARIM) >>> >>> Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) >>> Visiting address: Experimental Vascular Pathology group, Dept of >>> Pathology >>> - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 >>> HX >>> Maastricht >>> >>> E-mail: m.manca@maastrichtuniversity.nl >>> Office telephone: +31(0)433874633 >>> Personal mobile: +31(0)626441205 >>> Twitter: @markomanka >>> >>> >>> >>> >>> ****************************************************************** *************************************************** >>> >>> This email and any files transmitted with it are confidential and solely >>> for the use of the intended recipient. >>> >>> It may contain material protected by privacy or attorney-client >>> privilege. >>> If you are not the intended recipient or the person responsible for >>> >>> delivering to the intended recipient, be advised that you have received >>> this email in error and that any use is STRICTLY PROHIBITED. >>> >>> If you have received this email in error please notify us by telephone on >>> +31626441205 Dr Marco MANCA >>> >>> >>> >>> ****************************************************************** *************************************************** >>> ________________________________________ >>> Da: bioconductor-bounces@r-project.org [ >>> bioconductor-bounces@r-project.org] >>> per conto di viritha [viritha.k@gmail.com] >>> Inviato: martedì 8 febbraio 2011 18.14 >>> A: bioconductor@stat.math.ethz.ch >>> Oggetto: Re: [BioC] Time course experiment.... >>> >>> Naomi Altman<naomi@...> writes: >>> >>> >>>> I did not get the original posting, but doesn't Sohail just need >>>> "TopTable" for this? >>>> >>>> --Naomi >>>> >>> >>> Hi Naomi, >>> I had a similiar issue and used TopTable but I am getting only first 10 >>> but >>> how >>> should I write the code to get selected probesets based on an adjusted >>> p-value >>> less than or equal to 0.05. >>> >>> code: >>> design<- model.matrix(~factor(rep(1:2, each=3))) >>> fit<- lmFit(eset, design) >>> fit2<- eBayes(fit) >>> topTable(fit2,coef=2,adjust="BH") >>> >>> I have 3 control and 3 treated samples.The eset has the probenames and >>> the >>> expression matrix. >>> So my question is how do I get the data of all the probesets which pass >>> this >>> criteria? >>> Thanks, >>> Viritha >>> >>> >>>> At 08:51 AM 1/2/2006, Gordon Smyth wrote: >>>> >>>>> Dear Sohail, >>>>> >>>>> Well, there are lots of ways to generate such a table. Perhaps the >>>>> >>>> simplest >>> is >>> >>>> >>>>> fitsel<- fit2[sel.dif, ] >>>>> as.data.frame( fitsel ) >>>>> >>>>> Best wishes >>>>> Gordon >>>>> >>>>> Date: Tue, 20 Dec 2005 14:03:43 -0500 >>>>>> From: "Khan, Sohail"<khan@...> >>>>>> Subject: [BioC] Time course experiment.... >>>>>> To:<bioconductor@...> >>>>>> >>>>>> Dear List, >>>>>> >>>>>> I have performed a time course analysis using limma, as described in >>>>>> "Bioinformatics and Computational Biology Solutions ........". >>>>>> How can I get a list of differentially expressed genes? I've tried >>>>>> the code below: >>>>>> sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr")<0.05 >>>>>> This produces a logical vector, right?. I would like a table of >>>>>> differentially expressed genes with p vales etc. Sorry, if I missed >>>>>> this in the limma user's guide. Thanks for any suggestions. >>>>>> >>>>>> >>>>>> Sohail Khan >>>>>> Scientific Programmer >>>>>> COLD SPRING HARBOR LABORATORY >>>>>> Genome Research Center >>>>>> 500 Sunnyside Boulevard >>>>>> Woodbury, NY 11797 >>>>>> (516)422-4076 >>>>>> >>>>> >>>>> _______________________________________________ >>>>> Bioconductor mailing list >>>>> Bioconductor@... >>>>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>>>> >>>> >>>> Naomi S. Altman 814-865-3791 (voice) >>>> Associate Professor >>>> Dept. of Statistics 814-863-7114 (fax) >>>> Penn State University 814-865-1348 (Statistics) >>>> University Park, PA 16802-2111 >>>> >>>> >>> _______________________________________________ >>> Bioconductor mailing list >>> Bioconductor@r-project.org >>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>> Search the archives: >>> http://news.gmane.org/gmane.science.biology.informatics.conductor >>> >>> >> [[alternative HTML version deleted]] >> >> >> >> >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor@r-project.org >> https://stat.ethz.ch/mailman/listinfo/bioconductor >> Search the archives: >> http://news.gmane.org/gmane.science.biology.informatics.conductor >> > [[alternative HTML version deleted]]
ADD REPLY
0
Entering edit mode
Sorry, I had overlooked that part of the question... I thought you just wanted help in setting the threshold so I just added that line to your code. As James has already mentioned to have the full list of genes in your results you have to set something like: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) # If you want the first 10000 results... but usually I put "number=THE_NUMBER_OF_GENES_ON_MY_CHIP" top.list<-topTable(fit2,number=10000,coef=2,adjust="BH") your.results <- top.list[top.list$adj.P.Val < 0.05, ] ...I can only say to my defense "you only see what you look for" =P I should have given a more complete answer since the beginning. All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX Maastricht E-mail: m.manca at maastrichtuniversity.nl Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confidential and solely for the use of the intended recipient. It may contain material protected by privacy or attorney-client privilege. If you are not the intended recipient or the person responsible for delivering to the intended recipient, be advised that you have received this email in error and that any use is STRICTLY PROHIBITED. If you have received this email in error please notify us by telephone on +31626441205 Dr Marco MANCA ********************************************************************** *********************************************** ________________________________________ Da: viritha kaza [viritha.k at gmail.com] Inviato: marted? 8 febbraio 2011 19.01 A: Manca Marco (PATH) Cc: bioconductor at stat.math.ethz.ch Oggetto: Re: [BioC] Time course experiment.... Hi Marco, Thanks for your reply. But topTable by default only gives top 10.So how do I know the number in topTable fuction which satisfies the criteria of 0.05 that I need to give. waiting for your reply... Thanks, Viritha On Tue, Feb 8, 2011 at 12:35 PM, Manca Marco (PATH) <m.manca at="" maastrichtuniversity.nl<mailto:m.manca="" at="" maastrichtuniversity.nl="">> wrote: Hi Viritha, just do code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) top.list<-topTable(fit2,coef=2,adjust="BH") your.results <- top.list[top.list$adj.P.Val < 0.05, ] ...it should do the trick... All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX Maastricht E-mail: m.manca at maastrichtuniversity.nl<mailto:m.manca at="" maastrichtuniversity.nl=""> Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confidential and solely for the use of the intended recipient. It may contain material protected by privacy or attorney-client privilege. If you are not the intended recipient or the person responsible for delivering to the intended recipient, be advised that you have received this email in error and that any use is STRICTLY PROHIBITED. If you have received this email in error please notify us by telephone on +31626441205 Dr Marco MANCA ********************************************************************** *********************************************** ________________________________________ Da: bioconductor-bounces at r-project.org<mailto:bioconductor-bounces at="" r-project.org=""> [bioconductor-bounces at r-project.org<mailto :bioconductor-bounces="" at="" r-project.org="">] per conto di viritha [viritha.k at gmail.com<mailto:viritha.k at="" gmail.com="">] Inviato: marted? 8 febbraio 2011 18.14 A: bioconductor at stat.math.ethz.ch<mailto:bioconductor at="" stat.math.ethz.ch=""> Oggetto: Re: [BioC] Time course experiment.... Naomi Altman <naomi at="" ...=""> writes: > > I did not get the original posting, but doesn't Sohail just need > "TopTable" for this? > > --Naomi Hi Naomi, I had a similiar issue and used TopTable but I am getting only first 10 but how should I write the code to get selected probesets based on an adjusted p-value less than or equal to 0.05. code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) topTable(fit2,coef=2,adjust="BH") I have 3 control and 3 treated samples.The eset has the probenames and the expression matrix. So my question is how do I get the data of all the probesets which pass this criteria? Thanks, Viritha > > At 08:51 AM 1/2/2006, Gordon Smyth wrote: > >Dear Sohail, > > > >Well, there are lots of ways to generate such a table. Perhaps the simplest is > > > > fitsel <- fit2[sel.dif, ] > > as.data.frame( fitsel ) > > > >Best wishes > >Gordon > > > > >Date: Tue, 20 Dec 2005 14:03:43 -0500 > > >From: "Khan, Sohail" <khan at="" ...=""> > > >Subject: [BioC] Time course experiment.... > > >To: <bioconductor at="" ...=""> > > > > > >Dear List, > > > > > >I have performed a time course analysis using limma, as described in > > >"Bioinformatics and Computational Biology Solutions ........". > > >How can I get a list of differentially expressed genes? I've tried > > >the code below: > > >sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr") <0.05 > > >This produces a logical vector, right?. I would like a table of > > >differentially expressed genes with p vales etc. Sorry, if I missed > > >this in the limma user's guide. Thanks for any suggestions. > > > > > > > > >Sohail Khan > > >Scientific Programmer > > >COLD SPRING HARBOR LABORATORY > > >Genome Research Center > > >500 Sunnyside Boulevard > > >Woodbury, NY 11797 > > >(516)422-4076 > > > >_______________________________________________ > >Bioconductor mailing list > >Bioconductor at ... > >https://stat.ethz.ch/mailman/listinfo/bioconductor > > Naomi S. Altman 814-865-3791 (voice) > Associate Professor > Dept. of Statistics 814-863-7114 (fax) > Penn State University 814-865-1348 (Statistics) > University Park, PA 16802-2111 > _______________________________________________ Bioconductor mailing list Bioconductor at r-project.org<mailto:bioconductor at="" r-project.org=""> https://stat.ethz.ch/mailman/listinfo/bioconductor Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD REPLY
0
Entering edit mode
Hi Marco, Thanks... That happens. I wanted to include one more condition on the top.table list of logFC to consider only those which are greater than or equal to1 and less than or equal to -1. I just tried to expand the code > your.results1<-your.results[your.results$logFC < -1,] > your.results2<-your.results1[your.results1$logFC > 1,] > write.table(your.results2,"Panc1_final_list.txt",sep='\t') but just got a blank txt file with headings.I also tried using & and | operators, but got errors. I hope you could help on this. Thanks, Viritha On Tue, Feb 8, 2011 at 3:32 PM, Manca Marco (PATH) < m.manca@maastrichtuniversity.nl> wrote: > Sorry, I had overlooked that part of the question... I thought you just > wanted help in setting the threshold so I just added that line to your code. > > As James has already mentioned to have the full list of genes in your > results you have to set something like: > > design <- model.matrix(~factor(rep(1:2, each=3))) > fit <- lmFit(eset, design) > fit2 <- eBayes(fit) > # If you want the first 10000 results... but usually I put > "number=THE_NUMBER_OF_GENES_ON_MY_CHIP" > top.list<-topTable(fit2,number=10000,coef=2,adjust="BH") > your.results <- top.list[top.list$adj.P.Val < 0.05, ] > > ...I can only say to my defense "you only see what you look for" =P I > should have given a more complete answer since the beginning. > > All the best, Marco > > > -- > Marco Manca, MD > University of Maastricht > Faculty of Health, Medicine and Life Sciences (FHML) > Cardiovascular Research Institute (CARIM) > > Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) > Visiting address: Experimental Vascular Pathology group, Dept of Pathology > - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX > Maastricht > > E-mail: m.manca@maastrichtuniversity.nl > Office telephone: +31(0)433874633 > Personal mobile: +31(0)626441205 > Twitter: @markomanka > > > > ******************************************************************** ************************************************* > > This email and any files transmitted with it are confidential and solely > for the use of the intended recipient. > > It may contain material protected by privacy or attorney-client privilege. > If you are not the intended recipient or the person responsible for > > delivering to the intended recipient, be advised that you have received > this email in error and that any use is STRICTLY PROHIBITED. > > If you have received this email in error please notify us by telephone on > +31626441205 Dr Marco MANCA > > > ******************************************************************** ************************************************* > ________________________________________ > Da: viritha kaza [viritha.k@gmail.com] > Inviato: martedì 8 febbraio 2011 19.01 > A: Manca Marco (PATH) > Cc: bioconductor@stat.math.ethz.ch > Oggetto: Re: [BioC] Time course experiment.... > > Hi Marco, > Thanks for your reply. > But topTable by default only gives top 10.So how do I know the number in > topTable fuction which satisfies the criteria of 0.05 that I need to give. > waiting for your reply... > Thanks, > Viritha > > On Tue, Feb 8, 2011 at 12:35 PM, Manca Marco (PATH) < > m.manca@maastrichtuniversity.nl<mailto:m.manca@maastrichtuniversity. nl="">> > wrote: > Hi Viritha, > > just do > > code: > design <- model.matrix(~factor(rep(1:2, each=3))) > fit <- lmFit(eset, design) > fit2 <- eBayes(fit) > top.list<-topTable(fit2,coef=2,adjust="BH") > your.results <- top.list[top.list$adj.P.Val < 0.05, ] > > ...it should do the trick... > > All the best, Marco > > > -- > Marco Manca, MD > University of Maastricht > Faculty of Health, Medicine and Life Sciences (FHML) > Cardiovascular Research Institute (CARIM) > > Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) > Visiting address: Experimental Vascular Pathology group, Dept of Pathology > - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX > Maastricht > > E-mail: m.manca@maastrichtuniversity.nl<mailto:> m.manca@maastrichtuniversity.nl> > Office telephone: +31(0)433874633 > Personal mobile: +31(0)626441205 > Twitter: @markomanka > > > > ******************************************************************** ************************************************* > > This email and any files transmitted with it are confidential and solely > for the use of the intended recipient. > > It may contain material protected by privacy or attorney-client privilege. > If you are not the intended recipient or the person responsible for > > delivering to the intended recipient, be advised that you have received > this email in error and that any use is STRICTLY PROHIBITED. > > If you have received this email in error please notify us by telephone on > +31626441205 Dr Marco MANCA > > > ******************************************************************** ************************************************* > ________________________________________ > Da: bioconductor-bounces@r-project.org<mailto:> bioconductor-bounces@r-project.org> [bioconductor- bounces@r-project.org > <mailto:bioconductor-bounces@r-project.org>] per conto di viritha [ > viritha.k@gmail.com<mailto:viritha.k@gmail.com>] > Inviato: martedì 8 febbraio 2011 18.14 > A: bioconductor@stat.math.ethz.ch<mailto:bioconductor@stat.math.ethz.ch> > Oggetto: Re: [BioC] Time course experiment.... > > Naomi Altman <naomi@...> writes: > > > > > I did not get the original posting, but doesn't Sohail just need > > "TopTable" for this? > > > > --Naomi > > Hi Naomi, > I had a similiar issue and used TopTable but I am getting only first 10 but > how > should I write the code to get selected probesets based on an adjusted > p-value > less than or equal to 0.05. > > code: > design <- model.matrix(~factor(rep(1:2, each=3))) > fit <- lmFit(eset, design) > fit2 <- eBayes(fit) > topTable(fit2,coef=2,adjust="BH") > > I have 3 control and 3 treated samples.The eset has the probenames and the > expression matrix. > So my question is how do I get the data of all the probesets which pass > this > criteria? > Thanks, > Viritha > > > > > At 08:51 AM 1/2/2006, Gordon Smyth wrote: > > >Dear Sohail, > > > > > >Well, there are lots of ways to generate such a table. Perhaps the > simplest > is > > > > > > fitsel <- fit2[sel.dif, ] > > > as.data.frame( fitsel ) > > > > > >Best wishes > > >Gordon > > > > > > >Date: Tue, 20 Dec 2005 14:03:43 -0500 > > > >From: "Khan, Sohail" <khan@...> > > > >Subject: [BioC] Time course experiment.... > > > >To: <bioconductor@...> > > > > > > > >Dear List, > > > > > > > >I have performed a time course analysis using limma, as described in > > > >"Bioinformatics and Computational Biology Solutions ........". > > > >How can I get a list of differentially expressed genes? I've tried > > > >the code below: > > > >sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr") <0.05 > > > >This produces a logical vector, right?. I would like a table of > > > >differentially expressed genes with p vales etc. Sorry, if I missed > > > >this in the limma user's guide. Thanks for any suggestions. > > > > > > > > > > > >Sohail Khan > > > >Scientific Programmer > > > >COLD SPRING HARBOR LABORATORY > > > >Genome Research Center > > > >500 Sunnyside Boulevard > > > >Woodbury, NY 11797 > > > >(516)422-4076 > > > > > >_______________________________________________ > > >Bioconductor mailing list > > >Bioconductor@... > > >https://stat.ethz.ch/mailman/listinfo/bioconductor > > > > Naomi S. Altman 814-865-3791 (voice) > > Associate Professor > > Dept. of Statistics 814-863-7114 (fax) > > Penn State University 814-865-1348 (Statistics) > > University Park, PA 16802-2111 > > > > _______________________________________________ > Bioconductor mailing list > Bioconductor@r-project.org<mailto:bioconductor@r-project.org> > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor > > [[alternative HTML version deleted]]
ADD REPLY
0
Entering edit mode
Hi Viritha, the first thing popping to my mind is that you could use the lfc condition in topTable (from limma's Reference Manual -> lfc specifies cutoff value for log2-fold-change. Only genes with larger fold changes are listed.) Thus your code would look like design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) # If you want the first 10000 results... but usually I put "number=THE_NUMBER_OF_GENES_ON_MY_CHIP" top.list<-topTable(fit2,number=10000,coef=2,lfc=1,adjust="BH") your.results <- top.list[top.list$adj.P.Val < 0.05, ] and you would have, at once, only genes whose logFC is higher than 1 or lower than -1 and with and adjusted p value of less than .05 In your code the mistake is this (look the line with the ash code for comment): > your.results1<-your.results[your.results$logFC < -1,] #you are selecting only genes with a logFC lower than -1 and saving them as your.results1 > your.results2<-your.results1[your.results1$logFC > 1,] #now from the previously saved list of genes with logFC < -1 you are asking to sort out those that have a logFC > 1 ...since the previous step was executed correctly, this results in a 0 list (no gene with logFC<-1 can also have logFC>1 of course) > write.table(your.results2,"Panc1_final_list.txt",sep='\t') #well... you are saving an empty object... I hope this helps... You may find other functions of some interest for you, as decideTest, or options for topTable if you read the reference manual -> http://www .bioconductor.org/packages/release/bioc/manuals/limma/man/limma.pdf All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX Maastricht E-mail: m.manca at maastrichtuniversity.nl Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confidential and solely for the use of the intended recipient. It may contain material protected by privacy or attorney-client privilege. If you are not the intended recipient or the person responsible for delivering to the intended recipient, be advised that you have received this email in error and that any use is STRICTLY PROHIBITED. If you have received this email in error please notify us by telephone on +31626441205 Dr Marco MANCA ********************************************************************** *********************************************** ________________________________________ Da: viritha kaza [viritha.k at gmail.com] Inviato: mercoled? 9 febbraio 2011 16.24 A: Manca Marco (PATH) Cc: james.reid at ifom-ieo-campus.it; bioconductor at stat.math.ethz.ch Oggetto: Re: [BioC] Time course experiment.... Hi Marco, Thanks... That happens. I wanted to include one more condition on the top.table list of logFC to consider only those which are greater than or equal to1 and less than or equal to -1. I just tried to expand the code > your.results1<-your.results[your.results$logFC < -1,] > your.results2<-your.results1[your.results1$logFC > 1,] > write.table(your.results2,"Panc1_final_list.txt",sep='\t') but just got a blank txt file with headings.I also tried using & and | operators, but got errors. I hope you could help on this. Thanks, Viritha On Tue, Feb 8, 2011 at 3:32 PM, Manca Marco (PATH) <m.manca at="" maastrichtuniversity.nl<mailto:m.manca="" at="" maastrichtuniversity.nl="">> wrote: Sorry, I had overlooked that part of the question... I thought you just wanted help in setting the threshold so I just added that line to your code. As James has already mentioned to have the full list of genes in your results you have to set something like: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) # If you want the first 10000 results... but usually I put "number=THE_NUMBER_OF_GENES_ON_MY_CHIP" top.list<-topTable(fit2,number=10000,coef=2,adjust="BH") your.results <- top.list[top.list$adj.P.Val < 0.05, ] ...I can only say to my defense "you only see what you look for" =P I should have given a more complete answer since the beginning. All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX Maastricht E-mail: m.manca at maastrichtuniversity.nl<mailto:m.manca at="" maastrichtuniversity.nl=""> Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confidential and solely for the use of the intended recipient. It may contain material protected by privacy or attorney-client privilege. If you are not the intended recipient or the person responsible for delivering to the intended recipient, be advised that you have received this email in error and that any use is STRICTLY PROHIBITED. If you have received this email in error please notify us by telephone on +31626441205 Dr Marco MANCA ********************************************************************** *********************************************** ________________________________________ Da: viritha kaza [viritha.k at gmail.com<mailto:viritha.k at="" gmail.com="">] Inviato: marted? 8 febbraio 2011 19.01 A: Manca Marco (PATH) Cc: bioconductor at stat.math.ethz.ch<mailto:bioconductor at="" stat.math.ethz.ch=""> Oggetto: Re: [BioC] Time course experiment.... Hi Marco, Thanks for your reply. But topTable by default only gives top 10.So how do I know the number in topTable fuction which satisfies the criteria of 0.05 that I need to give. waiting for your reply... Thanks, Viritha On Tue, Feb 8, 2011 at 12:35 PM, Manca Marco (PATH) <m.manca at="" maastrichtuniversity.nl<mailto:m.manca="" at="" maastrichtuniversity.nl=""><mailto:m.manca at="" maastrichtuniversity.nl<mailto:m.manca="" at="" maastrichtuniversity.nl="">>> wrote: Hi Viritha, just do code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) top.list<-topTable(fit2,coef=2,adjust="BH") your.results <- top.list[top.list$adj.P.Val < 0.05, ] ...it should do the trick... All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Netherlands) Visiting address: Experimental Vascular Pathology group, Dept of Pathology - Room5.08, Maastricht University Medical Center, P. Debyelaan 25, 6229 HX Maastricht E-mail: m.manca at maastrichtuniversity.nl<mailto:m.manca at="" maastrichtuniversity.nl=""><mailto:m.manca at="" maastrichtuniversity.nl<mailto:m.manca="" at="" maastrichtuniversity.nl="">> Office telephone: +31(0)433874633 Personal mobile: +31(0)626441205 Twitter: @markomanka ********************************************************************** *********************************************** This email and any files transmitted with it are confidential and solely for the use of the intended recipient. It may contain material protected by privacy or attorney-client privilege. If you are not the intended recipient or the person responsible for delivering to the intended recipient, be advised that you have received this email in error and that any use is STRICTLY PROHIBITED. If you have received this email in error please notify us by telephone on +31626441205 Dr Marco MANCA ********************************************************************** *********************************************** ________________________________________ Da: bioconductor-bounces at r-project.org<mailto:bioconductor-bounces at="" r-project.org=""><mailto:bioconductor-bounces at="" r-project.org<mailto="" :bioconductor-bounces="" at="" r-project.org="">> [bioconductor-bounces at r-project.org<mailto:bioconductor-bounces at="" r-project.org=""><mailto :bioconductor-bounces="" at="" r-project.org<mailto:bioconductor-bounces="" at="" r-project.org="">>] per conto di viritha [viritha.k at gmail.com<mailto:viritha.k at="" gmail.com=""><mailto:viritha.k at="" gmail.com<mailto:viritha.k="" at="" gmail.com="">>] Inviato: marted? 8 febbraio 2011 18.14 A: bioconductor at stat.math.ethz.ch<mailto:bioconductor at="" stat.math.ethz.ch=""><mailto:bioconductor at="" stat.math.ethz.ch<mailto:bioconductor="" at="" stat.math.ethz.ch="">> Oggetto: Re: [BioC] Time course experiment.... Naomi Altman <naomi at="" ...=""> writes: > > I did not get the original posting, but doesn't Sohail just need > "TopTable" for this? > > --Naomi Hi Naomi, I had a similiar issue and used TopTable but I am getting only first 10 but how should I write the code to get selected probesets based on an adjusted p-value less than or equal to 0.05. code: design <- model.matrix(~factor(rep(1:2, each=3))) fit <- lmFit(eset, design) fit2 <- eBayes(fit) topTable(fit2,coef=2,adjust="BH") I have 3 control and 3 treated samples.The eset has the probenames and the expression matrix. So my question is how do I get the data of all the probesets which pass this criteria? Thanks, Viritha > > At 08:51 AM 1/2/2006, Gordon Smyth wrote: > >Dear Sohail, > > > >Well, there are lots of ways to generate such a table. Perhaps the simplest is > > > > fitsel <- fit2[sel.dif, ] > > as.data.frame( fitsel ) > > > >Best wishes > >Gordon > > > > >Date: Tue, 20 Dec 2005 14:03:43 -0500 > > >From: "Khan, Sohail" <khan at="" ...=""> > > >Subject: [BioC] Time course experiment.... > > >To: <bioconductor at="" ...=""> > > > > > >Dear List, > > > > > >I have performed a time course analysis using limma, as described in > > >"Bioinformatics and Computational Biology Solutions ........". > > >How can I get a list of differentially expressed genes? I've tried > > >the code below: > > >sel.dif<-p.adjust(fit2$F.p.vlaue,method="fdr") <0.05 > > >This produces a logical vector, right?. I would like a table of > > >differentially expressed genes with p vales etc. Sorry, if I missed > > >this in the limma user's guide. Thanks for any suggestions. > > > > > > > > >Sohail Khan > > >Scientific Programmer > > >COLD SPRING HARBOR LABORATORY > > >Genome Research Center > > >500 Sunnyside Boulevard > > >Woodbury, NY 11797 > > >(516)422-4076 > > > >_______________________________________________ > >Bioconductor mailing list > >Bioconductor at ... > >https://stat.ethz.ch/mailman/listinfo/bioconductor > > Naomi S. Altman 814-865-3791 (voice) > Associate Professor > Dept. of Statistics 814-863-7114 (fax) > Penn State University 814-865-1348 (Statistics) > University Park, PA 16802-2111 > _______________________________________________ Bioconductor mailing list Bioconductor at r-project.org<mailto:bioconductor at="" r-project.org=""><mailto:bioconductor at="" r-project.org<mailto:bioconductor="" at="" r-project.org="">> https://stat.ethz.ch/mailman/listinfo/bioconductor Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
ADD REPLY
0
Entering edit mode
@gordon-smyth
Last seen 2 hours ago
WEHI, Melbourne, Australia
Dear Viritha, To get all genes that pass the FDR<0.05 criterion, you can use topTable(fit2, coef=2, n=Inf, p=0.05) To get all genes that have FDR<0.05 and have log2-fold-change greater than 1, topTable(fit2, coef=2, n=Inf, p=0.05, lfc=1) See ?topTable. Note that 'n' is the upper limit on the number of genes to list, not necessarily the actual number to list. Best wishes Gordon > Date: Tue, 8 Feb 2011 17:14:42 +0000 > From: viritha <viritha.k at="" gmail.com=""> > To: <bioconductor at="" stat.math.ethz.ch=""> > Subject: Re: [BioC] Time course experiment.... > Message-ID: <loom.20110208t174927-793 at="" post.gmane.org=""> > Content-Type: text/plain; charset="us-ascii" > > > Naomi Altman <naomi at="" ...=""> writes: >> >> I did not get the original posting, but doesn't Sohail just need >> "TopTable" for this? >> >> --Naomi > > Hi Naomi, > I had a similiar issue and used TopTable but I am getting only first 10 but how > should I write the code to get selected probesets based on an adjusted p-value > less than or equal to 0.05. > > code: > design <- model.matrix(~factor(rep(1:2, each=3))) > fit <- lmFit(eset, design) > fit2 <- eBayes(fit) > topTable(fit2,coef=2,adjust="BH") > > I have 3 control and 3 treated samples.The eset has the probenames and the > expression matrix. > So my question is how do I get the data of all the probesets which pass this > criteria? > Thanks, > Viritha ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}
ADD COMMENT
0
Entering edit mode
Hi group, Thanks for helping me on this. Viritha On Wed, Feb 9, 2011 at 5:14 PM, Gordon K Smyth <smyth@wehi.edu.au> wrote: > Dear Viritha, > > To get all genes that pass the FDR<0.05 criterion, you can use > > topTable(fit2, coef=2, n=Inf, p=0.05) > > To get all genes that have FDR<0.05 and have log2-fold-change greater than > 1, > > topTable(fit2, coef=2, n=Inf, p=0.05, lfc=1) > > See ?topTable. Note that 'n' is the upper limit on the number of genes to > list, not necessarily the actual number to list. > > Best wishes > Gordon > > > Date: Tue, 8 Feb 2011 17:14:42 +0000 >> From: viritha <viritha.k@gmail.com> >> To: <bioconductor@stat.math.ethz.ch> >> Subject: Re: [BioC] Time course experiment.... >> Message-ID: <loom.20110208t174927-793@post.gmane.org> >> Content-Type: text/plain; charset="us-ascii" >> >> >> Naomi Altman <naomi@...> writes: >> >>> >>> I did not get the original posting, but doesn't Sohail just need >>> "TopTable" for this? >>> >>> --Naomi >>> >> >> Hi Naomi, >> I had a similiar issue and used TopTable but I am getting only first 10 >> but how >> should I write the code to get selected probesets based on an adjusted >> p-value >> less than or equal to 0.05. >> >> code: >> design <- model.matrix(~factor(rep(1:2, each=3))) >> fit <- lmFit(eset, design) >> fit2 <- eBayes(fit) >> topTable(fit2,coef=2,adjust="BH") >> >> I have 3 control and 3 treated samples.The eset has the probenames and the >> expression matrix. >> So my question is how do I get the data of all the probesets which pass >> this >> criteria? >> Thanks, >> Viritha >> > > ______________________________________________________________________ > The information in this email is confidential and inte...{{dropped:10}}
ADD REPLY

Login before adding your answer.

Traffic: 534 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