nlmixr2 / nlmixr2est

nlmixr2 estimation routines
https://nlmixr2.github.io/nlmixr2est/
GNU General Public License v2.0
7 stars 3 forks source link

`print=0` doesn't seem to work #424

Closed mattfidler closed 6 months ago

mattfidler commented 6 months ago
library(nlmixr2)
#> Loading required package: nlmixr2data

set.seed(42)

lr.mod <- function() {
  ini({
    tb0 <- 0.9 
    tb1 <- 0.1
    eta ~ 0.01
  })
  model({
    intercept<-tb0
    slope<-tb1
    p <- expit(intercept + slope * time + eta)
    n <- 1
    err ~ dbinom(n, p) 
  })
}

d <- et(10, 100, by=0.5) %>%
  et(id=1:10)

d <- rxSolve(lr.mod, d)

d <- d %>%
  dplyr::select(-DV, -intercept, -slope) %>%
  dplyr::rename(DV=sim)

f1 <- lr.mod %>%
  model(p <- expit(intercept + slope * time)) %>%
  nlmixr(d, "nlminb", control=list(print=0))
#> ! remove between subject variability `eta`
#> Error in nlminbControl(print = 0): Assertion on 'print' failed: Element 1 is not >= 1.

Created on 2024-03-11 with reprex v2.1.0

mattfidler commented 6 months ago

Neither does covMethod="" at least for nlminbControl()