bmr_obj <- bootstrap_metaregression(oxybenzone, 15, 100)
slope_pod <- slope_pod_analysis(bmr_obj, 0.0001, 10, 0.1)
pod_and_threshold <- pod_envelope_analysis(bmr_obj, slope_pod, 10,
min(oxybenzone$Concentration), max(oxybenzone$Concentration), 0.1)
plot_metaregression_confidence_envelope(bmr_obj, graph_pod = TRUE,
pod = pod_and_threshold$pod, pod_threshold=pod_and_threshold$threshold)
Run the code above in your browser using DataLab