library(regsim)

data(mtcars)
mtcars$cyl <- as.factor(mtcars$cyl)
model <- lm(mpg ~ wt + disp + cyl, mtcars)
summary(model)
set.seed(666)
x <- list(wt = seq(1, 5), cyl = levels(mtcars$cyl))

sim <- regsim(model, x)
summary(sim)
plot(sim, ~ wt + cyl, lines.col = "Set1", lines.lwd = 1.5)


altaf-ali/regsim documentation built on May 9, 2019, 9:57 p.m.