Jamovi and R plotting different graphs
Posted: Wed Jun 12, 2019 1:59 pm
Hi,
I am trying to conduct a linear mixed effects analysis in Jamovi. Here is the code for the model that I am fitting in R.
sad.df.model13 <- lmer(response_time ~ language*target_intensity*z_AoA + task_order*language + language*target_intensity*z_PPVT
+ (1|subject) + (1|face), data = df_sad, REML = FALSE,
control= lmerControl(optimizer = "bobyqa"))
language, target_intensity and task_order are all factors with two levels. z_PPVT and z_AoA are centered continuous variables.
The interaction between language and z_AoA is significant so I plotted the interaction with interact_plot function to visualize the relationship.
interact_plot(sad.df.model13, pred = z_AoA, modx = language, interval = TRUE,
colors = "Set2", x.label = "Age of Acquisition", y.label = "Reaction Time (msecs)",
legend.main = "Language")
I did the same analysis with Jamovi and all of the model results and parameter estimates are the same.
However, for some reason, when I tried to plot the same interaction in Jamovi, the resulting graph is totally different than the one in R (even though the parameter estimate for the interaction is the same)
So I am wondering what might the underlying problem? I am attaching both graphs as well as my data.
Thanks in advance!
I am trying to conduct a linear mixed effects analysis in Jamovi. Here is the code for the model that I am fitting in R.
sad.df.model13 <- lmer(response_time ~ language*target_intensity*z_AoA + task_order*language + language*target_intensity*z_PPVT
+ (1|subject) + (1|face), data = df_sad, REML = FALSE,
control= lmerControl(optimizer = "bobyqa"))
language, target_intensity and task_order are all factors with two levels. z_PPVT and z_AoA are centered continuous variables.
The interaction between language and z_AoA is significant so I plotted the interaction with interact_plot function to visualize the relationship.
interact_plot(sad.df.model13, pred = z_AoA, modx = language, interval = TRUE,
colors = "Set2", x.label = "Age of Acquisition", y.label = "Reaction Time (msecs)",
legend.main = "Language")
I did the same analysis with Jamovi and all of the model results and parameter estimates are the same.
However, for some reason, when I tried to plot the same interaction in Jamovi, the resulting graph is totally different than the one in R (even though the parameter estimate for the interaction is the same)
So I am wondering what might the underlying problem? I am attaching both graphs as well as my data.
Thanks in advance!