Hello,
I have a multi-year, multi-season northern quoll trapping data set. Each year includes three trapping sessions and around 22 occasions for each session (however this does vary slightly).
To begin with, I am running simple single session SCR models for one year. The model runs fine for one session but for the other two I am getting NAs. The output is:
> fmaskpre<-make.mask(traps(fpre), buffer= 4 * initialsigmaqfpre, spacing=10, type= "trapbuffer", poly= quollarea)
> fit.HNpref<-secr.fit(fpre, detectfn = "HN", mask=fmaskpre)
Checking data
Preparing detection design matrices
Preparing density design matrix
Finding initial parameter values...
Initial values D = 0.29225, g0 = 0.13642, sigma = 110.3693
Maximizing likelihood...
Eval Loglik D g0 sigma
1 NA -1.2301 -1.8454 4.7038
2 NA -1.2301 -1.8454 4.7038
3 NA -1.2301 -1.8454 4.7038
4 NA -1.2301 -1.8454 4.7038
5 NA -1.2301 -1.8454 4.7038
6 NA -1.2300 -1.8454 4.7038
7 NA -1.2301 -1.8453 4.7038
8 NA -1.2301 -1.8454 4.7043
9 NA -1.2299 -1.8454 4.7038
10 NA -1.2300 -1.8453 4.7038
11 NA -1.2300 -1.8454 4.7043
12 NA -1.2301 -1.8452 4.7038
13 NA -1.2301 -1.8453 4.7043
14 NA -1.2301 -1.8454 4.7048
Completed in 395.47 seconds at 14:15:13 28 Nov 2017
Detector type single
Detector number 200
Average spacing 59.77485 m
x-range 652481.6 654025.7 m
y-range 8468958 8470848 m
N animals : 40
N detections : 187
N occasions : 23
Mask area : 370.37 ha
As suggested in another post I tried included starting values for density, g0 and sigma but I got the same result. I have a mask and have tried variations to buffer width, starting values and detection functions but no success. All my data is formatted the same way and I have had success for all three sessions using the following year. I have double checked my data many many times so I am hoping it is not a silly human error. I am completely stumped and any help or guidance would be very much appreciated.
Thank-you!
Jaime