gmfamm

Running examples in ‘gmfamm-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: trafficfam3
> ### Title: Draft of family for traffic example
> ### Aliases: trafficfam3
> 
> ### ** Examples
> 
> # Construct data
> set.seed(123)
> # Number of subjects
> n <- 10
> 
> # Number of observations
> ni <- 3
> 
> # Covariate vector
> x <- rep(rnorm(n), each = ni)
> t <- rep(c(0, 0.5, 1), times = n)
> 
> # Additive predictor
> eta_1 <- eta_2 <- eta_3 <- eta_4 <- t + 0.5*x
> 
> # Outcomes
> y1 <- rnbinom(n*ni, exp(eta_1), 0.3)
> y2 <- rnbinom(n*ni, exp(eta_2), 0.4)
> y3 <- rgamma(n*ni, shape = 0.3, scale = exp(eta_3) / 0.3)
> y4 <- rgamma(n*ni, shape = 0.4, scale = exp(eta_4) / 0.4)
> 
> # Data format
> dat <- data.frame(
+    id = factor(rep(seq_len(n), each = ni)),
+    y = c(y1, y2, y3, y4),
+    dim = factor(rep(c(1:4), each = n*ni)),
+    t = t,
+    x = x,
+    fpc = 1
+ )
> 
> # Specify formula
> f <- list(
+   gm(y, dim) ~ t + x,
+   sigma1 ~ 1,
+   mu2 ~ t + x,
+   sigma2 ~ 1,
+   mu3 ~ t + x,
+   sigma3 ~ 1,
+   mu4 ~ t + x,
+   sigma4 ~ 1,
+   Lambda ~ -1 + s(id, by = fpc, bs = "re")
+ )
> 
> # Model
> b <- bamlss(f, family = trafficfam3, n.iter = 20, burnin = 10,
+             data = dat)
Error in a$specials[[i]] : subscript out of bounds
Calls: bamlss ... assign.design -> drop.terms.bamlss -> delete.response -> off
Execution halted

Running ‘testthat.R’ [5s/5s]
Running the tests in ‘tests/testthat.R’ failed.
Complete output:
  > # This file is part of the standard setup for testthat.
  > # It is recommended that you do not modify it.
  > #
  > # Where should you do additional test configuration?
  > # Learn more about the roles of various files in:
  > # * https://r-pkgs.org/testing-design.html#sec-tests-files-overview
  > # * https://testthat.r-lib.org/articles/special-files.html
  > 
  > library(testthat)
  > library(gmfamm)
  Loading required package: bamlss
  Loading required package: coda
  Loading required package: colorspace
  Loading required package: distributions3
  
  Attaching package: 'distributions3'
  
  The following object is masked from 'package:stats':
  
      Gamma
  
  The following object is masked from 'package:grDevices':
  
      pdf
  
  Loading required package: mgcv
  Loading required package: nlme
  This is mgcv 1.9-4. For overview type '?mgcv'.
  -
  For citation info, use citation("
...[truncated]...
mlss:::predict.bamlss(...) at gmfamm/R/predict.R:93:5
   5.       └─bamlss:::term.labels2(...)
   6.         └─bamlss:::drop.terms.bamlss(...)
   7.           └─stats::delete.response(tx)
   8.             └─stats (local) off(a$specials[[i]])
  ── Error ('test-predict.R:184:3'): prediction for single model terms works ─────
  <subscriptOutOfBoundsError/error/condition>
  Error in `a$specials[[i]]`: subscript out of bounds
  Backtrace:
      ▆
   1. ├─stats::predict(b, term = "x1", type = "link", intercept = TRUE) at test-predict.R:184:3
   2. └─bamlss:::predict.bamlss(b, term = "x1", type = "link", intercept = TRUE)
   3.   └─family$predict(...)
   4.     └─bamlss:::predict.bamlss(...) at gmfamm/R/predict.R:93:5
   5.       └─bamlss:::term.labels2(...)
   6.         └─bamlss:::drop.terms.bamlss(...)
   7.           └─stats::delete.response(tx)
   8.             └─stats (local) off(a$specials[[i]])
  
  [ FAIL 4 | WARN 0 | SKIP 0 | PASS 3 ]
  Error:
  ! Test failures.
  Execution halted