含交互项的有序回归模型结果效应图绘制求助

我的论文已被期刊接收，但审稿人要求基于含交互项的有序回归模型绘制效应图，我从未做过这类操作，不知道怎么实现，求建议或指引方向。以下是我运行有序回归的代码及输出：

library(ordinal)
library(rcompanion) 
library(MASS) 
library(brant) 
library(AER)
library(car)

options(scipen = 999)

Ordinaldf <- read.delim(file.choose()) 
View(Ordinaldf)

Ordinaldf$AoA_scaled <- c(scale(Ordinaldf$AoA))
Ordinaldf$Ra_scaled <- c(scale(Ordinaldf$Range))
Ordinaldf$Fr_scaled <- c(scale(Ordinaldf$Frequency))
Ordinaldf$Div_scaled <- c(scale(Ordinaldf$Diversity))

Ordinaldf$Grade<-as.factor(Ordinaldf$Grade) 
Ordinaldf$Text<-as.factor(Ordinaldf$Text)

modelnull <- clm(as.factor(Grade) ~ 1,
                 data = Ordinaldf,
                 link = "logit")

model1 <- clm(Grade ~ Text *(AoA_scaled + Ra_scaled + Fr_scaled + Div_scaled),
              data = Ordinaldf,
              link = "logit")

anova(modelnull, model1)
nagelkerke(fit = model1, null = modelnull)  
summary(model1)
confint(model1)
exp(coef(model1))
exp(confint(model1))

model1 <- clm(Grade ~ Text + AoA_scaled + Ra_scaled + Fr_scaled + Div_scaled + 
                Text*AoA_scaled + Text*Ra_scaled + Text*Fr_scaled + Text*Div_scaled, 
              data = Ordinaldf,
              link = "logit")

模型输出结果：

link  threshold nobs logLik  AIC    niter max.grad cond.H 
     logit flexible  165  -261.80 551.60 5(0)  1.51e-12 6.9e+01

Coefficients:
                        Estimate Std. Error z value  Pr(>|z|)
TextTeacher            -0.008491   0.284899  -0.030    0.9762
AoA_scaled             -1.374471   0.311885  -4.407 0.0000105
Ra_scaled              -0.039952   0.318091  -0.126    0.9000
Fr_scaled              -0.056617   0.307116  -0.184    0.8537
Div_scaled             -0.367834   0.259090  -1.420    0.1557
TextTeacher:AoA_scaled  0.772882   0.375991   2.056    0.0398
TextTeacher:Ra_scaled  -0.065245   0.428996  -0.152    0.8791
TextTeacher:Fr_scaled  -0.116927   0.418322  -0.280    0.7799
TextTeacher:Div_scaled -0.343463   0.356337  -0.964    0.3351
                       
TextTeacher               
AoA_scaled             ***
Ra_scaled                 
Fr_scaled                 
Div_scaled                
TextTeacher:AoA_scaled *  
TextTeacher:Ra_scaled     
TextTeacher:Fr_scaled     
TextTeacher:Div_scaled    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Threshold coefficients:
    Estimate Std. Error z value
A|B -2.12352    0.29146  -7.286
B|C -1.09145    0.24938  -4.377
C|D -0.05468    0.23512  -0.233
D|E  0.89250    0.24752   3.606
E|F  2.36155    0.30808   7.665

内容的提问来源于stack exchange，提问作者Christian