Nothing
context("Uniform scaled beta (USB) distribution")
set.seed(5)
sq <- seq(-1, 1, by = 0.001)
test_that("Density integrates to one", {
expect_equal(integrate(dusb, -1, 1)$value, 1)
expect_equal(integrate(dusb, -1, 1, theta = 0.4)$value, 1)
expect_equal(integrate(dusb, -1, 1, theta = 0.4, nu = 10)$value, 1)
})
test_that("Simulation works", {
expect_lt(max(abs(rusb(500))), 0.5)
expect_gt(max(abs(rusb(500))), 0.48)
expect_lt(max(abs(rusb(500, theta = 0.2))), 0.6)
expect_gt(max(abs(rusb(500, theta = 0.1))), 0.5)
expect_lt(max(abs(rusb(500, theta = 0.99, nu = 10))), 1)
expect_gt(max(abs(rusb(500, theta = 0.99, nu = 1))), 0.94)
expect_lt(max(abs(rusb(500, quasi = TRUE))), 0.5)
expect_gt(max(abs(rusb(500, quasi = TRUE))), 0.48)
expect_lt(max(abs(rusb(500, quasi = TRUE, theta = 0.2))), 0.6)
expect_gt(max(abs(rusb(500, quasi = TRUE, theta = 0.1))), 0.5)
expect_lt(max(abs(rusb(500, quasi = TRUE, theta = 0.99, nu = 10))), 1)
expect_gt(max(abs(rusb(500, quasi = TRUE, theta = 0.99, nu = 1))), 0.94)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.