Я выполняю простую общую линейную смешанную модель, используя lmer.Я хотел бы узнать, есть ли какие-либо различия в обилии цветочных посетителей между участками обработки питательными веществами.У меня есть 7 питательных участков, каждый из которых повторяется один раз в 3 блока.Обработки: control, N, P, K, NPK, f_control, f_NPK.Блоки: A, B, C.
Мои данные выглядят так:
totvisobs<-plots$tot_vis_obs_withesc
Treatment<-plots$Treatment
Block<-plots$Block
Block<-as.factor(Block)
С:
Block: Factor w/ 3 levels "A","B","C": 1 1 1 1 1 1 1 2 2 2 ...
Treatment: Factor w/ 7 levels "control","f_control",..: 1 4 7 5 6 2 3 1 4 7 ...
totvisobs: int [1:21] 63 57 38 150 97 38 83 53 69 47 ...
Моя модель:
totvisobs_glmm<-lmer(totvisobs~Treatment+(1|blocks), REML = TRUE, data = plots)
Anova(totvisobs_glmm, test.statistic = "F")
Вот что я получаю:
Analysis of Deviance Table (Type II Wald F tests with Kenward-Roger df)
Response: totbeeobs
F Df Df.res Pr(>F)
Treatment 5.102 6 12 0.008054 **
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Предупреждающее сообщение: В бета - бетаH: большая длина объекта не кратна короткой длине объекта
Я не знаю, должен ли я игнорировать это предупреждение или нет.Я пытался все, чтобы решить эту проблему, но я не знаю, что делать, так как мои данные выглядят хорошо.Кто-нибудь знает, как с этим бороться?
Заранее спасибо!
Лизанна