Resolving the '‘==’ only defined for equally-sized data frames' Error in Generalized Additive Models with gratia in R