Ошибка модели GAM

Мой фрейм данных выглядит так:

head(bush_status)
distance  status count
       0 endemic   844
       1 exotic     8
       5 native     3
      10 endemic    5
      15 endemic    4
      20 endemic    3

Данные подсчета распределены ненормально. Я пытаюсь подогнать обобщенную аддитивную модель к своим данным двумя способами, поэтому я могу использовать anova, чтобы увидеть, поддерживает ли p-значение m2.

m1 <- gam(count ~ s(distance) + status, data=bush_status, family="nb")
m2 <- gam(count ~ s(distance, by=status) + status, data=bush_status, family="nb")

m1 работает нормально, но m2 отправляет сообщение об ошибке:

"Error in smoothCon(split$smooth.spec[[i]], data, knots, absorb.cons, 
scale.penalty = scale.penalty,  : 
  Can't find by variable"

Это не для меня, поэтому, если бы кто-нибудь мог дать какой-либо совет, был бы очень признателен!