context("Custom C++ probability routines")
test_that("Evaluating probabilities for intervals of normal r.v.s", {
set.seed(2018)
mu = rnorm(1)
sigma = runif(1)
breaks = seq(0, 1, length.out = 4)
expect_equal(
as.numeric(.Call(`_telefit_r_qintnorm`,
qnorm(breaks[-c(1, length(breaks))], mu, sigma), mu, sigma)),
diff(breaks),
tolerance = 1e-3
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.