Nothing
context("analytical vs numerical gradient / hessian")
test_that("gradients", {
data("Interview", package = "mhurdle")
for (form in list(shows ~ educ | linc | age, shows ~ educ | linc, shows ~ 0 | linc | age)){
for (dist in c("n", "ln")){
for (h2 in c(TRUE, FALSE)){
for (corr in c(TRUE, FALSE)){
za <- mhurdle(form, data = Interview, h2 = h2, corr = corr, dist = dist, check_gradient = TRUE)$gradient
cat(paste(deparse(form), ", dist = ", dist, ", h2 = ", h2, ", corr = ", corr, "\n"))
expect_true(all(za[, "rel_diff"] < 1E-04))
}
}
}
}
})
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.