# NOT RUN {
# Extract data and crawl inputs from miExample
obsData <- miExample$obsData
inits <- miExample$inits
err.model <- miExample$err.model
# Fit crawl to obsData
crwOut <- crawlWrap(obsData,ncores=1,theta=c(4,0),fixPar=c(1,1,NA,NA),
initial.state=inits,err.model=err.model)
# Fit four imputations
bPar <- miExample$bPar
HMMfits <- MIfitHMM(crwOut,nSims=4,ncores=1,poolEstimates=FALSE,
nbStates=2,dist=list(step="gamma",angle="vm"),
Par0=bPar$Par,beta0=bPar$beta,delta0=bPar$delta,
formula=~cov1+cos(cov2),
estAngleMean=list(angle=TRUE),
covNames=c("cov1","cov2"))
miHMM <- momentuHMM:::miHMM(list(miSum=MIpool(HMMfits,ncores=1),HMMfits=HMMfits))
print(miHMM)
# }
Run the code above in your browser using DataLab