Ошибка в glht.matrix (model = new ("glmerMod", resp = new ("glmResp", .xData = <environment>),: 'ncol (linfct)' не равно 'length (coef (model)) " - PullRequest
Новая
       80

Ошибка в glht.matrix (model = new ("glmerMod", resp = new ("glmResp", .xData = <environment>),: 'ncol (linfct)' не равно 'length (coef (model)) "

0 голосов
/ 08 апреля 2020

Я управлял 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 1-й кв .: 27,25 MP: 11 NM: 82 1-й кв .: 3,00 1005 * Медиана: 1 Медиана: 53,50 НММ: 51 Медиана: 10,50 Медиана: 70,00 М: 21 С2: 3
Среднее: 1 Среднее: 56,90 NMP: 31 Среднее: 23,35 Среднее: 63,06 СН1: 3
3-й кв. : 1 3-й кв .: 86,75 3-й кв .: 28,00 3-й кв .: 90,00 GS-1: 3
Макс. : 1 Макс. : 115,00 Макс. : 189,00 Макс. : 100,00 LK2: 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 ...*

После успешного запуска модели я хочу сделать сообщение ho c тест, чтобы узнать, значительно ли различаются скорости захвата различных морфов двух видов в двух типах местообитаний (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))’**

Я застрял на некоторое время. Любая помощь будет принята с благодарностью.

...