Hi rmarkers,

I've been trying to estimate daily survival of nightjars using sex and age(adults, juveniles) as covariates. We have captured about 150 individuals over 5 years and we have 311 occasions.

The top ranked model has age dependent survival and sex dependent recapture and looks like this (where M stand for males and H for females):

results.eanomalus$Phi.age.p.sex$results$real #get real estimates

estimate se lcl ucl

Phi gAH c1 a1 t1 0.9992911 0.0001352532 0.9989697 0.9995123

Phi gJH c1 a0 t1 0.9980639 0.0003490498 0.9972436 0.9986403

p gAH c1 a1.0027397260274 t2 0.0058037 0.0012220000 0.0038396 0.0087637

p gAM c1 a1.0027397260274 t2 0.0268471 0.0020530000 0.0231034 0.0311781

fixed note

Phi gAH c1 a1 t1

Phi gJH c1 a0 t1

p gAH c1 a1.0027397260274 t2

p gAM c1 a1.0027397260274 t2

I was trying to run a gof test using release to this dataset and I got the following error:

eanomalus.processed <- process.data(eanomalus, groups = c("age", "sex"),

time.intervals = time.intervals,

age.var = 1, age.unit = age.unit,

initial.ages = c(1,0),nocc = 311)

release.gof(eanomalus.processed)

Length of command exceeded MAXCHR/2 characters.

RELEAS ERROR TERMINATION

Error in (x3 + 4):length(out) : argument of length 0

Looking at a previous post it has been suggested that this may happen with nocc less or equal to 2.

Thus looking at the summary of my data, I can see that actually nocc is 1 and not 311.

data 6 data.frame list

model 1 -none- character

mixtures 1 -none- numeric

freq 4 data.frame list

nocc 1 -none- numeric

nocc.secondary 0 -none- NULL

time.intervals 310 -none- numeric

begin.time 1 -none- numeric

age.unit 1 -none- numeric

initial.ages 4 -none- numeric

group.covariates 2 data.frame list

Is there something I am doing obviously wrong? Thoughts are welcome!