diff --git a/articles/web_only/basic-intro.html b/articles/web_only/basic-intro.html index 2800d1ffd..6ef2e3f0b 100644 --- a/articles/web_only/basic-intro.html +++ b/articles/web_only/basic-intro.html @@ -406,8 +406,8 @@

Model diagnostics#> #> SAMPLING FOR MODEL 'tmb_generic' NOW (CHAIN 1). #> Chain 1: -#> Chain 1: Gradient evaluation took 0.010077 seconds -#> Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 100.77 seconds. +#> Chain 1: Gradient evaluation took 0.010127 seconds +#> Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 101.27 seconds. #> Chain 1: Adjust your expectations accordingly! #> Chain 1: #> Chain 1: @@ -432,9 +432,9 @@

Model diagnostics#> Chain 1: Iteration: 100 / 101 [ 99%] (Warmup) #> Chain 1: Iteration: 101 / 101 [100%] (Sampling) #> Chain 1: -#> Chain 1: Elapsed Time: 28.079 seconds (Warm-up) -#> Chain 1: 0.2 seconds (Sampling) -#> Chain 1: 28.279 seconds (Total) +#> Chain 1: Elapsed Time: 27.666 seconds (Warm-up) +#> Chain 1: 0.199 seconds (Sampling) +#> Chain 1: 27.865 seconds (Total) #> Chain 1: r <- residuals(m3, "mle-mcmc", mcmc_samples = samps) qqnorm(r) @@ -736,7 +736,6 @@

Time-varying effectsnd$depth_scaled2 <- nd$depth_scaled^2 p <- predict(m4, newdata = nd, se_fit = TRUE, re_form = NA) -#> Warning in sqrt(diag(cov)): NaNs produced ggplot(p, aes(depth_scaled, exp(est), ymin = exp(est - 1.96 * est_se), diff --git a/articles/web_only/basic-intro_files/figure-html/plot-cv-1.png b/articles/web_only/basic-intro_files/figure-html/plot-cv-1.png index f0bd0242a..36fd75b60 100644 Binary files a/articles/web_only/basic-intro_files/figure-html/plot-cv-1.png and b/articles/web_only/basic-intro_files/figure-html/plot-cv-1.png differ diff --git a/articles/web_only/basic-intro_files/figure-html/residuals-1.png b/articles/web_only/basic-intro_files/figure-html/residuals-1.png index b3ea44d49..6bae34c3b 100644 Binary files a/articles/web_only/basic-intro_files/figure-html/residuals-1.png and b/articles/web_only/basic-intro_files/figure-html/residuals-1.png differ diff --git a/articles/web_only/basic-intro_files/figure-html/residuals-2.png b/articles/web_only/basic-intro_files/figure-html/residuals-2.png index ec866f64f..ecc7177ad 100644 Binary files a/articles/web_only/basic-intro_files/figure-html/residuals-2.png and b/articles/web_only/basic-intro_files/figure-html/residuals-2.png differ diff --git a/articles/web_only/basic-intro_files/figure-html/residuals-mcmc-1.png b/articles/web_only/basic-intro_files/figure-html/residuals-mcmc-1.png index 63828f6cd..12687e223 100644 Binary files a/articles/web_only/basic-intro_files/figure-html/residuals-mcmc-1.png and b/articles/web_only/basic-intro_files/figure-html/residuals-mcmc-1.png differ diff --git a/articles/web_only/basic-intro_files/figure-html/tv-depth-eff-1.png b/articles/web_only/basic-intro_files/figure-html/tv-depth-eff-1.png index df7a7946b..873e4fd7b 100644 Binary files a/articles/web_only/basic-intro_files/figure-html/tv-depth-eff-1.png and b/articles/web_only/basic-intro_files/figure-html/tv-depth-eff-1.png differ diff --git a/articles/web_only/bayesian.html b/articles/web_only/bayesian.html index 872cbf56c..67c4ff16d 100644 --- a/articles/web_only/bayesian.html +++ b/articles/web_only/bayesian.html @@ -396,7 +396,7 @@

Passing the model to tmbstan#> omega_s[85] -0.29 0.01 0.21 -0.68 -0.44 -0.30 -0.15 0.11 1111 1.00 #> lp__ 136.06 0.93 9.26 118.11 129.87 136.38 142.46 152.88 100 1.02 #> -#> Samples were drawn using NUTS(diag_e) at Wed Apr 3 17:16:28 2024. +#> Samples were drawn using NUTS(diag_e) at Wed Apr 3 19:10:41 2024. #> For each parameter, n_eff is a crude measure of effective sample size, #> and Rhat is the potential scale reduction factor on split chains (at #> convergence, Rhat=1). diff --git a/articles/web_only/bayesian_files/figure-html/unnamed-chunk-9-1.png b/articles/web_only/bayesian_files/figure-html/unnamed-chunk-9-1.png index 2405e7ab8..06b21ee3d 100644 Binary files a/articles/web_only/bayesian_files/figure-html/unnamed-chunk-9-1.png and b/articles/web_only/bayesian_files/figure-html/unnamed-chunk-9-1.png differ diff --git a/articles/web_only/cross-validation.html b/articles/web_only/cross-validation.html index 62df661a4..4e36466d6 100644 --- a/articles/web_only/cross-validation.html +++ b/articles/web_only/cross-validation.html @@ -213,9 +213,9 @@

