library(pequod)

model1 <- lmres(sat~talk*per, centered =c("talk", "per"), data =dat)
summary(model1)

model2 <- simpleSlope(model1, pred="talk", mod1 = "per")
summary(model2)

PlotSlope(model2)

3要因の重回帰分析

model1 <- lmres(sat~talk*per*con,centered=c("talk","per","con"),data=dat)


model2 <- simpleSlope(model1,pred="talk",mod1 = "per",mod2="con")