Skip to content

Commit 40fd4dd

Browse files
committedMay 9, 2023
Add to vignette
1 parent 97cc183 commit 40fd4dd

File tree

1 file changed

+11
-0
lines changed

1 file changed

+11
-0
lines changed
 

‎vignettes/model-description.Rmd

+11
Original file line numberDiff line numberDiff line change
@@ -548,6 +548,17 @@ The mean is $(1-p) \cdot \mu_{1} + p \cdot \mu_{2}$ and the variance is $(1-p) ^
548548
Here, and for the other mixture distributions, the probability of the larger mean can be obtained from `plogis(fit$model$par[["logit_p_mix"]])` and the ratio of the larger mean to the smaller mean can be obtained from `1 + exp(fit$model$par[["log_ratio_mix"]])`.
549549
The standard errors are available in the TMB sdreport: `fit$sd_report`.
550550
551+
If you wish to fix the probability of a large (i.e., extreme) mean, which can be hard to estimate, you can use the `map` list. E.g.:
552+
553+
```r
554+
sdmTMB(...,
555+
control = sdmTMBcontrol(
556+
start = list(logit_p_mix = qlogis(0.05)), # 5% probability of 'mu2'
557+
map = list(logit_p_mix = factor(NA)) # don't estimate
558+
)
559+
)
560+
```
561+
551562
Example: `family = gamma_mix(link = "log")`. See also `family = delta_gamma_mix()` for an extension incorporating this distribution with delta models.
552563
553564
## Lognormal mixture

0 commit comments

Comments
 (0)
Please sign in to comment.