Проблема в журнале трансформированных контрастов - PullRequest
0 голосов
/ 01 апреля 2020

Я пытаюсь сделать ряд контрастов, чтобы сравнить поведение яйцекладки комаров в различных типах бассейнов. Я настроил свою модель и запустил свои данные следующим образом (в качестве примера предоставил только один контраст).

DMC<-read_excel('DMC.master.xlsx',sheet = 1)
DMC$Pool.Local.Comp <- as.factor(DMC$Pool.Local.Comp)
DMC$logSum <- log(DMC$Sum)

DMC$logSum[!is.finite(DMC$logSum)] <- 0
DMC.log.lmer <- lmer(logSum ~ Pool.Local.Comp + (1|Block), data = DMC)

levels(DMC$Pool.Local.Comp)
contest(DMC.log.lmer, c(-1,1,-1,1,0,0,0,0))

Проблема в том, что при запуске этого процесса я получаю p <0,001. Когда мой PI запускает это в SAS, он получает незначительное значение (не сказал мне, что именно). Кроме того, когда я рисую данные, это подтверждает его результат. Кто-нибудь знает, где я могу пойти не так? </p>

...