release.gof error
Posted: Tue Apr 27, 2021 12:35 pm
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!
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!