I am desperately trying to plot a Johnson-Neyman Plot for the following interaction:
nxsc_20 <-lm(meandec20 ~ centered_nep*centered_selfcontrol + factor(study)
