Hello ctmm community,
Wondering how I can add a map to my AKDE ID 4 plot so it looks like the one shown here?
This is what I have for my code so far:
==========================================
ds.tel.4 <- as.telemetry(df4)
GUESS <- ctmm.guess(ds.tel.4, interactive = FALSE)
FIT <- ctmm.fit(ds.tel.4, GUESS, method = 'pHREML')
# Develop and plot UD
UD <- occurrence(ds.tel.4,FIT)
class(UD)
plot(UD,col.level=NA)
title (main = "Utilization Distribution ID 4")
# Develop and plot AKDE
UD1 <- akde(ds.tel.4,FIT)
plot(ds.tel.4, UD = UD1, error = FALSE, pch = 20)
title (main = "AKDE ID 4")
dev.print(pdf, "UD_AKDE_4.pdf") # save plot to working directory as .pdf
==========================================