This function allows the point estimates and confidence intervals of an initial estimated movement model to be improved by parametric boostrap, as described in Fleming et al (2019).

ctmm.boot(data,CTMM,method=CTMM$method,AICc=FALSE,iterate=FALSE,robust=FALSE,error=0.01,
          cores=1,trace=TRUE,...)

Arguments

data

Timeseries data represented as a telemetry object.

CTMM

A ctmm movement-model object from the output of ctmm.fit containing the initial parameter estimates.

method

Fitting method to use: "ML", "HREML", "pREML", "pHREML", or "REML". See ctmm.fit for descriptions.

AICc

Run dual set of simulations to approximate AICc values via Kullback–Leibler divergence. Otherwise, only the AIC is updated.

iterate

Iteratively solve for the parameters such that the average estimate (of method) is that of the data, whereas with iterate=FALSE only the first-order correction is calculated from the initial estimate.

robust

Uses robust estimates of the average and covariation for debiasing. Useful when parameters are near boundaries.

error

Relative standard error target for bootstrap ensemble estimates and nonlinear iterations.

cores

Number of simulations to run in parallel. cores=NULL will use all cores, while cores<0 will reserve abs(cores).

trace

Report progress updates. Can be among 0:2 with increasing detail.

...

Further arguments passed to ctmm.fit.

Value

A model fit object with relatively unbiased estimates of location covariance, and autocorrelation timescales (and more accurate CIs than ctmm.fit). If AICc=TRUE, then, in addition to an updated AICc slot, the model fit object will also contain a VAR.AICc slot quantifying the numerical variance in the AICc estimate. This variance can be decreased by decreasing argument error.

References

C. H. Fleming, M. J. Noonan, E. P. Medici, J. M. Calabrese, ``Overcoming the challenge of small effective sample sizes in home-range estimation'', Methods in Ecology and Evolution 10:10, 1679-1689 (2019) doi: 10.1111/2041-210X.13270 .

Author

C. H. Fleming.

See also

Examples

# \donttest{ # Load package and data library(ctmm) data(gazelle) DATA <- gazelle[[3]] GUESS <- ctmm.guess(DATA,interactive=FALSE) FIT <- ctmm.select(DATA,GUESS) # some human-readable information summary(FIT)
#> $name #> [1] "OUf isotropic" #> #> $DOF #> mean area speed #> 3.719081 4.105042 18.886592 #> #> $CI #> low est high #> area (square kilometers) 7031.905002 25246.63770 54889.386848 #> τ (days) 2.614605 4.43780 7.532331 #> speed (kilometers/day) 9.048085 11.67099 14.288631 #>
# in general, you will want to set iterate=TRUE BOOT <- ctmm.boot(DATA,FIT,iterate=FALSE)
#> | | | 0% | |=== | 4% | |==== | 5% | |========== | 15% | |========== | 14% | |=========== | 15% | |=============== | 21% | |================ | 23% | |===================== | 30% | |========================= | 35% | |======================== | 34% | |======================= | 33% | |======================= | 34% | |======================== | 34% | |========================= | 36% | |======================== | 34% | |=========================== | 38% | |=========================== | 39% | |========================== | 38% | |========================== | 36% | |========================= | 35% | |========================= | 36% | |========================= | 35% | |========================= | 36% | |=========================== | 38% | |=========================== | 39% | |============================= | 42% | |============================= | 41% | |============================== | 43% | |================================= | 47% | |=================================== | 50% | |==================================== | 51% | |=================================== | 50% | |==================================== | 51% | |=================================== | 51% | |==================================== | 51% | |======================================= | 56% | |========================================== | 60% | |========================================= | 58% | |======================================================== | 79% | |=========================================================== | 84% | |====================================================== | 77% | |===================================================== | 76% | |==================================================== | 75% | |=================================================== | 73% | |===================================================== | 75% | |==================================================== | 74% | |====================================================== | 77% | |======================================================== | 80% | |========================================================= | 81% | |========================================================== | 82% | |======================================================== | 80% | |============================== | 42% | |==================================== | 51% | |========================================== | 60% | |============================================ | 63% | |=================================================== | 73% | |=============================================== | 68% | |============================================ | 63% | |================================================== | 72% | |===================================================== | 76% | |============================================= | 64% | |============================================ | 63% | |============================================= | 64% | |============================================ | 63% | |============================================== | 66% | |============================================== | 65% | |=============================================== | 67% | |============================================== | 66% | |================================================ | 68% | |================================================== | 71% | |==================================================== | 74% | |=================================================== | 73% | |==================================================== | 75% | |===================================================== | 76% | |========================================================= | 81% | |======================================================== | 81% | |========================================================= | 82% | |============================================================= | 88% | |================================================================ | 91% | |================================================================= | 93% | |================================================================ | 92% | |================================================================= | 92% | |=================================================================== | 95% | |================================================================== | 94%
# compare to the previous estimate summary(BOOT)
#> $name #> [1] "OUf isotropic" #> #> $DOF #> mean area speed #> 3.686173 4.144028 14.043887 #> #> $CI #> low est high #> area (square kilometers) 7370.757396 26254.711144 56913.517959 #> τ (days) 2.544622 4.497792 7.950155 #> speed (kilometers/day) 8.687447 11.742965 14.792744 #>
# }