improved CAR1 example

nicholasjclark · Apr 30, 2024 · e82f8b3 · e82f8b3
1 parent 6be041b
commit e82f8b3
Show file tree

Hide file tree

Showing 18 changed files with 97 additions and 87 deletions.
diff --git a/.github/workflows/R-CMD-check-rstan.yaml b/.github/workflows/R-CMD-check-rstan.yaml
@@ -47,6 +47,7 @@ jobs:
           dependencies: NA
           extra-packages: |
               BH
+              RcppEigen
               knitr
               extraDistr
               lubridate

diff --git a/R/RW.R b/R/RW.R
@@ -46,7 +46,7 @@
 #' seasonality <- beta1 * cov1 + beta2 * cov2
 #'
 #'# Take Gaussian observations with error and return
-#' data.frame(y = rnorm(n, mean = x + seasonality, sd = sigma),
+#' data.frame(y = rnorm(n, mean = x + seasonality, sd = sigma_obs),
 #'            season = rep(1:12, 20)[1:n],
 #'            time = cumsum(time_dis))
 #'}
@@ -60,13 +60,14 @@
 #'                              data.frame(series = 'series2'))) %>%
 #'       dplyr::mutate(series = as.factor(series))
 #'
-#'# mvgam with CAR(1) trends and series-level seasonal smooths
-#'mod <- mvgam(formula = y ~ s(season, bs = 'cc',
-#'                             k = 5, by = series),
+#'# mvgam with CAR(1) trends and series-level seasonal smooths; the
+#'# State-Space representation (using trend_formula) will be more efficient
+#'mod <- mvgam(formula = y ~ 1,
+#'             trend_formula = ~ s(season, bs = 'cc',
+#'                                 k = 5, by = trend),
 #'             trend_model = CAR(),
 #'             data = dat,
 #'             family = gaussian(),
-#'             burnin = 300,
 #'             samples = 300,
 #'             chains = 2)
 #'

diff --git a/docs/reference/RW-1.png b/docs/reference/RW-1.png
diff --git a/docs/reference/RW-2.png b/docs/reference/RW-2.png
diff --git a/docs/reference/RW-3.png b/docs/reference/RW-3.png
diff --git a/docs/reference/RW-4.png b/docs/reference/RW-4.png
diff --git a/docs/reference/RW-5.png b/docs/reference/RW-5.png
diff --git a/docs/reference/RW-6.png b/docs/reference/RW-6.png
diff --git a/docs/reference/RW.html b/docs/reference/RW.html
diff --git a/docs/reference/Rplot001.png b/docs/reference/Rplot001.png
diff --git a/docs/reference/Rplot002.png b/docs/reference/Rplot002.png
diff --git a/docs/reference/Rplot003.png b/docs/reference/Rplot003.png
diff --git a/docs/reference/Rplot004.png b/docs/reference/Rplot004.png
diff --git a/docs/reference/Rplot005.png b/docs/reference/Rplot005.png
diff --git a/docs/reference/Rplot006.png b/docs/reference/Rplot006.png
diff --git a/docs/reference/Rplot007.png b/docs/reference/Rplot007.png
diff --git a/man/RW.Rd b/man/RW.Rd
diff --git a/src/mvgam.dll b/src/mvgam.dll