confirm.melsm: Confirmatory Hypothesis Testing for 'melsm' Objects

In donaldRwilliams/hypMuVar: Bayesian Hypothesis Testing of Mean-Variance Relations

Description

Confirmatory hypotheses testing of scientific expectations.

Usage

## S3 method for class 'melsm'
confirm(object, hyp, ...)

Arguments

object	object of class melsm
hyp	hypothesis to test
...	currently ignored

Value

object of class confirm

• BF Bayes factor

• post_probs posterior probabilities

• prior_probs prior probabilities

• model model type (melsm)

• mixture mixture type

Examples

## 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)


# hypothesis (rho_01 is null and rho_02 is negative )
# vs compliment (not h1)
hyp <- list(h1 = c("rho_01 = k1", "rho_02 = k3"))
confirm(fit, hyp)

# hypothesis  (both positive) vs both null
hyp <- list(h1 = c("rho_01 = k2", "rho_02 = k2"),
h2 = c("rho_01 = k1", "rho_02 = k1"))
confirm(fit, hyp)


## End(Not run)

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