From addd9738806806f850b360ea069f29908c326577 Mon Sep 17 00:00:00 2001
From: Felix Cremer <felix.cremer@bgc-jena.mpg.de>
Date: Tue, 11 Feb 2025 10:14:29 +0100
Subject: [PATCH 1/2] Check for nan in tau_rr

---
 src/rqatrend.jl | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/src/rqatrend.jl b/src/rqatrend.jl
index ad1b68e..8978a17 100644
--- a/src/rqatrend.jl
+++ b/src/rqatrend.jl
@@ -65,8 +65,9 @@ function rqatrend_impl(data; thresh=2, border=10, theiler=1, metric=CheckedEucli
     y_mean = 0.0
     xy_mean = 0.0
     for x in xs
-        n += 1.0
         y = tau_rr(data, x; thresh, metric)
+        isnan(y) && continue
+        n += 1.0
         y_mean = smooth(y_mean, y, inv(n))
         xy_mean = smooth(xy_mean, x*y, inv(n))
     end
@@ -120,6 +121,7 @@ function tau_rr(y, d; thresh=2, metric=CheckedEuclidean())
             nominator += evaluate(metric, y[i], y[i+d]) <= thresh
             denominator += 1
         end
+        @show denominator
         return nominator/denominator
     end
 end

From 9eae887d772bc2157211c65bf9a4ec037b4cead0 Mon Sep 17 00:00:00 2001
From: Felix Cremer <felix.cremer@bgc-jena.mpg.de>
Date: Tue, 11 Feb 2025 10:14:52 +0100
Subject: [PATCH 2/2] Increase Rasters version

---
 Project.toml | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/Project.toml b/Project.toml
index ecf19ab..5ef319d 100644
--- a/Project.toml
+++ b/Project.toml
@@ -60,7 +60,7 @@ LoggingExtras = "1"
 Missings = "1"
 NetCDF = "0.12"
 Proj = "1"
-Rasters = "0.12,0.13"
+Rasters = "0.12,0.13, 0.14"
 RecurrenceAnalysis = "2"
 Statistics = "1.10"
 TestItems = "1.0.0"