From 569827936dde3e3b65d6842ab7a10e2b9d1c9269 Mon Sep 17 00:00:00 2001
From: Nicholas Clark <uqnclar2@uq.edu.au>
Date: Sun, 24 Mar 2024 22:12:45 +1000
Subject: [PATCH] small bug fix in variance calculation for t

---
 R/compute_edf.R |   5 ++++-
 src/mvgam.dll   | Bin 1084416 -> 1084416 bytes
 2 files changed, 4 insertions(+), 1 deletion(-)

diff --git a/R/compute_edf.R b/R/compute_edf.R
index 7d7c9797..9b532b0f 100644
--- a/R/compute_edf.R
+++ b/R/compute_edf.R
@@ -79,7 +79,10 @@ compute_edf = function(mgcv_model, object, rho_names, sigma_raw_names){
     # Calculate variance using family's mean-variance relationship
     mu_variance <- predict(object,
                            process_error = FALSE,
-                           type = 'variance')[best_draw, 1:length(eta)]
+                           type = 'variance')[best_draw, ]
+    if(length(mu_variance) > 1){
+      mu_variance <- mu_variance[1:length(eta)]
+    }
     if(any(mu_variance == 0)){
       mu_variance[which(mu_variance == 0)] <-
         mu[which(mu_variance == 0)]
diff --git a/src/mvgam.dll b/src/mvgam.dll
index 5cea10f45707cfc5bde5c5c7e5cff3dbf05ffdd4..325558e078036a9f195159a0b4a81e4eff7cb231 100644
GIT binary patch
delta 82
zcmV-Y0ImOkm_mS<LXe0A=@S5Ek&H|Nau$Kwh1vmy+5&~z1BKcIh1vy$+6IN%2Zh=Q
oh1v;)+6sl*3x(PYh1v~;+75-<54GA5Dlr7;698tn2r>}(S<y%%I{*Lx

delta 82
zcmV-Y0ImOkm_mS<LXe0AP)zz|k&H|NG>3uOh1vmy+5&~z1BKcIh1vy$+6IN%2Zh=Q
oh1v;)+6sl*3x(PYh1v~;+75-<54GA5Dlr64O!{TF2r>}(SqmT}pa1{>