I am having an odd problem. Using secr ver. 2.9.5, I fit the following model and got the following results:
BC0 <- secr.fit (CHbc, mask = bcmask, model = list (D=~1, g0 ~ 1, sigma ~ 1),
hcov="Sex", verify = FALSE)
Fitted (real) parameters evaluated at base levels of covariates
session = 1, h2 = F
link estimate SE.estimate lcl ucl
D log 1.246207e-03 1.503725e-04 9.845815e-04 1.577353e-03
g0 logit 2.339183e-01 2.715903e-02 1.849200e-01 2.912602e-01
sigma log 1.353240e+03 8.680046e+01 1.193528e+03 1.534325e+03
pmix logit 3.671874e-01 4.260695e-02 2.882925e-01 4.539032e-01
session = 1, h2 = M
link estimate SE.estimate lcl ucl
D log 1.246207e-03 1.503725e-04 9.845815e-04 1.577353e-03
g0 logit 2.339183e-01 2.715903e-02 1.849200e-01 2.912602e-01
sigma log 1.353240e+03 8.680046e+01 1.193528e+03 1.534325e+03
pmix logit 6.328126e-01 4.260695e-02 5.460968e-01 7.117075e-01
Fitting the same model after upgrading to ver 2.10.2, I get different results:
Fitted (real) parameters evaluated at base levels of covariates
session = 1, h2 = F
link estimate SE.estimate lcl ucl
D log 5.620850e-04 6.775755e-05 4.441833e-04 7.112819e-04
g0 logit 1.398790e-01 1.937748e-02 1.060286e-01 1.823323e-01
sigma log 2.672849e+03 1.664143e+02 2.366078e+03 3.019395e+03
pmix logit 3.671893e-01 4.260742e-02 2.882935e-01 4.539060e-01
session = 1, h2 = M
link estimate SE.estimate lcl ucl
D log 5.620850e-04 6.775755e-05 4.441833e-04 7.112819e-04
g0 logit 1.398790e-01 1.937748e-02 1.060286e-01 1.823323e-01
sigma log 2.672849e+03 1.664143e+02 2.366078e+03 3.019395e+03
pmix logit 6.328107e-01 4.260742e-02 5.460940e-01 7.117065e-01
I have not tried to install the previous version to see if I could recreate the old results but I was wondering if anyone had experienced similar problems. The problem does not occur with all my models. I have only identified this issue with 1 model. It is possible that there could be differences in the CH or mask when I ran the models with the old and new versions, but I have double checked and they all seem the same. The only difference I can see is that the 2 versions used different starting values in the MLE, but I would think that would be expected. Thanks for any advice.
Joe