Okay, done! Thanks very much.
Nicole
mask <- make.mask(traps(mouse17), buffer = 80, type = 'trapbuffer')
fit17a <- secr.fit(mouse17, model = list(D~1, g0~1, sigma~1), buffer = 80,
detectfn = 'EX', mask = mask)
## doesn't converge well...
## default mask spacing is wider than trap spacing - disconcerting and possibly problematic
## try a finer mask
mask5 <- make.mask(traps(mouse17), buffer = 80, type = 'trapbuffer', spacing = 5)
nrow(mask5)
## this will be slow...
fit17b <- secr.fit(mouse17, model = list(D~1, g0~1, sigma~1), buffer=80,
detectfn = 'EX', mask = mask5)
## no SE for sigma-hat, so recompute
fit17c <- secr.fit(mouse17, model = list(D~1, g0~1, sigma~1), buffer=80,
detectfn = 'EX', mask = mask5, start=fit17b, method = 'none')
predict(fit17a)
# link estimate SE.estimate lcl ucl
# D log 2.1939802 0.408944556 1.5273081 3.15165579
# g0 logit 0.0550488 0.007925989 0.0414255 0.07281199
# sigma log 18.1718790 NaN NaN NaN
predict(fit17b)
# link estimate SE.estimate lcl ucl
# D log 2.4744063 0.430321629 1.76417096 3.47057428
# g0 logit 0.0584761 0.008517707 0.04385152 0.07758227
# sigma log 18.3024897 NaN NaN NaN
predict(fit17c)
# link estimate SE.estimate lcl ucl
# D log 2.4744063 0.423730612 1.77320388 3.4528948
# g0 logit 0.0584761 0.001954107 0.05476171 0.0624258
# sigma log 18.3024897 0.186117364 17.94132677 18.6709229
Users browsing this forum: No registered users and 7 guests