statistical test for time course data
5
0
Entering edit mode
chris Jhon ▴ 260
@chris-jhon-5047
Last seen 9.6 years ago
Hi Richard, Thank you for help. In my data ,i have one point which i think it is different from other points and i would like to test statistical significance of the difference of this point. Your suggestion means that there is no direct function in R that i can use,i have to use a package which implement an algorithm. If so, i think edgeR can do the same analysis too,Am i right? Best Reagards, Chris On Thu, Jan 31, 2013 at 11:53 PM, Richard Friedman < friedman@cancercenter.columbia.edu> wrote: > Dear Chris, > > Limma can be used to test between time points > treating each time point as a categorical variable. > The program "EDGE" from the Storey lab, can test whether > there is significant change over a whole time course. > > http://www.ncbi.nlm.nih.gov/pubmed/16357033 > > with hopes that the above helps, > Rich > Richard A. Friedman, PhD > Associate Research Scientist, > Biomedical Informatics Shared Resource > Herbert Irving Comprehensive Cancer Center (HICCC) > Lecturer, > Department of Biomedical Informatics (DBMI) > Educational Coordinator, > Center for Computational Biology and Bioinformatics (C2B2)/ > National Center for Multiscale Analysis of Genomic Networks (MAGNet)/ > Columbia Initiative in Systems Biology > Room 824 > Irving Cancer Research Center > Columbia University > 1130 St. Nicholas Ave > New York, NY 10032 > (212)851-4765 (voice) > friedman@cancercenter.columbia.edu > http://cancercenter.columbia.edu/~friedman/ > > "Complex numbers! Ha! Ha! There is nothing weirder > than imaginary numbers. Architects don't need to know > complex numbers. Whenever I get a negative root for > an area, I throw it out. And don't talk to me about > quaternions. I am not going into computer animation." > -Rose Friedman, age 16 > > > On Jan 30, 2013, at 11:43 PM, chris Jhon wrote: > > Hi All, > > I have data at different time points for time course experiment. > I have a response for each time point and i would like to test whether the > difference between response of two time points is statistically significant > or not. > my data is linear plot where response on y axis and time on x axis. > > what statistical test shall i use? > > > I appreciate any help. > > Best Regards, > Chris > > [[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]]
Cancer limma edgeR Cancer limma edgeR • 2.1k views
ADD COMMENT
0
Entering edit mode
@richard-friedman-513
Last seen 9.6 years ago
Dear Chris, The F-test in Limma will tell you if at least one point is different from the others. It won't tell you which one. The program EDGE (not EdgeR which is a different program entirely) will do the same based on an model which explicitly takes temporal variation into account (this is why I mention it). There is also a t-test that deals with whether a single measurement is a member of the same normal distribution as other measurements. I am not sure off hand, how to do this in R. If you implement this in R and do it for all the rows, it may be the test that you want. If you do this, you will be forgoing the empirical Bayesian increase in accuracy in Limma or the explicit temporal modeling in EDGE, but that is up to you. If you do this, you should current for false discoveries. With hopes that this help, Best wishes, Rich ------------------------------------------------------------ Richard A. Friedman, PhD Associate Research Scientist, Biomedical Informatics Shared Resource Herbert Irving Comprehensive Cancer Center (HICCC) Lecturer, Department of Biomedical Informatics (DBMI) Educational Coordinator, Center for Computational Biology and Bioinformatics (C2B2)/ National Center for Multiscale Analysis of Genomic Networks (MAGNet) Room 824 Irving Cancer Research Center Columbia University 1130 St. Nicholas Ave New York, NY 10032 (212)851-4765 (voice) friedman at cancercenter.columbia.edu http://cancercenter.columbia.edu/~friedman/ In Memoriam, Hymie Simon On Jan 31, 2013, at 10:19 PM, chris Jhon wrote: > Hi Richard, > > Thank you for help. > In my data ,i have one point which i think it is different from other points and i would like to test statistical significance of the difference of this point. > Your suggestion means that there is no direct function in R that i can use,i have to use a package which implement an algorithm. > If so, i think edgeR can do the same analysis too,Am i right? > > Best Reagards, > Chris > > On Thu, Jan 31, 2013 at 11:53 PM, Richard Friedman <friedman at="" cancercenter.columbia.edu=""> wrote: > Dear Chris, > > Limma can be used to test between time points > treating each time point as a categorical variable. > The program "EDGE" from the Storey lab, can test whether > there is significant change over a whole time course. > > http://www.ncbi.nlm.nih.gov/pubmed/16357033 > > with hopes that the above helps, > Rich > Richard A. Friedman, PhD > Associate Research Scientist, > Biomedical Informatics Shared Resource > Herbert Irving Comprehensive Cancer Center (HICCC) > Lecturer, > Department of Biomedical Informatics (DBMI) > Educational Coordinator, > Center for Computational Biology and Bioinformatics (C2B2)/ > National Center for Multiscale Analysis of Genomic Networks (MAGNet)/ > Columbia Initiative in Systems Biology > Room 824 > Irving Cancer Research Center > Columbia University > 1130 St. Nicholas Ave > New York, NY 10032 > (212)851-4765 (voice) > friedman at cancercenter.columbia.edu > http://cancercenter.columbia.edu/~friedman/ > > "Complex numbers! Ha! Ha! There is nothing weirder > than imaginary numbers. Architects don't need to know > complex numbers. Whenever I get a negative root for > an area, I throw it out. And don't talk to me about > quaternions. I am not going into computer animation." > -Rose Friedman, age 16 > > > On Jan 30, 2013, at 11:43 PM, chris Jhon wrote: > >> Hi All, >> >> I have data at different time points for time course experiment. >> I have a response for each time point and i would like to test whether the >> difference between response of two time points is statistically significant >> or not. >> my data is linear plot where response on y axis and time on x axis. >> >> what statistical test shall i use? >> >> >> I appreciate any help. >> >> Best Regards, >> Chris >> >> >> [[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 COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 9 hours ago
WEHI, Melbourne, Australia
Dear Chris, Yes, limma can easily test for a difference at one point, or can test for a significant change over the whole time course like EDGE. I don't understand your experiment well enough to give more specific advice. You would need to tell us your experimental design, in terms of the targets frame, and exactly what you want to test. Best wishes Gordon > Date: Fri, 1 Feb 2013 12:19:26 +0900 > From: chris Jhon <cjhon217 at="" gmail.com=""> > To: Richard Friedman <friedman at="" cancercenter.columbia.edu=""> > Cc: Bioconductor mailing list <bioconductor at="" r-project.org=""> > Subject: Re: [BioC] statistical test for time course data > > Hi Richard, > > Thank you for help. > In my data ,i have one point which i think it is different from other > points and i would like to test statistical significance of the difference > of this point. > Your suggestion means that there is no direct function in R that i can > use,i have to use a package which implement an algorithm. > If so, i think edgeR can do the same analysis too,Am i right? > > Best Reagards, > Chris > > On Thu, Jan 31, 2013 at 11:53 PM, Richard Friedman < > friedman at cancercenter.columbia.edu> wrote: > >> Dear Chris, >> >> Limma can be used to test between time points >> treating each time point as a categorical variable. >> The program "EDGE" from the Storey lab, can test whether >> there is significant change over a whole time course. >> >> http://www.ncbi.nlm.nih.gov/pubmed/16357033 >> >> with hopes that the above helps, >> Rich >> Richard A. Friedman, PhD >> Associate Research Scientist, >> Biomedical Informatics Shared Resource >> Herbert Irving Comprehensive Cancer Center (HICCC) >> Lecturer, >> Department of Biomedical Informatics (DBMI) >> Educational Coordinator, >> Center for Computational Biology and Bioinformatics (C2B2)/ >> National Center for Multiscale Analysis of Genomic Networks (MAGNet)/ >> Columbia Initiative in Systems Biology >> Room 824 >> Irving Cancer Research Center >> Columbia University >> 1130 St. Nicholas Ave >> New York, NY 10032 >> (212)851-4765 (voice) >> friedman at cancercenter.columbia.edu >> http://cancercenter.columbia.edu/~friedman/ >> >> "Complex numbers! Ha! Ha! There is nothing weirder >> than imaginary numbers. Architects don't need to know >> complex numbers. Whenever I get a negative root for >> an area, I throw it out. And don't talk to me about >> quaternions. I am not going into computer animation." >> -Rose Friedman, age 16 >> >> >> On Jan 30, 2013, at 11:43 PM, chris Jhon wrote: >> >> Hi All, >> >> I have data at different time points for time course experiment. >> I have a response for each time point and i would like to test whether the >> difference between response of two time points is statistically significant >> or not. >> my data is linear plot where response on y axis and time on x axis. >> >> what statistical test shall i use? >> >> >> I appreciate any help. >> >> Best Regards, >> Chris >> ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}
ADD COMMENT
0
Entering edit mode
Hi Gordon, I've been curious about this as well. Can limma model the following situation? We have two treatments and multiple time points. We are interested in if the mean profile for each treatment differs over time, treating time continuously. The measurements are over time for each individual so we have to account for correlations within individuals. Ideally, I would like to allow for a random intercept and a random slope (possible quadratic) if needed. EDGE uses splines, so that would be nice as well. I am aware of duplicateCorrelation, Is this the way to proceed? I'll post a reproducible example for guidance if you indicate that the above is possible. Thanks for your time and for your work on limma, which I use frequently. On Fri, Feb 1, 2013 at 6:46 PM, Gordon K Smyth <smyth at="" wehi.edu.au=""> wrote: > Dear Chris, > > Yes, limma can easily test for a difference at one point, or can test for a > significant change over the whole time course like EDGE. > > I don't understand your experiment well enough to give more specific advice. > You would need to tell us your experimental design, in terms of the targets > frame, and exactly what you want to test. > > Best wishes > Gordon > >> Date: Fri, 1 Feb 2013 12:19:26 +0900 >> From: chris Jhon <cjhon217 at="" gmail.com=""> >> To: Richard Friedman <friedman at="" cancercenter.columbia.edu=""> >> Cc: Bioconductor mailing list <bioconductor at="" r-project.org=""> >> Subject: Re: [BioC] statistical test for time course data >> >> Hi Richard, >> >> Thank you for help. >> In my data ,i have one point which i think it is different from other >> points and i would like to test statistical significance of the difference >> of this point. >> Your suggestion means that there is no direct function in R that i can >> use,i have to use a package which implement an algorithm. >> If so, i think edgeR can do the same analysis too,Am i right? >> >> Best Reagards, >> Chris >> >> On Thu, Jan 31, 2013 at 11:53 PM, Richard Friedman < >> friedman at cancercenter.columbia.edu> wrote: >> >>> Dear Chris, >>> >>> Limma can be used to test between time points >>> treating each time point as a categorical variable. >>> The program "EDGE" from the Storey lab, can test whether >>> there is significant change over a whole time course. >>> >>> http://www.ncbi.nlm.nih.gov/pubmed/16357033 >>> >>> with hopes that the above helps, >>> Rich >>> Richard A. Friedman, PhD >>> Associate Research Scientist, >>> Biomedical Informatics Shared Resource >>> Herbert Irving Comprehensive Cancer Center (HICCC) >>> Lecturer, >>> Department of Biomedical Informatics (DBMI) >>> Educational Coordinator, >>> Center for Computational Biology and Bioinformatics (C2B2)/ >>> National Center for Multiscale Analysis of Genomic Networks (MAGNet)/ >>> Columbia Initiative in Systems Biology >>> Room 824 >>> Irving Cancer Research Center >>> Columbia University >>> 1130 St. Nicholas Ave >>> New York, NY 10032 >>> (212)851-4765 (voice) >>> friedman at cancercenter.columbia.edu >>> http://cancercenter.columbia.edu/~friedman/ >>> >>> "Complex numbers! Ha! Ha! There is nothing weirder >>> than imaginary numbers. Architects don't need to know >>> complex numbers. Whenever I get a negative root for >>> an area, I throw it out. And don't talk to me about >>> quaternions. I am not going into computer animation." >>> -Rose Friedman, age 16 >>> >>> >>> On Jan 30, 2013, at 11:43 PM, chris Jhon wrote: >>> >>> Hi All, >>> >>> I have data at different time points for time course experiment. >>> I have a response for each time point and i would like to test whether >>> the >>> difference between response of two time points is statistically >>> significant >>> or not. >>> my data is linear plot where response on y axis and time on x axis. >>> >>> what statistical test shall i use? >>> >>> >>> I appreciate any help. >>> >>> Best Regards, >>> Chris >>> > > ______________________________________________________________________ > The information in this email is confidential and intend...{{dropped:4}} > > _______________________________________________ > 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 Juliet, On Mon, Feb 11, 2013 at 10:39 PM, Juliet Hannah <juliet.hannah at="" gmail.com=""> wrote: > Hi Gordon, > > I've been curious about this as well. Can limma model the following > situation? [snip] I think you'll find a thread that appeared a bit after the one you are replying to: http://thread.gmane.org/gmane.science.biology.informatics.conductor/46 188 If you look in Section 8.6 of the limmaUserGuide (in development): http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc/us ersguide.pdf You'll find a section on time course analysis. 8.6.2, in particular, uses splines. 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
@gordon-smyth
Last seen 9 hours ago
WEHI, Melbourne, Australia
Dear Rich, > Date: Fri, 1 Feb 2013 10:30:16 -0500 > From: Richard Friedman <friedman at="" cancercenter.columbia.edu=""> > To: chris Jhon <cjhon217 at="" gmail.com=""> > Cc: Bioconductor mailing list <bioconductor at="" r-project.org=""> > Subject: Re: [BioC] statistical test for time course data > > Dear Chris, > > The F-test in Limma will tell you if at least one point is different > from the others. It won't tell you which one. The time course example in the limma User's Guide does use an F-test, but limma isn't limited to that approach. > The program EDGE (not EdgeR which is a different program entirely) will > do the same based on an model which explicitly takes temporal variation > into account (this is why I mention it). limma can easily do explicit temporal modelling by fitting time course trends, for example as polynomials or as spline curves. It's just a matter of defining the design matrix. > There is also a t-test that deals with whether a single measurement is a > member of the same normal distribution as other measurements. I am not > sure off hand, how to do this in R. If you implement this in R and do it > for all the rows, it may be the test that you want. If you do this, you > will be forgoing the empirical Bayesian increase in accuracy in Limma or > the explicit temporal modeling in EDGE, but that is up to you. If you do > this, you should current for false discoveries. limma can test for a single observation as an outlier by creating a design matrix column specific for that observation. This produces an empirical Bayes version of the classical t-test for an outlier. Best wishes Gordon > With hopes that this help, > > Best wishes, > Rich > ------------------------------------------------------------ > Richard A. Friedman, PhD > Associate Research Scientist, > Biomedical Informatics Shared Resource > Herbert Irving Comprehensive Cancer Center (HICCC) > Lecturer, > Department of Biomedical Informatics (DBMI) > Educational Coordinator, > Center for Computational Biology and Bioinformatics (C2B2)/ > National Center for Multiscale Analysis of Genomic Networks (MAGNet) > Room 824 > Irving Cancer Research Center > Columbia University > 1130 St. Nicholas Ave > New York, NY 10032 > (212)851-4765 (voice) > friedman at cancercenter.columbia.edu > http://cancercenter.columbia.edu/~friedman/ > > In Memoriam, Hymie Simon > On Jan 31, 2013, at 10:19 PM, chris Jhon wrote: > >> Hi Richard, >> >> Thank you for help. >> In my data ,i have one point which i think it is different from other points and i would like to test statistical significance of the difference of this point. >> Your suggestion means that there is no direct function in R that i can use,i have to use a package which implement an algorithm. >> If so, i think edgeR can do the same analysis too,Am i right? >> >> Best Reagards, >> Chris >> >> On Thu, Jan 31, 2013 at 11:53 PM, Richard Friedman <friedman at="" cancercenter.columbia.edu=""> wrote: >> Dear Chris, >> >> Limma can be used to test between time points >> treating each time point as a categorical variable. >> The program "EDGE" from the Storey lab, can test whether >> there is significant change over a whole time course. >> >> http://www.ncbi.nlm.nih.gov/pubmed/16357033 >> >> with hopes that the above helps, >> Rich >> Richard A. Friedman, PhD >> Associate Research Scientist, >> Biomedical Informatics Shared Resource >> Herbert Irving Comprehensive Cancer Center (HICCC) >> Lecturer, >> Department of Biomedical Informatics (DBMI) >> Educational Coordinator, >> Center for Computational Biology and Bioinformatics (C2B2)/ >> National Center for Multiscale Analysis of Genomic Networks (MAGNet)/ >> Columbia Initiative in Systems Biology >> Room 824 >> Irving Cancer Research Center >> Columbia University >> 1130 St. Nicholas Ave >> New York, NY 10032 >> (212)851-4765 (voice) >> friedman at cancercenter.columbia.edu >> http://cancercenter.columbia.edu/~friedman/ >> >> "Complex numbers! Ha! Ha! There is nothing weirder >> than imaginary numbers. Architects don't need to know >> complex numbers. Whenever I get a negative root for >> an area, I throw it out. And don't talk to me about >> quaternions. I am not going into computer animation." >> -Rose Friedman, age 16 >> >> >> On Jan 30, 2013, at 11:43 PM, chris Jhon wrote: >> >>> Hi All, >>> >>> I have data at different time points for time course experiment. >>> I have a response for each time point and i would like to test whether the >>> difference between response of two time points is statistically significant >>> or not. >>> my data is linear plot where response on y axis and time on x axis. >>> >>> what statistical test shall i use? >>> >>> >>> I appreciate any help. >>> >>> Best Regards, >>> Chris ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}
ADD COMMENT
0
Entering edit mode
Dear Gordon, Thank you very much for the clarification. Now that I think of it, the one-against all is straightforward. However, If there are any worked examples you could point me towards for polynomial and spline modeling of the time series I would greatly appreciate it. I am especially interested in testing the hypothesis that the temporal behavior of 2 treatments are different. Best wishes, Rich On Feb 2, 2013, at 6:50 PM, Gordon K Smyth wrote: > Dear Rich, > >> Date: Fri, 1 Feb 2013 10:30:16 -0500 >> From: Richard Friedman <friedman at="" cancercenter.columbia.edu=""> >> To: chris Jhon <cjhon217 at="" gmail.com=""> >> Cc: Bioconductor mailing list <bioconductor at="" r-project.org=""> >> Subject: Re: [BioC] statistical test for time course data >> >> Dear Chris, >> >> The F-test in Limma will tell you if at least one point is different from the others. It won't tell you which one. > > The time course example in the limma User's Guide does use an F-test, but limma isn't limited to that approach. > >> The program EDGE (not EdgeR which is a different program entirely) will do the same based on an model which explicitly takes temporal variation into account (this is why I mention it). > > limma can easily do explicit temporal modelling by fitting time course trends, for example as polynomials or as spline curves. It's just a matter of defining the design matrix. > >> There is also a t-test that deals with whether a single measurement is a member of the same normal distribution as other measurements. I am not sure off hand, how to do this in R. If you implement this in R and do it for all the rows, it may be the test that you want. If you do this, you will be forgoing the empirical Bayesian increase in accuracy in Limma or the explicit temporal modeling in EDGE, but that is up to you. If you do this, you should current for false discoveries. > > limma can test for a single observation as an outlier by creating a design matrix column specific for that observation. This produces an empirical Bayes version of the classical t-test for an outlier. > > Best wishes > Gordon > >> With hopes that this help, >> >> Best wishes, >> Rich >> ------------------------------------------------------------ >> Richard A. Friedman, PhD >> Associate Research Scientist, >> Biomedical Informatics Shared Resource >> Herbert Irving Comprehensive Cancer Center (HICCC) >> Lecturer, >> Department of Biomedical Informatics (DBMI) >> Educational Coordinator, >> Center for Computational Biology and Bioinformatics (C2B2)/ >> National Center for Multiscale Analysis of Genomic Networks (MAGNet) >> Room 824 >> Irving Cancer Research Center >> Columbia University >> 1130 St. Nicholas Ave >> New York, NY 10032 >> (212)851-4765 (voice) >> friedman at cancercenter.columbia.edu >> http://cancercenter.columbia.edu/~friedman/ >> >> In Memoriam, Hymie Simon >> On Jan 31, 2013, at 10:19 PM, chris Jhon wrote: >> >>> Hi Richard, >>> >>> Thank you for help. >>> In my data ,i have one point which i think it is different from other points and i would like to test statistical significance of the difference of this point. >>> Your suggestion means that there is no direct function in R that i can use,i have to use a package which implement an algorithm. >>> If so, i think edgeR can do the same analysis too,Am i right? >>> >>> Best Reagards, >>> Chris >>> >>> On Thu, Jan 31, 2013 at 11:53 PM, Richard Friedman <friedman at="" cancercenter.columbia.edu=""> wrote: >>> Dear Chris, >>> >>> Limma can be used to test between time points >>> treating each time point as a categorical variable. >>> The program "EDGE" from the Storey lab, can test whether >>> there is significant change over a whole time course. >>> >>> http://www.ncbi.nlm.nih.gov/pubmed/16357033 >>> >>> with hopes that the above helps, >>> Rich >>> Richard A. Friedman, PhD >>> Associate Research Scientist, >>> Biomedical Informatics Shared Resource >>> Herbert Irving Comprehensive Cancer Center (HICCC) >>> Lecturer, >>> Department of Biomedical Informatics (DBMI) >>> Educational Coordinator, >>> Center for Computational Biology and Bioinformatics (C2B2)/ >>> National Center for Multiscale Analysis of Genomic Networks (MAGNet)/ >>> Columbia Initiative in Systems Biology >>> Room 824 >>> Irving Cancer Research Center >>> Columbia University >>> 1130 St. Nicholas Ave >>> New York, NY 10032 >>> (212)851-4765 (voice) >>> friedman at cancercenter.columbia.edu >>> http://cancercenter.columbia.edu/~friedman/ >>> >>> "Complex numbers! Ha! Ha! There is nothing weirder >>> than imaginary numbers. Architects don't need to know >>> complex numbers. Whenever I get a negative root for >>> an area, I throw it out. And don't talk to me about >>> quaternions. I am not going into computer animation." >>> -Rose Friedman, age 16 >>> >>> >>> On Jan 30, 2013, at 11:43 PM, chris Jhon wrote: >>> >>>> Hi All, >>>> >>>> I have data at different time points for time course experiment. >>>> I have a response for each time point and i would like to test whether the >>>> difference between response of two time points is statistically significant >>>> or not. >>>> my data is linear plot where response on y axis and time on x axis. >>>> >>>> what statistical test shall i use? >>>> >>>> >>>> I appreciate any help. >>>> >>>> Best Regards, >>>> Chris > > ______________________________________________________________________ > The information in this email is confidential and inte...{{dropped:6}}
ADD REPLY
0
Entering edit mode
@gordon-smyth
Last seen 9 hours ago
WEHI, Melbourne, Australia
Dear Rich, I have added a time course example using splines to the limma User's Guide, see page 48: http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc/us ersguide.pdf Best wishes Gordon ------------------ original message ------------------ [BioC] statistical test for time course data Richard Friedman friedman at cancercenter.columbia.edu Sun Feb 3 20:18:03 CET 2013 Dear Gordon, Thank you very much for the clarification. Now that I think of it, the one-against all is straightforward. However, If there are any worked examples you could point me towards for polynomial and spline modeling of the time series I would greatly appreciate it. I am especially interested in testing the hypothesis that the temporal behavior of 2 treatments are different. Best wishes, Rich ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}
ADD COMMENT
0
Entering edit mode
On Feb 5, 2013, at 4:02 AM, Gordon K Smyth <smyth at="" wehi.edu.au=""> wrote: > Dear Rich, > > I have added a time course example using splines to the limma User's > Guide, see page 48: > > http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc/ usersguide.pdf > Does this work also for RNA-seq data? d > Best wishes > Gordon > > ------------------ original message ------------------ > [BioC] statistical test for time course data > Richard Friedman friedman at cancercenter.columbia.edu > Sun Feb 3 20:18:03 CET 2013 > > Dear Gordon, > > Thank you very much for the clarification. Now that I think of it, > the one-against all is straightforward. However, If there are any worked > examples you could point me towards for polynomial and spline modeling of > the time series I would greatly appreciate it. I am especially interested > in testing the hypothesis that the temporal behavior of 2 treatments are > different. > > Best wishes, > Rich > > ______________________________________________________________________ > The information in this email is confidential and intend...{{dropped:4}} > > _______________________________________________ > 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 /* Davide Cittaro, PhD Coordinator of Bioinformatics Core Center for Translational Genomics and Bioinformatics Ospedale San Raffaele Via Olgettina 58 20132 Milano Italy Office: +39 02 26439140 Mail: cittaro.davide at hsr.it Skype: daweonline */
ADD REPLY
0
Entering edit mode
Gordon, Thank you. The example is clear. Is there a minimum # of time points for which this works? Or to phrase it somewhat differently, what is the relationship between the n (the degree of the spline) and t, the minimum number of time points necessary to parameterize the spline? Also, how do I know how many parameters is in the model and which coefficients to use? Can you please point me at a reference which treats this kind of model? Thanks and best wishes, Rich On Feb 4, 2013, at 10:02 PM, Gordon K Smyth wrote: > Dear Rich, > > I have added a time course example using splines to the limma User's Guide, see page 48: > > http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc/ usersguide.pdf > > Best wishes > Gordon > > ------------------ original message ------------------ > [BioC] statistical test for time course data > Richard Friedman friedman at cancercenter.columbia.edu > Sun Feb 3 20:18:03 CET 2013 > > Dear Gordon, > > Thank you very much for the clarification. Now that I think of it, the one-against all is straightforward. However, If there are any worked examples you could point me towards for polynomial and spline modeling of the time series I would greatly appreciate it. I am especially interested in testing the hypothesis that the temporal behavior of 2 treatments are different. > > Best wishes, > Rich > > ______________________________________________________________________ > The information in this email is confidential and inte...{{dropped:6}}
ADD REPLY
0
Entering edit mode
Hi All, Thank you Gordon and Richard very much. In my data,for each time point i have the number of expressed genes and i would like to find if the number of expressed genes at t1 is different from number of expressed genes at t2 ,or is different from all other time point using statistical test. the data look like this: time t1 t2 .... tn expressed genes # # ......# I have only one group,Shall i use same design matrix ? shall i use df=5 as in example?? Best Regards, Chris On Tue, Feb 5, 2013 at 12:02 PM, Gordon K Smyth <smyth@wehi.edu.au> wrote: > Dear Rich, > > I have added a time course example using splines to the limma User's > Guide, see page 48: > > http://bioconductor.org/**packages/2.12/bioc/vignettes/** > limma/inst/doc/usersguide.pdf<http: bioconductor.org="" packages="" 2.12="" bioc="" vignettes="" limma="" inst="" doc="" usersguide.pdf=""> > > Best wishes > Gordon > > ------------------ original message ------------------ > [BioC] statistical test for time course data > Richard Friedman friedman at cancercenter.columbia.edu > Sun Feb 3 20:18:03 CET 2013 > > Dear Gordon, > > Thank you very much for the clarification. Now that I think of it, > the one-against all is straightforward. However, If there are any worked > examples you could point me towards for polynomial and spline modeling of > the time series I would greatly appreciate it. I am especially interested > in testing the hypothesis that the temporal behavior of 2 treatments are > different. > > Best wishes, > Rich > > ______________________________**______________________________**____ ______ > The information in this email is confidential and inte...{{dropped:13}}
ADD REPLY
0
Entering edit mode
Dear Chris, For the questions you are asking I recommend not using splines. For the comparison of t1 vs t2 us a design matrix which makes every point a different time point and then do t2 vs t1, For t1 compared to all other points, I would label t1 A, and all other points B. If anyone on the list has a different opinion in the matter I would appreciate hearing from them. With hopes that this helps, Rich On Feb 5, 2013, at 10:09 AM, chris Jhon wrote: > Hi All, > > Thank you Gordon and Richard very much. > > In my data,for each time point i have the number of expressed genes and i would like to find if the number of expressed genes at t1 is different from number of expressed genes at t2 ,or is different from all other time point using statistical test. > > the data look like this: > > time t1 t2 .... tn > expressed genes # # ......# > > I have only one group,Shall i use same design matrix ? shall i use df=5 as in example?? > > > Best Regards, > Chris > > On Tue, Feb 5, 2013 at 12:02 PM, Gordon K Smyth <smyth@wehi.edu.au> wrote: > Dear Rich, > > I have added a time course example using splines to the limma User's Guide, see page 48: > > http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc/ usersguide.pdf > > Best wishes > Gordon > > ------------------ original message ------------------ > [BioC] statistical test for time course data > Richard Friedman friedman at cancercenter.columbia.edu > Sun Feb 3 20:18:03 CET 2013 > > Dear Gordon, > > Thank you very much for the clarification. Now that I think of it, the one-against all is straightforward. However, If there are any worked examples you could point me towards for polynomial and spline modeling of the time series I would greatly appreciate it. I am especially interested in testing the hypothesis that the temporal behavior of 2 treatments are different. > > Best wishes, > Rich > > ______________________________________________________________________ > The information in this email is confidential and inte...{{dropped:14}}
ADD REPLY
0
Entering edit mode
On Feb 5, 2013, at 4:02 AM, Gordon K Smyth <smyth at="" wehi.edu.au=""> wrote: > Dear Rich, > > I have added a time course example using splines to the limma User's > Guide, see page 48: > > http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc/ usersguide.pdf > What if multiple conditions are available? d
ADD REPLY
0
Entering edit mode
@gordon-smyth
Last seen 9 hours ago
WEHI, Melbourne, Australia
> Date: Tue, 5 Feb 2013 07:23:32 +0100 > From: Davide Cittaro <cittaro.davide at="" hsr.it=""> > To: Bioconductor mailing list <bioconductor at="" r-project.org=""> > Subject: Re: [BioC] statistical test for time course data > > > On Feb 5, 2013, at 4:02 AM, Gordon K Smyth <smyth at="" wehi.edu.au=""> wrote: > >> Dear Rich, >> >> I have added a time course example using splines to the limma User's >> Guide, see page 48: >> >> http://bioconductor.org/packages/2.12/bioc/vignettes/limma/inst/doc /usersguide.pdf > > > Does this work also for RNA-seq data? > > d Yes. The same approach works for edgeR glms with hardly any change. Best wishes Gordon ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:4}}
ADD COMMENT

Login before adding your answer.

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