我想做一个方差分析来了解哪里有意义。我已经回答了我的问题,但没有发现错误。
名:
[1] "Tier_ID" "species" "Klima" "Ressource" "Datum" "Gewicht" "IngestionRate"
数据框:
'data.frame': 70 obs. of 7 variables:
$ Tier_ID : Factor w/ 70 levels "Raupe1","Raupe10",..: 1 12 23 34 45 56 67 69 70 2 ...
$ species : Factor w/ 1 level "Agrotis exclamationis": 1 1 1 1 1 1 1 1 1 1 ...
$ Klima : Factor w/ 1 level "BL4": 1 1 1 1 1 1 1 1 1 1 ...
$ Ressource : Factor w/ 3 levels "Kontrolle","N",..: 2 2 2 2 2 2 2 2 2 2 ...
$ Datum : Factor w/ 1 level "06.08.2015": 1 1 1 1 1 1 1 1 1 1 ...
$ Gewicht : num 7.8 4.1 10.8 51.2 33.3 17.9 40.6 11.7 35.1 7.1 ...
$ IngestionRate: num 0.385 1.057 1.598 0.164 0.396 ...
我做了这样的子集:
K_NPK<-subset(Agro,Agro$Ressource!="N")
我的模特:
mod4 <- lmer(IngestionRate~Ressource+(1|Gewicht), data=K_NPK)
答案:
错误:每个分组因子的级别数必须<&lt;观察次数
但如果我做这个子集
N_K<-subset(Agro,Agro$Ressource!="NPK")
和这个模型
mod4 <- lmer(IngestionRate~Ressource+(1|Gewicht), data=N_K)
如果运行则没有错误。
我希望你明白我尝试做什么。 谁能告诉我什么错了?