I am running a CJS model with 4 samples per summer (June-Sept, 1 month apart) over 3 years, so 9 months between last summer survey of year 1 and first summer survey of year 2.
Here is my code
- Code: Select all
timint = c(1,1,1,9,1,1,1,9,1,1,1)
fil.proc=process.data(cjscmbW, model="CJS", time.intervals=timint)
fil.dll=make.design.data(fil.proc)
If I look at the structure in the .dll I am unclear of the interpretation. I understand that the estimates are monthly estimates, but I am not sure which of the values reflects the summer versus winter (monthly) estimates. (sorry, that is a mouthful I know :{ )
I would assume time 1 is the survival from June-July in year 1, time 2 is the survival from June-July in year 1, and then is time 4 the survival from Sept year 1 to June year 2? And 16 is the survival from Sept year 2 to June of year 3? And 27 is the survival from Aug year 3 to Sep year 3?
- > fil.dll$Phi$time
[1] 1 2 3 4 13 14 15 16 25 26 27 2 3 4 13 14 15 16 25 26 27 3 4 13 14 15 16 25 26 27 4
[32] 13 14 15 16 25 26 27 13 14 15 16 25 26 27 14 15 16 25 26 27 15 16 25 26 27 16 25 26 27 25 26
[63] 27 26 27 27
Levels: 1 2 3 4 13 14 15 16 25 26 27
I want to add a seasonal covariate and I am not sure which index to assign this to.
Would the winter monthly estimate be indexes 4 and 16 or 13 and 25?
My time.intervals coding seems intuitively correct, but am I missing something?
Thanks very much!