I am running a multilevel model with two-levels with one interaction term at the within-level.
I need to present two plots: a) simple slopes with slopes at -2SD, mean, +2SD of the moderator, with my x variable on the x-axis and three slopes (M) and b) a Johnson-Neyman (limits -2SD, +2SD)
From what I've read from examples and the forum, I have added to my MPLUS script:
MODEL CONSTRAINT: Plot(ss); Loop (mod,-2,2,0.1); ss = a+b*mod;
This script gives me a gh5 file for R (I am working on a Mac), but when I put it into R, it produces a plot with one slope plus LCI/UCI with the moderator on the X-axis. I think that it's an error in my Mplus part of the process but am not sure.