profile confidence intervals in RMark
Posted: Mon Jan 08, 2018 12:24 pm
I am using the Huggins CRDMS model in RMark. I have a cumbersome dataset - 22 primary periods, 3 states (1 of which is unobserable with p=0), sex as group, 'age model' based on age at marking - so working in the MARK GUI is not really practical.
I can fit a simple model [S~1, p=c~1, Psi~stratum], but get the following error when I set profile.int=TRUE for the same model:
I get this message regardless of whether I specify mlogit0=TRUE or options="SIMANNEAL". Some parameters are fixed (by adding the ddl$Psi$fix field to the design data), so not sure if that is affecting things.
Can anyone provide insight on this?
I can fit a simple model [S~1, p=c~1, Psi~stratum], but get the following error when I set profile.int=TRUE for the same model:
- Code: Select all
Error in wtable[i, cols] = real[pims[i, cols]] :
number of items to replace is not a multiple of replacement length
I get this message regardless of whether I specify mlogit0=TRUE or options="SIMANNEAL". Some parameters are fixed (by adding the ddl$Psi$fix field to the design data), so not sure if that is affecting things.
Can anyone provide insight on this?