data(cats, package = "MASS")
l1 <- linmod(Hwt ~ Bwt * Sex, data = cats)
l2 <- lm(Hwt ~ Bwt * Sex, data = cats)
test_that("same estimated coefficients as lm function", {
expect_equal(round(l1$coefficients, 3), round(l2$coefficients, 3))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.