Ошибка при вычислении Nested ANOVA R

Question

мои данные вложены на уровне группы. Есть пять разных процедур. В рамках каждого лечения четыре участника группируются. Речь идет о поведении доноров тех участников (зависимая переменная = пожертвование, метри c, в евро) во время соревнования (пояснительная переменная = лечение, порядковый номер). Структура данных выглядит следующим образом:

Treatment   Group   Donation
1             1         0,6
1             1         0 
1             1         0,3
1             1         1,2
1             2         0,1
1             2         0,5
1             2         0,1
1             2         0,8
2             1         0,4 
...

Я хочу запустить вложенную ANOVA на основе примера на этой странице: Rcompanion . Структура должна быть очень похожей: Лечение = Технология, Группа = Крыса, Донорство = Белок. Однако, если я пытаюсь запустить команду с моделью:

> model = lme(Donation ~ Treatment,
+             random = ~1|Group,
+             data=R_Data,
+             method = "REML")

, я получаю это сообщение об ошибке:

Error in `rownames<-`(`*tmp*`, value = rownames(Fitted) <- origOrder) : 
  attempt to set 'rownames' on an object with no dimensions
In addition: Warning message:
In Ops.factor(y[revOrder], Fitted) : ‘-’ not meaningful for factors

Я совершенно растерялся, что я делаю не так?