Thanks Kinsford. I thought it would be appropriate. As a follow up question: My first thought is to set up the data file with three columns: year, population (A,B), and nest success and then to input the following formula: success.gls=gls(success~year*population). This would allow me to test for the effect of year for each population and then also test for differences between the two populations. My questions are 1) have I specified the model right for those questions and 2) would the acf function calculate the autocorrelation correctly even though my 'year' in the data file is repeated twice (once for each value of nest success/population)? Thanks. Lisa (hope all is well with you)
Kingsford Jones wrote: > > The gls function in nlme fits a general linear model, so yes you can > have categorical predictors (the advantage over the lm function is the > error covariance matrix may have non-zero off-diagonals, such as with > an autocorrelation structure, and non-constant diagonals). > > hth, > Kingsford Jones > > On Fri, Jan 1, 2010 at 2:44 PM, LisaB <lisaba...@hotmail.com> wrote: >> >> Hello - >> >> I need to analyze some time series data in an ANOVA framework, but am >> unsure >> of how to go about it. I have data on nest success (response) over a 22 >> year period for two populations. For each year I have one value of nest >> success per population. I am interested in determining 1) whether there >> are >> differences in nest success over time between these two populations and >> 2) >> what are the trends for each population over time. My thought is to use >> GLS >> and model temporal autocorrelation if the acf function indicates this is >> an >> issue, but since population is a categorical variable I'm unsure if this >> is >> appropriate. Any advice would be much appreciated. Thank you. Lisa >> >> >> >> >> -- >> View this message in context: >> http://n2.nabble.com/Time-series-and-GLS-tp4240700p4240700.html >> Sent from the r-sig-ecology mailing list archive at Nabble.com. >> >> _______________________________________________ >> R-sig-ecology mailing list >> R-sig-ecology@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology >> > > _______________________________________________ > R-sig-ecology mailing list > R-sig-ecology@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-ecology > > -- View this message in context: http://n2.nabble.com/Time-series-and-GLS-tp4240700p4244082.html Sent from the r-sig-ecology mailing list archive at Nabble.com. _______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
