# \donttest{
# Load and process data first
nugis_data <- load_builtin_nugis("county", 2016)
processed_data <- agri_process_nugis(nugis_data, "county")
cropland_threshold <- 500 * 2.47105 # 500 ha in acres
# Standard analysis with default 0.5 efficiency
results_default <- agri_classify_nitrogen_custom(
processed_data, cropland_threshold = cropland_threshold, n_efficiency = 0.5
)
# Sensitivity analysis with higher efficiency (e.g., improved management)
results_high <- agri_classify_nitrogen_custom(
processed_data, cropland_threshold = cropland_threshold, n_efficiency = 0.7
)
# Conservative analysis with lower efficiency
results_low <- agri_classify_nitrogen_custom(
processed_data, cropland_threshold = cropland_threshold, n_efficiency = 0.3
)
# Compare classification changes across efficiency scenarios
table(results_default$N_class)
table(results_high$N_class)
table(results_low$N_class)
# }
Run the code above in your browser using DataLab