model_select.melsm: Select Model for 'melsm' Objects

Description Usage Arguments Value Examples

View source: R/model_select.melsm.R

Description

Model selection for mixed effects location scale models.

Usage

1
2
## S3 method for class 'melsm'
model_select(object, type = "hpm", ...)

Arguments

object

melsm object

type

model to select (currently the only option is the highest posterior model hpm)

...

currently ignored

Value

object of class model_select.melsm

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
## Not run: 
dat <- nlme::Orthodont
fit <- melsm(fixed_location = distance ~  age,
            random_location = ~ age | Subject,
            fixed_scale = sigma ~  1, k = 3,
            random_scale = ~ 1 | Subject,
            adapt = 5000,
            iter = 10000,
            rho_test = "all",
            mixture = "SSVS",
            data = dat)

# highest posterior model
hpm <- model_select(fit, type = "hpm")
summary(hpm, cred = 0.95)

## End(Not run)

donaldRwilliams/hypMuVar documentation built on Jan. 10, 2020, 9:45 a.m.