“Non-numeric argument to binary operator error” when making a grouped boxplot

Question

I am trying to create a grouped boxplot using a fairly simple, but long dataset called trch.

First I used trchtrue<- trch[complete.cases(trch),