From 74d6cac72db9f22cb9bb4315ca418d006af94117 Mon Sep 17 00:00:00 2001
From: Chitu Okoli <Chitu.Okoli@skema.edu>
Date: Mon, 10 Feb 2025 09:58:04 +0100
Subject: [PATCH] fix: Use method='fREML' by default

---
 NEWS.md     | 1 +
 R/autogam.R | 1 +
 2 files changed, 2 insertions(+)

diff --git a/NEWS.md b/NEWS.md
index 7e0ecd7..2cc87a7 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -1,5 +1,6 @@
 # autogam 0.1.0
 
+* Use bam(discrete=TRUE, method='fREML') by default.
 * Update required R version to >= 4.1.0 (mainly because it uses the |> pipe).
 * Create `autogam` S3 object with support for all `mgcv::gam` s3 methods. The `autogam` methods mostly call the `mgcv::gam` methods. For now, only the `print()` method is lightly modified to add performance measures.
 
diff --git a/R/autogam.R b/R/autogam.R
index 3ee1c13..f039065 100644
--- a/R/autogam.R
+++ b/R/autogam.R
@@ -72,6 +72,7 @@ autogam <- function(
   # Default to bam
   gam_fun <- bam
   args$discrete <- TRUE
+  args$method <- 'fREML'
 
   ## Detect interactions ---------------------------