test_that("binomial.sns.test works", {
test.function <- binomial.sns.test
# This is the example given in this ms:
# http://papers.ssrn.com/sol3/papers.cfm?abstract_id=2256237
p1 <- c(0.0, 0.0, 0.0, 0.0, 0.0)
expect_that(test.function(p1)$'higher', equals(0.0))
expect_that(test.function(p1)$'lower', equals(5.0))
expect_that(test.function(p1)$'p.rskew', equals(0.03125))
expect_that(test.function(p1)$'p.lskew', equals(1.00000))
# This is the counter example, just for completeness
p2 <- c(0.05, 0.05, 0.05, 0.05, 0.05)
expect_that(test.function(p2)$'higher', equals(5.0))
expect_that(test.function(p2)$'lower', equals(0.0))
expect_that(test.function(p2)$'p.rskew', equals(1.00000))
expect_that(test.function(p2)$'p.lskew', equals(0.03125))
# This one should return 1, because no data
p3 <- c(0.025, 0.025, 0.025, 0.025, 0.025, 0.025)
expect_that(test.function(p3)$'higher', equals(0.0))
expect_that(test.function(p3)$'lower', equals(0.0))
expect_that(test.function(p3)$'p.rskew', equals(1.0000))
expect_that(test.function(p3)$'p.lskew', equals(1.0000))
})
test_that("binomial.all.test works", {
test.function <- binomial.all.test
p1 <- c(0.0, 0.0, 0.0, 0.0, 0.0)
expect_that(test.function(p1)$'higher', equals(0.0))
expect_that(test.function(p1)$'lower', equals(5.0))
expect_that(test.function(p1)$'p', equals(0.0625))
# This is the counter example, just for completeness
p2 <- c(0.05, 0.05, 0.05, 0.05, 0.05)
expect_that(test.function(p2)$'higher', equals(5.0))
expect_that(test.function(p2)$'lower', equals(0.0))
expect_that(test.function(p2)$'p', equals(0.0625))
# This one should return 1, because no data
p3 <- c(0.025, 0.025, 0.025, 0.025, 0.025, 0.025)
expect_that(test.function(p3)$'higher', equals(0.0))
expect_that(test.function(p3)$'lower', equals(0.0))
expect_that(test.function(p3)$'p', equals(1))
})
test_that("binomial.bias.test works", {
test.function <- binomial.bias.test
p1 <- c(0.0, 0.0, 0.0, 0.0, 0.0)
expect_that(test.function(p1)$'lower', throws_error())
expect_that(test.function(p1)$'higher', throws_error())
expect_that(test.function(p1)$'p', throws_error())
p2 <- c(0.05, 0.05, 0.05, 0.05, 0.05)
expect_that(test.function(p2)$'lower', equals(0.0))
expect_that(test.function(p2)$'higher', equals(5.0))
expect_that(test.function(p2)$'p', equals(0.03125))
p3 <- c(0.05, 0.05, 0.05, 0.05, 0.039, 0.05, 0.04001, 0.041, 0.040)
expect_that(test.function(p3)$'lower', equals(1.0))
expect_that(test.function(p3)$'higher', equals(7.0))
expect_that(test.function(p3)$'p', equals(0.03515625))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.