Ошибка в glht.matrix(model = new("glmerMod", resp = new("glmResp", .xData = <environment>),: 'ncol (linfct)' не равно 'length(coef(model))'
Я запустил GLMM:
**m7 = glmer(CF ~ (Morph * Habitat)/Can+ (1|Subject), data = Muntjak2, family = poisson(link = "log"), nAQ=0, control=glmerControl(check.nlev.gtr.1="ignore"))**
**summary(Muntjak2)**
*Subject Sampling Morph Species CF HC Habitat Canopy Can Station
Мин.:1 мин.: 1.00 ММ: 5 Месяц:16 Мин.: 1.00 I:47 SWLS:51 мин.: 0.00 H:55 MBJ1: 4
1-я четверть: 1-я четверть: 27,25 MP:11 NM:82 1-я четверть: 3,00 II:51 TERR:47 1-я четверть: 50.00 L:22 BS4: 3
Медиана: 1 Медиана: 53,50 NMM:51
Среднее значение: 10,50 Среднее значение: 70,00 M:21 C2: 3 Среднее значение: 1 Среднее значение: 56,90 NMP:31 Среднее значение: 23,35 Среднее значение: 63,06 CH1: 3 3-
й квартал: 1 3- й квартал: 86,75 3- й квартал: 28,00 3-й квартал: 90.00 GS-1: 3
Макс.: 1 Макс.: 115.00 Макс.:189.00 Макс.: 100.00 ЛК2: 3
(Другое):79
**str(Muntjak2)**
'data.frame': 98 obs. of 10 variables:
$ Subject : int 1 1 1 1 1 1 1 1 1 1 ...
$ Sampling: int 1 2 3 4 5 6 7 8 9 10 ...
$ Morph : Factor w/ 4 levels "MM","MP","NMM",..: 2 2 2 2 2 2 2 2 2 2 ...
$ Species : Factor w/ 2 levels "M","NM": 1 1 1 1 1 1 1 1 1 1 ...
$ CF : int 4 24 5 1 2 3 1 1 6 10 ...
$ HC : Factor w/ 2 levels "I","II": 1 1 2 1 2 2 2 2 1 1 ...
$ Habitat : Factor w/ 2 levels "SWLS","TERR": 2 2 1 2 1 1 1 1 2 2 ...
$ Canopy : int 20 50 90 90 90 70 0 50 40 40 ...
$ Can : Factor w/ 3 levels "H","L","M": 2 3 1 1 1 1 2 3 3 3 ...
$ Station : Factor w/ 57 levels "BR1","BS1","BS2",..: 8 9 25 28 37 51 5 38 15 27 ...*
После успешного запуска модели я хочу провести апостериорный тест, чтобы узнать, значительно ли различается ли скорость захвата различных морфов двух видов в двух типах среды обитания (SWLS) и (TERR).
Я побежал: summary(glht(m7, lsm(pairwise ~ Morph * Habitat)), test=adjusted(type = "holm")) Ошибка:
**Error in glht.matrix(model = new("glmerMod", resp = new("glmResp", .xData = <environment>), :
‘ncol(linfct)’ is not equal to ‘length(coef(model))’**
Я застрял на какое-то время. Любая помощь будет высоко ценится.