Measuring model performance#> Set a parallel `future::plan()` to use parallel processing. m_cv$fold_loglik # fold log-likelihood -#> [1] -1669.255 -1546.654 -1569.066 -1813.842 +#> [1] -1529.945 -1792.257 -1545.111 -1766.426 m_cv$sum_loglik # total log-likelihood -#> [1] -6598.817 +#> [1] -6633.738 +#> [1] -605.4381 @@ -284,9 +284,9 @@

Comparing two or more models# Compare log-likelihoods -- higher is better! m1$sum_loglik -#> [1] -6737.907 +#> [1] -6721.528 m2$sum_loglik -#> [1] -6591.109 +#> [1] -6573.109

Model ensembling diff --git a/articles/web_only/delta-models.html b/articles/web_only/delta-models.html index a5ba9679b..7c27696d9 100644 --- a/articles/web_only/delta-models.html +++ b/articles/web_only/delta-models.html @@ -517,8 +517,8 @@

#> #> Dispersion parameter: 0.94 #> Matérn range: 0.01 -#> Spatial SD: 727.76 -#> Spatiotemporal IID SD: 2066.27 +#> Spatial SD: 727.59 +#> Spatiotemporal IID SD: 2065.79 #> ML criterion at convergence: 5102.136 #> #> See ?tidy.sdmTMB to extract these values as a data frame. diff --git a/articles/web_only/delta-models_files/figure-html/cv-1.png b/articles/web_only/delta-models_files/figure-html/cv-1.png index 8c30e86da..6e46d0b49 100644 Binary files a/articles/web_only/delta-models_files/figure-html/cv-1.png and b/articles/web_only/delta-models_files/figure-html/cv-1.png differ diff --git a/articles/web_only/delta-models_files/figure-html/unnamed-chunk-14-1.png b/articles/web_only/delta-models_files/figure-html/unnamed-chunk-14-1.png index f4f013147..0039810de 100644 Binary files a/articles/web_only/delta-models_files/figure-html/unnamed-chunk-14-1.png and b/articles/web_only/delta-models_files/figure-html/unnamed-chunk-14-1.png differ diff --git a/articles/web_only/delta-models_files/figure-html/unnamed-chunk-19-1.png b/articles/web_only/delta-models_files/figure-html/unnamed-chunk-19-1.png index cb6d5b4cf..ef700095e 100644 Binary files a/articles/web_only/delta-models_files/figure-html/unnamed-chunk-19-1.png and b/articles/web_only/delta-models_files/figure-html/unnamed-chunk-19-1.png differ diff --git a/articles/web_only/delta-models_files/figure-html/unnamed-chunk-20-1.png b/articles/web_only/delta-models_files/figure-html/unnamed-chunk-20-1.png index 29cc9e1c8..0bbb84a2d 100644 Binary files a/articles/web_only/delta-models_files/figure-html/unnamed-chunk-20-1.png and b/articles/web_only/delta-models_files/figure-html/unnamed-chunk-20-1.png differ diff --git a/articles/web_only/index-standardization_files/figure-html/residuals-1.png b/articles/web_only/index-standardization_files/figure-html/residuals-1.png index fd8b14e02..185329b82 100644 Binary files a/articles/web_only/index-standardization_files/figure-html/residuals-1.png and b/articles/web_only/index-standardization_files/figure-html/residuals-1.png differ diff --git a/articles/web_only/index-standardization_files/figure-html/residuals-2.png b/articles/web_only/index-standardization_files/figure-html/residuals-2.png index 4affa267c..6eea69322 100644 Binary files a/articles/web_only/index-standardization_files/figure-html/residuals-2.png and b/articles/web_only/index-standardization_files/figure-html/residuals-2.png differ diff --git a/articles/web_only/index-standardization_files/figure-html/residuals-3.png b/articles/web_only/index-standardization_files/figure-html/residuals-3.png index 3c09795ff..ca91d0c15 100644 Binary files a/articles/web_only/index-standardization_files/figure-html/residuals-3.png and b/articles/web_only/index-standardization_files/figure-html/residuals-3.png differ diff --git a/articles/web_only/residual-checking.html b/articles/web_only/residual-checking.html index b12fd5ca1..440888419 100644 --- a/articles/web_only/residual-checking.html +++ b/articles/web_only/residual-checking.html @@ -229,8 +229,8 @@

