Запустите ANOVA с подгруппами в R

Question

Я пытаюсь вычислить вложенный ANOVA, используя R. Для этого у меня есть Участники, которые жертвуют деньги на разных Сессиях на основе определенного c Обработки. Мои данные выглядят следующим образом:

Treatment     Session     Participant
3             1           0.5
3             1           1.0
3             1           1.1
3             2           0
3             2           0.7
3             2           1.6
4             1           0.6
4             1           0.4
4             1           0.4
4             2           0
4             2           0.6
4             2           0.2

Мой подход:

> Data$Session = as.factor(Data$Session)
> library(nlme)
> 
> model = lme(Participant ~ Treatment, random=~1|Session,
+             method="REML")
> anova.lme(model,
+           type="sequential",
+           adjustSigma = FALSE)

Но когда я запускаю это, я получаю такой вывод:

> model = lme(Participant ~ Treatment, random=~1|Session,
+             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
> 
> anova.lme(model,
+           type="sequential",
+           adjustSigma = FALSE)
Error in anova.lme(model, type = "sequential", adjustSigma = FALSE) : 
  Object has to inherit from class "lme". 

Есть ли какая-либо ошибка в рассуждениях с моей стороны? Заранее спасибо!