Линейная гистограмма Лайкерта с расслоением - проблемы совместимости «заголовка» и «гистограммы»

Question

Я использую библиотеку 'likert', чтобы показать некоторые взаимосвязи между некоторыми данными likert. Ниже немного MWE:

library(likert)

mydata <- read.table(header=TRUE, text="
Q1 Q2
Always Never
Always Often
Always Sometimes
Always Sometimes
Always Always
Always Sometimes
Always Infrequently
Always Often
Always Infrequently
Always Often
Always Often
Always Sometimes
Always Often
Always Sometimes
Always Sometimes
Infrequently Never
Infrequently Never
Infrequently Sometimes
Infrequently Never
Infrequently Infrequently
Infrequently Never
Infrequently Never
Often Sometimes
Often Sometimes
Often Often
Often Sometimes
Often Sometimes
Often Infrequently
Often Infrequently
Sometimes Often
Sometimes Infrequently
Sometimes Infrequently
Sometimes Infrequently
Sometimes Sometimes
Sometimes Infrequently
")

# refining the data
groupLevels <- c("Never", "Infrequently","Sometimes","Often","Always")
mydata$Q1 <- ordered(mydata$Q1, levels = groupLevels)
mydata$Q2 <- ordered(mydata$Q2, levels = groupLevels)

# now I plot my data
plot(likert(mydata$Q2, grouping=mydata$Q1), include.histogram=TRUE, group.order=rev(groupLevels))

пока все работает нормально. Тем не менее, в тот момент, когда я пытаюсь добавить заголовок, как это plot (likert (mydata $ Q2, grouping = mydata $ Q1), include.histogram = TRUE, group.order = rev (groupLevels)) + ggtitle ("TEST")

Я получаю следующий код ошибки:

Error in plot(likert(mydata$Q2, grouping = mydata$Q1), include.histogram = TRUE,  : 
non-numeric argument to binary operator

Пока мне стало ясно, что установка заголовка работает до тех пор, пока я не использую гистограмму. Таким образом, следующее работает, а выше не работает:

    plot(likert(mydata$Q2, grouping=mydata$Q1), group.order=rev(groupLevels)) + ggtitle("TEST")

У кого-нибудь есть идея, как я могу использовать гистограмму и одновременно добавить заголовок к диаграмме?