MCMC-based randomized-quantile #> #> SAMPLING FOR MODEL 'tmb_generic' NOW (CHAIN 1). #> Chain 1: -#> Chain 1: Gradient evaluation took 0.000969 seconds -#> Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 9.69 seconds. +#> Chain 1: Gradient evaluation took 0.000958 seconds +#> Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 9.58 seconds. #> Chain 1: Adjust your expectations accordingly! #> Chain 1: #> Chain 1: @@ -247,9 +247,9 @@

MCMC-based randomized-quantile #> Chain 1: Iteration: 720 / 800 [ 90%] (Sampling) #> Chain 1: Iteration: 800 / 800 [100%] (Sampling) #> Chain 1: -#> Chain 1: Elapsed Time: 6.685 seconds (Warm-up) -#> Chain 1: 4.746 seconds (Sampling) -#> Chain 1: 11.431 seconds (Total) +#> Chain 1: Elapsed Time: 6.679 seconds (Warm-up) +#> Chain 1: 4.757 seconds (Sampling) +#> Chain 1: 11.436 seconds (Total) #> Chain 1: mcmc_res <- residuals(fit_nb2, type = "mle-mcmc", mcmc_samples = samps) qqnorm(mcmc_res) @@ -500,8 +500,8 @@

The need for one-sample residuals#> #> SAMPLING FOR MODEL 'tmb_generic' NOW (CHAIN 1). #> Chain 1: -#> Chain 1: Gradient evaluation took 0.002174 seconds -#> Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 21.74 seconds. +#> Chain 1: Gradient evaluation took 0.002196 seconds +#> Chain 1: 1000 transitions using 10 leapfrog steps per transition would take 21.96 seconds. #> Chain 1: Adjust your expectations accordingly! #> Chain 1: #> Chain 1: @@ -518,11 +518,11 @@

The need for one-sample residuals#> Chain 1: Iteration: 360 / 400 [ 90%] (Sampling) #> Chain 1: Iteration: 400 / 400 [100%] (Sampling) #> Chain 1: -#> Chain 1: Elapsed Time: 30.974 seconds (Warm-up) -#> Chain 1: 30.443 seconds (Sampling) -#> Chain 1: 61.417 seconds (Total) +#> Chain 1: Elapsed Time: 29.354 seconds (Warm-up) +#> Chain 1: 32.877 seconds (Sampling) +#> Chain 1: 62.231 seconds (Total) #> Chain 1: -#> Warning: The largest R-hat is 1.08, indicating chains have not mixed. +#> Warning: The largest R-hat is 1.09, indicating chains have not mixed. #> Running the chains for more iterations may help. See #> https://mc-stan.org/misc/warnings.html#r-hat #> Warning: Bulk Effective Samples Size (ESS) is too low, indicating posterior means and medians may be unreliable. @@ -540,7 +540,7 @@

The need for one-sample residuals#> Asymptotic one-sample Kolmogorov-Smirnov test #> #> data: r3 -#> D = 0.024522, p-value = 0.5845 +#> D = 0.036042, p-value = 0.1488 #> alternative hypothesis: two-sided

Here that gets us something similar and we would (correctly) fail to reject the hypothesis that \(r_3 \sim @@ -774,7 +774,10 @@

How do those randomized variable and turn it into a uniform-distributed variable.

Now we need to do the same with the Poisson equivalent: ppois(). The ppois(y) function gives us the -cumulative probability density up to value y.

+cumulative probability density up to value y. Say we have a +Poisson variable with a mean (lambda) of 5 and we have an observation +with a value of 3. We can calculate the density up to the value of 3 +as:

 lambda <- 5
 ppois(3, lambda)
diff --git a/articles/web_only/residual-checking_files/figure-html/unnamed-chunk-17-1.png b/articles/web_only/residual-checking_files/figure-html/unnamed-chunk-17-1.png
index 50ad2538a..fb3b74569 100644
Binary files a/articles/web_only/residual-checking_files/figure-html/unnamed-chunk-17-1.png and b/articles/web_only/residual-checking_files/figure-html/unnamed-chunk-17-1.png differ
diff --git a/articles/web_only/spatial-trend-models_files/figure-html/unnamed-chunk-13-1.png b/articles/web_only/spatial-trend-models_files/figure-html/unnamed-chunk-13-1.png
index 8142a55f3..5a996bac8 100644
Binary files a/articles/web_only/spatial-trend-models_files/figure-html/unnamed-chunk-13-1.png and b/articles/web_only/spatial-trend-models_files/figure-html/unnamed-chunk-13-1.png differ
diff --git a/news/index.html b/news/index.html
index 7e9e82510..d3d68cb04 100644
--- a/news/index.html
+++ b/news/index.html
@@ -64,7 +64,8 @@
 
     

sdmTMB (development version)

-