1. For the return behavior, identify the three stable phases A, B, and A, and calculate AKDE for each. Since A is calculated twice and overlaps to a large extent, I cannot simply consider A+B+A as the annual home range. Therefore, I plan to import ABA into QGIS, merge them together, eliminate overlapping areas, and retain the outside boundary to serve as the annual home range.
2. For the commute behavior, I will try to identify relatively stable A and B. If distinguishable, I will calculate the area as in step 1. If it is difficult to distinguish (e.g., frequent and irregular commuting), I will not segment the data and will calculate it directly, as long as the variograms tend to asymptote.
3. Besides calculating the annual stable home range, I am also interested in understanding the frequency and duration of exploration behaviors. Is it feasible to overlay the home range shapefiles with the tracking positions in GIS to assess occurrences and durations outside the home range?
4. I am not sure if the above methods are optimal, and considering my study involves 41 individuals, each with 2-4 years of data, manually and visually segmenting or fitting models for each individual's annual spatial behavior would be very time-consuming. I hope there is a faster and more reliable method, or if there are any reference R codes available.