Nothing
test_that("hfr", {
data("Boston", package = "MASS")
x <- model.matrix(medv ~ 0 + . + crim*tax, Boston)
y <- Boston$medv
expect_no_error(mod_1 <- hfr(x, y, kappa = 1))
expect_no_error(mod_2 <- hfr(x, y, kappa = 0))
expect_no_error(mod_3 <- hfr(x, y, kappa = 0.5))
expect_no_error(coef(mod_1))
expect_no_error(coef(mod_2))
expect_no_error(coef(mod_3))
expect_no_error(predict(mod_1, x))
expect_no_error(predict(mod_2, x))
expect_no_error(predict(mod_3, x))
expect_no_error(fitted(mod_1))
expect_no_error(fitted(mod_2))
expect_no_error(fitted(mod_3))
expect_no_error(se.avg(mod_1))
expect_no_error(se.avg(mod_2))
expect_no_error(se.avg(mod_3))
expect_no_error(plot(mod_1))
expect_no_error(plot(mod_3))
expect_false(any(is.na(coef(mod_1))))
expect_false(any(is.na(coef(mod_2))))
expect_false(any(is.na(coef(mod_3))))
expect_false(any(is.na(predict(mod_1, x))))
expect_false(any(is.na(predict(mod_2, x))))
expect_false(any(is.na(predict(mod_3, x))))
expect_false(any(is.na(fitted(mod_1))))
expect_false(any(is.na(fitted(mod_2))))
expect_false(any(is.na(fitted(mod_3))))
expect_false(any(is.na(se.avg(mod_1))))
expect_false(any(is.na(se.avg(mod_2))))
expect_false(any(is.na(se.avg(mod_3))))
})
test_that("hfr-linear-dep", {
data("Boston", package = "MASS")
x <- model.matrix(medv ~ 0 + . + crim*tax, Boston)
y <- Boston$medv
x[,1] <- x[,2]
expect_warning(mod_1 <- hfr(x, y, kappa = 1))
expect_warning(mod_2 <- hfr(x, y, kappa = 0))
expect_warning(mod_3 <- hfr(x, y, kappa = 0.5))
expect_no_error(coef(mod_1))
expect_no_error(coef(mod_2))
expect_no_error(coef(mod_3))
expect_no_error(predict(mod_1, x))
expect_no_error(predict(mod_2, x))
expect_no_error(predict(mod_3, x))
expect_no_error(fitted(mod_1))
expect_no_error(fitted(mod_2))
expect_no_error(fitted(mod_3))
expect_no_error(se.avg(mod_1))
expect_no_error(se.avg(mod_2))
expect_no_error(se.avg(mod_3))
expect_no_error(plot(mod_1))
expect_no_error(plot(mod_3))
expect_false(any(is.na(coef(mod_1))))
expect_false(any(is.na(coef(mod_2))))
expect_false(any(is.na(coef(mod_3))))
expect_false(any(is.na(predict(mod_1, x))))
expect_false(any(is.na(predict(mod_2, x))))
expect_false(any(is.na(predict(mod_3, x))))
expect_false(any(is.na(fitted(mod_1))))
expect_false(any(is.na(fitted(mod_2))))
expect_false(any(is.na(fitted(mod_3))))
expect_false(any(is.na(se.avg(mod_1))))
expect_false(any(is.na(se.avg(mod_2))))
expect_false(any(is.na(se.avg(mod_3))))
})
test_that("hfr-const-col", {
data("Boston", package = "MASS")
x <- model.matrix(medv ~ 0 + . + crim*tax, Boston)
y <- Boston$medv
x[,1] <- 1
expect_no_error(mod_1 <- hfr(x, y, kappa = 1))
expect_no_error(mod_2 <- hfr(x, y, kappa = 0))
expect_no_error(mod_3 <- hfr(x, y, kappa = 0.5))
expect_no_error(coef(mod_1))
expect_no_error(coef(mod_2))
expect_no_error(coef(mod_3))
expect_no_error(predict(mod_1, x))
expect_no_error(predict(mod_2, x))
expect_no_error(predict(mod_3, x))
expect_no_error(fitted(mod_1))
expect_no_error(fitted(mod_2))
expect_no_error(fitted(mod_3))
expect_no_error(se.avg(mod_1))
expect_no_error(se.avg(mod_2))
expect_no_error(se.avg(mod_3))
expect_warning(plot(mod_1))
expect_warning(plot(mod_3))
expect_true(any(is.na(coef(mod_1))))
expect_true(any(is.na(coef(mod_2))))
expect_true(any(is.na(coef(mod_3))))
expect_false(any(is.na(predict(mod_1, x))))
expect_false(any(is.na(predict(mod_2, x))))
expect_false(any(is.na(predict(mod_3, x))))
expect_false(any(is.na(fitted(mod_1))))
expect_false(any(is.na(fitted(mod_2))))
expect_false(any(is.na(fitted(mod_3))))
expect_true(any(is.na(se.avg(mod_1))))
expect_true(any(is.na(se.avg(mod_2))))
expect_true(any(is.na(se.avg(mod_3))))
})
test_that("hfr-high-dim", {
x <- matrix(rnorm(10*15), 10)
y <- rnorm(10)
expect_warning(mod_1 <- hfr(x, y, kappa = 1))
expect_warning(mod_2 <- hfr(x, y, kappa = 0))
expect_warning(mod_3 <- hfr(x, y, kappa = 0.5))
expect_no_error(coef(mod_1))
expect_no_error(coef(mod_2))
expect_no_error(coef(mod_3))
expect_no_error(predict(mod_1, x))
expect_no_error(predict(mod_2, x))
expect_no_error(predict(mod_3, x))
expect_no_error(fitted(mod_1))
expect_no_error(fitted(mod_2))
expect_no_error(fitted(mod_3))
expect_no_error(se.avg(mod_1))
expect_no_error(se.avg(mod_2))
expect_no_error(se.avg(mod_3))
expect_no_error(plot(mod_1))
expect_no_error(plot(mod_3))
expect_false(any(is.na(coef(mod_1))))
expect_false(any(is.na(coef(mod_2))))
expect_false(any(is.na(coef(mod_3))))
expect_false(any(is.na(predict(mod_1, x))))
expect_false(any(is.na(predict(mod_2, x))))
expect_false(any(is.na(predict(mod_3, x))))
expect_false(any(is.na(fitted(mod_1))))
expect_false(any(is.na(fitted(mod_2))))
expect_false(any(is.na(fitted(mod_3))))
expect_true(any(is.na(se.avg(mod_1))))
expect_true(any(is.na(se.avg(mod_2))))
expect_true(any(is.na(se.avg(mod_3))))
})
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.