# NOT RUN {
# Using a fitted model as formula input
fiti <- lm(Income ~ Frost + Murder*Illiteracy,
data=as.data.frame(state.x77))
probe_interaction(model = fiti, pred = Murder, modx = Illiteracy,
modxvals = "plus-minus")
# 3-way interaction
fiti3 <- lm(Income ~ Frost*Murder*Illiteracy,
data=as.data.frame(state.x77))
probe_interaction(model = fiti3, pred = Murder, modx = Illiteracy,
mod2 = Frost, mod2vals = "plus-minus")
# With svyglm
library(survey)
data(api)
dstrat <- svydesign(id=~1,strata=~stype, weights=~pw, data=apistrat, fpc=~fpc)
regmodel <- svyglm(api00~ell*meals + sch.wide, design = dstrat)
probe_interaction(model = regmodel, pred = ell, modx = meals,
modxvals = "plus-minus", cond.int = TRUE) # conditional intercepts
# 3-way with survey and factor input
regmodel3 <- svyglm(api00~ell*meals*sch.wide, design = dstrat)
probe_interaction(model = regmodel3, pred = ell, modx = meals, mod2 = sch.wide)
# Can try different configurations of 1st vs 2nd mod
probe_interaction(model = regmodel3, pred = ell, modx = sch.wide, mod2 = meals)
# }
Run the code above in your browser using DataLab