Эффект повторения в парных множественных сравнениях в модели GLM - PullRequest
Новая
       39

Эффект повторения в парных множественных сравнениях в модели GLM

0 голосов
/ 26 февраля 2019

Я хотел бы знать разницу и интерпретацию парных множественных сравнений в модели GLM.

Если мне нужно учесть повторения в модели.

См. Мой пример ниже:

Сначала я создаю некоторые переменные 

set.seed(123)
d<-NULL
N<-54
d$Production <- rgamma(N,10)
d$Feature <- ifelse(d$Production >7 & d$Production<10, c("A"),ifelse(d$Production>11, 
c("B"), c("C")))
d$Feature<-sort(d$Feature)  
#d$Temp<-rnorm(N,20,5)
d$Temp<-rep(1:3,18)
d$Rep<-sort(rep(1:3,18),decreasing=TRUE)
d<-as.data.frame(d)
#

Во-вторых, я делаю парные множественные сравнения без повторений: 

# Gamma glm model and pair-wise comparison without repetition 
m1<- glm(Production ~ Feature + Temp, family= Gamma, data = d)
cld(glht(m1, linfct = mcp(Feature = "Tukey")),decreasing=TRUE) 
#
  A   B   C 
"a" "a" "a"

В-третьих, я рассматриваю повторы: 

# Gamma glm model and pair-wise comparison with repetition 
m2<- glm(Production ~ Feature + Temp + Rep, family= Gamma, data = d)
cld(glht(m2, linfct = mcp(Feature = "Tukey")),decreasing=TRUE)
# 
  A   B   C 
"a" "ab" "b"

Мой вопрос: каков правильный статистический подход, если это решение меняет результаты?

Заранее спасибо

...