Compare Space Use Results Across Multiple Probability Thresholds

Analyzes how space use areas change across different probability thresholds and compares them to reference track-based space use estimates.

Usage

compare_space_use_thresholds(
  multi_threshold_results,
  reference_tracks,
  reference_space_use,
  methods = NULL,
  overlap_method = "intersection",
  create_plots = TRUE,
  plot_by_fish = TRUE
)

Arguments

multi_threshold_results

Output from calculate_space_use with multiple prob_thresholds

reference_tracks

Track data (e.g., fish_simulation$tracks) for overlap analysis

reference_space_use

Track-based space use results for area comparison

methods

Character vector of methods to compare (default: all available)

overlap_method

Method for calculating track overlap (default: "intersection")

create_plots

Whether to generate visualization plots (default: TRUE)

plot_by_fish

Whether to include individual fish patterns in plots (default: TRUE)

Value

List containing comparison results and plots:

• area_comparison: Data frame with area ratios by threshold

• track_overlap: Data frame with track overlap analysis

• summary_stats: Summary statistics and optimal thresholds

• plots: List of ggplot objects for visualization

Details

This function compares positioning-based space use areas (calculated with different probability thresholds) against reference track-based space use.

Key analyses:

• Area ratio: positioning_area / track_area for each threshold

• Track overlap: percentage of positioning area containing track points

• Threshold sensitivity: how areas change across thresholds

• Optimal threshold: which threshold produces areas closest to track-based

Examples

if (FALSE) {
# Calculate positioning results with multiple thresholds
multi_results <- calculate_space_use(
  track_data = position_points_coords,
  prob_column = "probability",
  prob_thresholds = c(0.1, 0.2, 0.3, 0.4, 0.5),
  methods = c("constrained_convex_hull", "grid_cell_count")
)

# Calculate reference track-based space use
track_results <- calculate_space_use(
  track_data = fish_simulation$tracks,
  methods = c("constrained_convex_hull", "grid_cell_count")
)

# Compare results
comparison <- compare_space_use_thresholds(
  multi_threshold_results = multi_results,
  reference_tracks = fish_simulation$tracks,
  reference_space_use = track_results
)

# View area ratio plots
print(comparison$plots$area_ratio_boxplot)
print(comparison$plots$area_ratio_by_fish)
}