From 3f85e8fa94fe81c09953d3f8a8d014dd0f93107b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Timoth=C3=A9e=20Poisot?= <timothee.poisot@umontreal.ca>
Date: Tue, 15 Nov 2022 08:57:12 -0500
Subject: [PATCH] Draft of vignette on stats

---
 docs/src/vignettes/07_statistics_on_layers.jl | 78 +++++++++++++++++++
 1 file changed, 78 insertions(+)
 create mode 100644 docs/src/vignettes/07_statistics_on_layers.jl

diff --git a/docs/src/vignettes/07_statistics_on_layers.jl b/docs/src/vignettes/07_statistics_on_layers.jl
new file mode 100644
index 000000000..9ccc38074
--- /dev/null
+++ b/docs/src/vignettes/07_statistics_on_layers.jl
@@ -0,0 +1,78 @@
+# # Statistics on layers
+
+using SpeciesDistributionToolkit
+using CairoMakie
+using Statistics
+
+# ...
+
+spatial_extent = (left = 3.032, bottom = 55.178, right = 19.687, top = 64.717)
+
+# 
+
+rangifer = taxon("Rangifer tarandus tarandus"; strict = false)
+query = [
+    "occurrenceStatus" => "PRESENT",
+    "hasCoordinate" => true,
+    "country" => "NO",
+    "country" => "SE",
+    "country" => "FI",
+    "country" => "DE",
+]
+presences = occurrences(rangifer, query...)
+
+# layers
+
+dataprovider = RasterData(CHELSA1, BioClim)
+temperature = 0.1SimpleSDMPredictor(dataprovider; layer = "BIO1", spatial_extent...)
+precipitation = 1.0SimpleSDMPredictor(dataprovider; layer = "BIO12", spatial_extent...)
+
+#
+
+mean(temperature)
+
+# 
+
+std(temperature)
+
+# 
+
+z_temperature = (temperature - mean(temperature)) / std(temperature)
+
+# 
+
+fig, ax, hm = heatmap(
+    sprinkle(z_temperature)...;
+    colormap = :roma,
+    colorrange = (-2, 2),
+    figure = (; resolution = (800, 400)),
+    axis = (; aspect = DataAspect()),
+)
+Colorbar(fig[:, end + 1], hm)
+current_figure()
+
+# 
+
+fig, ax, hm = heatmap(
+    sprinkle(rescale(precipitation, (0.0, 1.0)))...;
+    colormap = :bamako,
+    figure = (; resolution = (800, 400)),
+    axis = (; aspect = DataAspect()),
+)
+Colorbar(fig[:, end + 1], hm)
+current_figure()
+
+# 
+
+fig, ax, hm = heatmap(
+    sprinkle(rescale(precipitation, collect(0.0:0.05:1.0)))...;
+    colormap = :bamako,
+    figure = (; resolution = (800, 400)),
+    axis = (; aspect = DataAspect()),
+)
+Colorbar(fig[:, end + 1], hm)
+current_figure()
+
+#
+
+quantile(temperature, [0.05, 0.95])