使用两个连续预测变量可视化线性混合效应模型

Question

我从以下线性混合效果模型中获得了重要的互动：

library(lmer4)
library(ggplot2)
library(lmerTest)
library(effects)
library(extrafont)

mod_father_son <- lmer(AIP_s_child.z ~ AIP_s_parent.z*Q_mean.z + 
                  (1 + AIP_s_parent.z:Q_mean.z || Family_number), 
                data = data_father_son, REML = FALSE)

我的数据与此类似（Q_mean和AIP_s是z分数）：

Id_parent   Family_number   AIP_s_parent.z   Q_mean.z   Child_id   AIP_s_child.z

A1          1               -.008            -0.5       B1         .005
A1          1               -.008            -0.5       B2         .04
C1          2                .06             -.006      D1         -.007
E1          3               -.1              0.02       F1         -.06

我希望以最佳方式将其可视化，我正在考虑一个线图，显示Q_mean.z的+1,0，-1处的交互。我当前图表的代码是：

fs <- as.data.frame(Effect(c("AIP_s_parent.z","Q_mean.z"),mod = mod_father_son))

fs <- mutate(fs, Q_mean.z = as.character(Q_mean.z))

plot_fs <- ggplot(fs, aes(AIP_s_parent.z,fit, group=Q_mean.z, linetype = Q_mean.z))+
 geom_line() + 
 labs(x = "Gender stereotypes father", y = "Gender stereotypes son", linetype = "Questionnaire score")

plot_fs + theme_classic() +
 scale_color_gradient(low = "black", high = "gray90") +
 theme(text=element_text(size=12,  family="Times New Roman"))

图表如下所示：

这是否正确，因为图表的解释对模型的输出没有意义？有人告诉我，我应该使用带有ggplot的predict（），但我不知道该怎么做。另外，我想在我的情节中只有三行（Q_mean.z的+1,0，-1），因为五行令人困惑且难以阅读。如果折线图不是可视化交互的最佳方式，我愿意接受建议。