17.2 probmod package

• Not recommend: package has serious problem with subscript.

install.packages("probemod")

library(probemod)

myModel <-
    lm(loss ~ hours * gender, data = dat %>% 
           select(loss, hours, gender))
jnresults <- jn(myModel,
                dv = 'loss',
                iv = 'hours',
                mod = 'gender')


pickapoint(
    myModel,
    dv = 'loss',
    iv = 'hours',
    mod = 'gender',
    alpha = .01
)

plot(jnresults)