context("Logit Integrand")
test_that("Value of Logit Integrand equals", {
# Check it works when n_i = 1
expect_equal(logit_integrand(0,
X = cbind(1, .5),
A = 1,
parameters = c(.1, .1),
randomization = .5), 0.107201744028612)
# No random effect
expect_equal(logit_integrand(0,
X = cbind(rep(1, 2), c(.5, .6)),
A = c(0, 1),
parameters = c(.1, .1),
randomization = .5), 0.078757529007955)
# With random effect
expect_equal(logit_integrand(c(0, 1),
X = cbind(rep(1, 2), c(.5, .6)),
A = c(0, 1),
parameters = c(.1, .1, 1),
randomization = .5), c(0.0787575290079550,
0.0571307956171674))
# Random effect with allocation within integrand
expect_equal(logit_integrand(c(0, 1),
X = cbind(rep(1, 2), c(.5, .6)),
A = c(0, 1),
parameters = c(.1, .1, 1),
randomization = .5,
allocation = .5 ), c(0.315030116031820,
0.228523182468669))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.