#test for variances
vtest1 <- bartlett_test(InsectSprays$count, InsectSprays$spray)
hypothesis <- "A=B=F>C=D=E"
set.seed(123)
BF1 <- BF(x=vtest1,hypothesis)
BF11 <- BF(x=vtest1,prior.hyp.explo = 3:4)
#check results exploratory test
test_that("BF.bartlett_htest exploratory hypotheses correctly evaluated", {
expect_equivalent(
round(BF1$PHP_exploratory,7),c(0.0044175,0.9955825)
)
expect_equivalent(
unname(round(BF11$PHP_exploratory,6)),
unname(round(BF1$BFtu_exploratory * (3:4) / sum(BF1$BFtu_exploratory * (3:4)),6))
)
})
#check results confirmatory test
test_that("BF.bartlett_htest confirmatory hypotheses correctly evaluated", {
expect_equivalent(
round(BF1$PHP_confirmatory,7),c(0.9911905,0.0088095)
)})
BF1a <- BF(x=vtest1,hypothesis,log=TRUE)
test_that("BF.bartlett_htest exploratory hypotheses correctly evaluated on log scale", {
expect_equivalent(
round(exp(BF1a$BFtu_exploratory),5),round(BF1$BFtu_exploratory,5)
)})
hypothesis <- "A=B=F>C=D=E; A=B=F>C>D>E"
set.seed(123)
BF1 <- BF(x=vtest1,hypothesis,complement = F, log = TRUE)
#check results confirmatory test
test_that("BF.bartlett_htest confirmatory hypotheses correctly evaluated log(BF)", {
expect_equivalent(
round(BF1$BFtu_confirmatory,3),c(4.723,4.245)
)})
test_that("BF.bartlett_htest confirmatory hypotheses correctly evaluated log(BF)", {
expect_equivalent(
round(BF1$BFtu_exploratory[1],4),-5.4178
)})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.