GLHT дает только один линейный тест - PullRequest
Новая
       5

GLHT дает только один линейный тест

0 голосов
/ 10 декабря 2018

Я перезапускаю код, который у меня был для теста glht, и я получаю меньше вывода, чем раньше.Я использую модель с двумя терминами и взаимодействием.Мои данные выглядят так: 

type   pop chemoheterotrophy
1 larvae Bibel       0.929379395
2 larvae Bibel       0.945769649
3 larvae Bibel       0.991292099
4 larvae Bibel       0.947357992
5 larvae Bibel       0.926360018
6 larvae Bibel       0.960807347

Есть два уровня типа и 5 уровней поп-музыки.Моя модель выглядит следующим образом: 

k <- matrix(c(28.53,17.72,31.07,7.95, 7.65, 28.53,17.72,31.07,7.95, 7.65),1)

test = glht(glm(chemoheterotrophy ~ type + pop + type*pop, data=boxdata, family=quasibinomial(link="logit")),linfct=k)

Когда я запускал это с summary(), я получал бы p-значения для intercept, pop, type и pop*type.Теперь я просто получаю: 

Simultaneous Tests for General Linear Hypotheses

Fit: glm(formula = chemoheterotrophy ~ type + pop + type * pop, family = quasibinomial(link = "logit"), 
    data = boxdata)

Linear Hypotheses:
       Estimate Std. Error z value Pr(>|z|)
1 == 0   -5.107     18.154  -0.281    0.778
(Adjusted p values reported -- single-step method)

Есть идеи, что происходит?Код и данные не изменились, но теперь я получаю другой вывод.

...