Почему я получаю эти сообщения об ошибках при вычислении post-hoc теста на glm - PullRequest
Новая
       79

Почему я получаю эти сообщения об ошибках при вычислении post-hoc теста на glm

0 голосов
/ 25 января 2019

Я хочу вычислить post-hoc тест на glm, но продолжаю сталкиваться с сообщениями об ошибках.Код, который я запускаю, выглядит следующим образом: 

  # Read data
  read.table("adult_pref.txt")

  # Assign a name to the data
  epi <- read.table("adult_pref.txt", 
                     ,col.names=c('generation','preference','A','B'),
                     header=TRUE, sep="\t", na.strings="NA", dec=".", strip.white=TRUE)

  glm <- glm(cbind(A,B)~generation, data=epi, family=quasibinomial(link="logit"))
  anova(glm, test="Chisq")
  library(multcomp)

  pref.test <- glht(glm, linfct=mcp(preference="Tukey"))
summary(pref.test)

Когда я запускаю pref.test <- glht(glm, linfct=mcp(preference="Tukey")), я вижу следующую ошибку: 

    Error in mcp2matrix(model, linfct = linfct) : 
  Variable(s) ‘preference’ have been specified in ‘linfct’ but cannot be found in ‘model’!

Кажется, что R не может найти переменнуюпредпочтение в модели.Просто чтобы дать вам представление о том, как выглядят мои данные: 

             V1           V2          V3          V4
1   generation   preference      A            B
2           F0        -0.43 28.47457627 71.52542373
3           F0        -0.82 9.150326797  90.8496732
4           F0           -1           0         100
5           F0        -0.38 30.83832335 69.16167665
6           F0           -1           0         100

Кто-нибудь знает, что я могу делать неправильно?

...