Я получил набор данных и коды для оценки дисперсионной составляющей. В наборе данных «у» - ответ, «топливо» зафиксировано, а «день» имеет случайные эффекты в линейной смешанной модели. Я могу оценить случайный эффект из-за «дня», и нет никакой возможности иметь p-значение, соответствующее оценке случайного эффекта. Как я могу оценить значение p? Есть ли другой пакет для этого? Любой ответ будет высоко ценится.
day<-c(1, 1, 2, 2, 3, 3) fuel<-c(1,1,1,2,2,2) y<-c(32.32048, 32.28583, 31.97912, 33.61024, 33.33915, 33.63342) dataTB<-data.frame(day,fuel,y) library(nlme) contrastTB.lme<-lme(y~fuel, random= ~1 | day,data=dataTB) summary(contrastTB.lme)