Problems with Smoothing Factor Estimations Function

[Zach-Kellogg]

I am getting "NA" for 9 of 37 of my smoothing factors after running this code. It may be do to getting a warning from my proj4string code. The warning says "Warning message: In proj4string(xy) : CRS object has comment, which is lost in output"

I have made sure to create the functions for reptile home range prior to running this code as well.

tracking.df <- read.csv("All_Turtle_Points_Used_in_Analysis.csv",
stringsAsFactors = FALSE, header = T)
head(tracking.df)
tail(tracking.df)

tracking.df <- tracking.df[!is.na(tracking.df$x) & !is.na(tracking.df$y),]

tracking.sp <- tracking.df

library(sp)
coordinates(tracking.sp) <- c("x", "y")

proj4string(tracking.sp) <- CRS( "+proj=utm +zone=16 +datum=NAD83 +units=m +no_defs" )

library(adehabitatHR)

tracking.mcp <- mcp(tracking.sp[,1], percent = 100)
print(tracking.mcp)

system.time(h.pred <- reptile.kd.h.fn(tracking.sp)) # The sp object should contain at least x,y, & id

h.pred

[jamesepaterson]

The error message occurs when an object doesn't have a CRS, probably because the kernel fitting function is failing with certain smoothing factors on the 7 animals the function returns NA for. If you send a reproducible example (e.g. scaled coordinates to mask site location) then I can try to diagnose why you're getting this error and update the function.

[a-janssen]

Hi James, I am having the same problem for 39 of my 51 lizard home ranges. I didn't find out yet how to solve the CRS-warning and (maybe therefore?) am not able to use your functions properly. Could you find out something yet?
Maybe it has got to do with the new rgdal version? https://r-spatial.org/r/2020/03/17/wkt.html