Как построить одиночную числовую ковариату, используя emmeans (или другой пакет) из модели?

Использование sjPlot :

plot_model(m1, terms = "disp [all]", type = "pred")

дает ту же графику.

Использование emmeans :

em1 <- ref_grid(m1, at = list(disp = seq(min(mtcars$disp), max(mtcars$disp), 1)))
emmip(em1, ~disp, CIs = T)

возвращает графику с небольшой разницей в макете. Альтернативой является добавление результата к объекту и построение графика так, как я хочу:

d1 <- emmip(em1, ~disp, CIs = T, plotit = F)

Я подозреваю, что ваша проблема связана с тем, что по умолчанию ковариаты сводятся к их среднему значению в emmeans . Вы можете использовать аргументы at или cov.reduce, чтобы указать большее количество значений. См. Документацию для ref_grid и vignette(“basics”, “emmeans”) или указатель тем виньетки.