Hello list, I'm trying to estimate average home range across animals per season. I use here a subset of my data, 3 individuals, 7 relocations each, tracked across 3 seasons (actual dataset is ~ 100K, ~30 animals). Not all fish are present in all seasons, and their locations differ between seasons. I need to estimate average home range area across fish for each season and plot it.
I think that Anthony Fischbach's reply from the thread below is what I need to do, I'm just stuck partway there (mainly because I'm very new to spatial analysis and adehabitatHR). http://r-sig-geo.2731867.n2.nabble.com/averaging-home-range-contours-td7499301.html#a7503972 Where I'm at now: ## data subset: df <- structure(list(Season = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Spring", "Summer", "Fall"), class = "factor"), Easting = c(425952.746437786, 425930.124103401, 425904.838430173, 425897.396103185, 425859.77203381, 425769.155466566, 425698.934422453, 425172.162553713, 425201.54784946, 425265.188960297, 425293.907156855, 425334.519543017, 425351.98115297, 425341.657291998, 426499.521534761, 426491.617459401, 426489.21804795, 426458.425998741, 426444.572705444, 426435.13014907, 426428.521796826, 425206.363932968, 425214.055516556, 425214.686198032, 425219.462846981, 425224.610973811, 425225.509004604, 425231.130704745, 425295.438309455, 425295.033258023, 425298.385156621, 425286.605300153, 425292.736058183, 425296.563442322, 425289.068691518, 425262.047903812, 425245.264536305, 425236.068286685, 425246.961026478, 425244.291558159, 425234.858476333, 425220.329721481), Northing = c(5630038.14127036, 5630042.68938361, 5630053.17010795, 5630050.05067675, 5630050.03021924, 5630044.42581987, 5630080.68394433, 5631660.65042482, 5631622.41326505, 5631580.90371775, 5631564.80939618, 5631533.30748204, 5631498.46768363, 5631476.0382073, 5629521.34564057, 5629530.46592294, 5629545.1806574, 5629555.1802463, 5629557.15541335, 5629572.08088964, 5629598.08820948, 5630913.54144105, 5630904.64457731, 5630904.41307655, 5630898.56103167, 5630889.14462961, 5630887.90831366, 5630886.93779994, 5630645.78035996, 5630647.00958702, 5630649.85318956, 5630643.90518234, 5630648.93328633, 5630650.5466496, 5630648.65226268, 5631661.36063509, 5631635.46509059, 5631632.70549606, 5631635.88559859, 5631631.47516021, 5631626.93945528, 5631610.91014489), ID = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c("01", "04", "06"), class = "factor")), .Names = c("Season", "Easting", "Northing", "ID"), row.names = c(NA, 42L), class = "data.frame") library(adehabitatHR) ## to calculate UD for each fish-season combination: df$IDS <- paste(df$ID, df$Season, sep = "") sp_df <- SpatialPoints(subset(df, select = c(Easting, Northing)), proj4string = CRS("+proj=utm +zone=11 +ellps=WGS84")) sp_df <- SpatialPointsDataFrame(coords = sp_df, data = subset(df, select = c(IDS))) kud_full <- kernelUD(sp_df, grid = 300, extent = 0.4, same4all=TRUE) kud_spdf <- estUDm2spixdf(kud_full) mimage(kud_spdf) Now that I have separate UDs, how do I proceed to calculating seasonal averages of the LOCATION (not area) of home ranges? I can calculate area values and average them out, not a problem, but I can't figure out how to estimate the average polygons so that I can then plot them. Completely stuck here... Thank you so much! -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Average-home-range-across-animals-tp7586534.html Sent from the R-sig-geo mailing list archive at Nabble.com. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo