diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2019-02-26 20:28:15 +0100 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2019-02-26 20:33:32 +0100 |
| commit | 396e40d2c45d81a8e74c41a2f632665021e980a1 (patch) | |
| tree | c7ab9b5356361b1939e9b00183a7f31071782f64 /R/mkinfit.R | |
| parent | 64e636d11311be57978ef394c48c9b43fe219059 (diff) | |
NAFTA SOP Attachment vignette
Also:
- Change rounding in print.nafta
- Add dots argument to nafta()
- Use cores=1 in examples
- Restrict N in IORE model to values > 0
Diffstat (limited to 'R/mkinfit.R')
| -rw-r--r-- | R/mkinfit.R | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/R/mkinfit.R b/R/mkinfit.R index 2dec3a72..63455300 100644 --- a/R/mkinfit.R +++ b/R/mkinfit.R @@ -374,6 +374,10 @@ mkinfit <- function(mkinmod, observed, lower <- rep(-Inf, length(c(state.ini.optim, transparms.optim)))
upper <- rep(Inf, length(c(state.ini.optim, transparms.optim)))
names(lower) <- names(upper) <- c(names(state.ini.optim), names(transparms.optim))
+
+ # IORE exponentes are not transformed, but need a lower bound of zero
+ index_N <- grep("^N", names(lower))
+ lower[index_N] <- 0
if (!transform_rates) {
index_k <- grep("^k_", names(lower))
lower[index_k] <- 0
|