Как переупорядочить факторный термин в plot_model при отображении эффектов упорядоченной логистической регрессии (функция polr в MASS))

Question

Я строю эффекты упорядоченной логистической регрессии, используя функцию sjplot plot_model.

код для подгонки модели:

mult_4_rank_best <- polr(SpRank ~ RodentRank + SpName + day_start_scale + Lat_scale +
                           Lon2_scale + Lon_scale + RodentRank:SpName + 
                           day_start_scale:SpName + Lon2_scale:SpName + 
                           Lon_scale:SpName + RodentRank:SpName:day_start_scale, 
                         data = data_mult,
                         na.action = na.pass, 
                         Hess = TRUE)

Несмотря на то, что есть много предикторов, я хочу увидеть влияние Rodent Rank и SpName на результат и использовать этот код:

plot_model(mult_4_rank_best, type = "eff", terms = c("RodentRank", "SpName"))

Несмотря на то, что переменные в наборе данных были упорядочены таким образом до подгонки модели

data_mult$RodentRank <- factor(data_mult$RodentRank, 
                               levels = c("low", "average", "high"), 
                               ordered=TRUE)

Я продолжаю получать итоговую картину, где уровни RodentRank расположены в алфавитном порядке. order.terms = c() не помогает.

Я не мог гуглить ответ на эту проблему. Если кто-то знает, как с этим справиться, я был бы очень признателен.