Robust design with Hugging random effect models

[sergi_20]

Hi all,

I am running a robust design analysis in Rmark with the three types of emigration models (random, markovian and no emigration) and I would like to test the RD - Hugging model with random effects, which fit individual random effect for p.

These models are present in Mark but I can't find them in Rmark. Do you know if they will be available in Rmark soon?

Thanks
Sergi

[jlaake]

You are going to have to be more specific about which model you want. Provide me the number from the list in the MARK interface which is provided in Help/Data Types. Is it

146. RDHugginsRE: Robust Design with Huggins' p and c with Random Effects

-jeff

[sergi_20]

Yes, it is the number 146 RDHugginsRE: Robust Design with Huggins' p and c with Random Effects

from the new markmodel list https://github.com/jlaake/RMark/blob/ma ... Models.pdf

Is this available on the latest RMark version? I ouldn't find it on the https://cran.r-project.org/web/packages/RMark/RMark.pdf

Thanks,
Sergi

[jlaake]

No it is not yet but I can add it. Model 145 is available but not yet on the CRAN version. It is available at
[url]https://drive.google.com/open?id=0B77g1ScdUwVeVU9rUE5lVS1kZ2M. I keep pre-cran updated versions on my google drive. The most current is version 2.1.14 at that location.

[sergi_20]

Jeff, I installed the version 2.1.14 (from your google drive), but the model number 146 Robust design Huggins with Random effects is not available.

Error in setup.model(model, nocc, mixtures) :
Invalid type of model = RDHugginsRE Valid types are
CJS
Recovery
Burnham
Known
Closed
Robust
Barker
Multistrata
Brownie
Jolly
Huggins
RDHuggins
Pradel
Pradsen
Pradlambda
Pradrec
RDBarker
POPAN
MSLiveDead
HetClosed
FullHet
Nest
HugHet
HugFullHet
Occupancy
RDOccupPE
RDOccupPG
RDOccupEG
LinkBarker
ORDMS
CRDMS
HCRDMS
HetRDMS
FHetRDMS
HHetRDMS
HFHetRDMS
RDHet
RDFullHet
RDHHet
RDHFHet
RDBarkHug
RDBarkHet
RDBarkFHet
RDBarkHHet
RDBarkHFHet
RDPdGClosed
RDPdGHuggins
RDPdGHet
RDPdGFullHet
RDPdGHugHet
RDPdGHugFullHet
RDPdLClosed
RDPdLHuggins
RDPdLHet
RDPdLFullHet
RDPdLHugHet
RDPdLHugFullHet
RDPdfClosed
RDPdfHuggins
RDPdfHet
RDPdfFullHet
RDPdfHugHet
RDPdfHugFullHet
ORDPdMS
CRDMSOHug
CRDMSOHet
CRDMSOFHet
OccupHet
RDOccupHetPE
RDOccupHetPG
RDOccupHetEG
OccupRNPoisson
OccupRNNegBin
2SpecOccup
LogitNormalMR
PoissonMR
MSOccupancy
OccupRPoisson
OccupRNegBin
ORDMSState
IELogitNormalMR
RDMSMisClass
MultScalOcc
RDMSOccRepro
RDMSOpenMisClass
CJSMixture
PradSenMix
PradLambdaMix
PradelRecMix
LinkBark

[jlaake]

If you reread my message I said model 145 is in that version. It will take some time before I can get to adding 146.

[sergi_20]

I though that was number 146. Thanks for your help Jeff
Sergi

[jlaake]

I just put a new copy of version 2.1.14 on my google drive. It should be dated 7 Dec 2015. It now contains model 146 which is the one that you requested.

I tried the following default model and it worked fine.

Code: Select all

library(RMark)
data(robust)
time.intervals=c(0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0)
mark(robust,model="RDHugginsRE",time.intervals=time.intervals)

--jeff