Displaying 2 results from an estimated 2 matches for "0.20211".
Did you mean:
0.0211
2006 Mar 04
1
replicated time series - lme?
Dear R-helpers,
I have a time series analysis problem in R:
I want to analyse the output of my simulation model which is proportional
cover of shrubs in a savanna plot for each of 500 successive years. I have
run the model (which includes stochasticity, especially in the initial
conditions) 17 times generating 17 time series of shrub cover.
I am interested in a possible periodicity of shrub
2006 May 20
0
No subject
#short example version of input file with 2 runs and 5 time steps (instead
of 17 runs and 500 time steps)
run t cover
1 1 0.234306
1 2 0.188896
1 3 0.198193
1 4 0.213959
1 5 0.184952
2 1 0.189316
2 2 0.185631
2 3 0.20211
2 4 0.216064
2 5 0.216064
#calculate the correlation of lag 1 over 17 replicates
a<-0
for (i in 1:17)
{
c<-ts( cover[run==i] )
d<-acf( c, lag=1, plot=F)
a<-a+d$acf[2]