tests/testthat/test_polya_marginal.r
In kylebittinger/polyafit: Fit observed data to a multivariate Polya distribution
test_that("marginal probability of each category is beta binomial", {
observed <- ppolya_marginal(c(1, 5, 6), c(8, 15, 16))
expect_equal(observed[1], pbetabinom(1, 12, 8, 31)) # 15 + 16 = 31
expect_equal(observed[2], pbetabinom(5, 12, 15, 24)) # 8 + 16 = 24
expect_equal(observed[3], pbetabinom(6, 12, 16, 23)) # 8 + 15 = 23
})
kylebittinger/polyafit documentation built on Jan. 11, 2023, 8:53 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
test_that("marginal probability of each category is beta binomial", {
observed <- ppolya_marginal(c(1, 5, 6), c(8, 15, 16))
expect_equal(observed[1], pbetabinom(1, 12, 8, 31)) # 15 + 16 = 31
expect_equal(observed[2], pbetabinom(5, 12, 15, 24)) # 8 + 16 = 24
expect_equal(observed[3], pbetabinom(6, 12, 16, 23)) # 8 + 15 = 23
})
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.