www.phidot.org

Hello again

I've run into an issue with running some of the secr models. Some of my data is a bit sparse so secr was struggling to run the models;

Error in secr.fit(HYBRIDCH, model = list(D ~ 1, g0 ~ 1), mask = HYBRIDhabitatmask, :
too few values session 1 to determine start;
set manually or select different session with details autoini

I tried the autoini:

autoini(HYBRIDCH, HYBRIDhabitatmask)
Error in autoini(HYBRIDCH, HYBRIDhabitatmask) :
too few values for autoini

So then I tried setting the start values e.g.:

HYBRIDsecrBASE<- secr.fit(HYBRIDCH, model=list(D=0.01,g0=0.01, sigma=1000), mask=HYBRIDhabitatmask, method = "Nelder-Mead",CL=TRUE, detectfn = 'HN',trace = FALSE)

This worked fine, I ran the other models, they seemed to run fine as well, but when I looked at the AIC output all the models came out exactly the same:

AIC(fits)
model detectfn npar logLik AIC AICc dAICc AICcwt
BASE g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
TRP g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
TRN g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
SRSP g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
STRP g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
STTR g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667

So I then re-ran the models using a previous model that had worked e.g:

HYBRIDsecrTRP <- secr.fit(HYBRIDCH, start=ALLsecrBASE, mask=HYBRIDhabitatmask, method =
+ "Nelder-Mead",CL=TRUE, detectfn = 'HN', trace=FALSE)

AIC(fits)
model detectfn npar logLik AIC AICc dAICc AICcwt
BASE g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
TRP g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
TRN g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
SRSP g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
STRP g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667
STTR g0~1 sigma~1 halfnormal 2 -190.1319 384.264 396.264 0 0.1667

Which gave the same problem....I'm obviously doing something wrong but I can't see what.

Any help would be much appreciated!

Thanks

Kerry

Hi Kerry
That does look odd, but we can't see how you generated the models in the 'fits' object. Are you sure the models were different and fitted successfully? I presume you meant

Code: Select all

HYBRIDsecrBASE<- secr.fit(HYBRIDCH, start=list(D=0.01,g0=0.01, sigma=1000), mask=HYBRIDhabitatmask, method = "Nelder-Mead",CL=TRUE, detectfn = 'HN',trace = FALSE)


(and if CL = TRUE then D=0.01 is not used).
Murray