dat <- data_test1
lm1 <- lm(y ~ x1 + x2, dat)
lm2 <- lm(y ~ x1 + x2 + x3 + x4, dat)
lm3 <- lm(y ~ x1 + cat1 + cat2 + x2 + x3 + x4, dat)
lm4 <- lm(y ~ x1 + x2*x3 + x4, dat)
hierarchical_lm(lm1, lm3, lm2)
hierarchical_lm(lm1, lm2, lm4)
Run the code above in your browser using DataLab