context("test-misc fixest")
test_that("weird regression model with r^2",
{
mean = rnorm(10)
error = rnorm(10, sd = 0.3)
y = mean + error
constant = rep(2, 10)
group = as.factor(c(rep(1:5, 2)))
factor = as.factor(1:10)
# multicol_model = lm(y ~ group + constant)
# expect_warning(partial_r2(multicol_model))
data = model.matrix(y ~ factor)[,-1]
data = cbind(y, data)
fml <- reformulate(paste(colnames(data)[-1], collapse = " + "), response = "y")
saturated_model = fixest::feols(fml = fml, data = data)
expect_equal(TRUE, all(is.nan(partial_r2(saturated_model))))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.