Функция bargraph.CI
принимает необработанные данные, а не агрегированные данные. Поэтому попробуйте следующее:
bargraph.CI(parameter_variants, response=score, group=scale, data=dat,
main="likeability", legend=TRUE)
Это должно дать вам один "двусторонний" сюжет. Если вам не нравится его внешний вид, есть много аргументов, которые вносят поверхностные коррективы. Подробности смотрите на странице справки.
Чтобы получить отдельные графики для каждой из четырех шкал, я думаю, вы можете сделать что-то вроде этого:
library(dplyr)
dat %>%
filter(scale=="like") %>% # change the value here.
bargraph.CI(parameter_variants, response=score, data=., main="likeability")
Решение Base R:
with(subset(dat, subset=scale=="like"),
bargraph.CI(parameter_variants, response=score, main="likeability")
)