Стандартные ошибки случайных моделей перехвата одинаковы

В R I есть модель случайного перехвата, использующая набор данных ветчины из lmerTest(). Существует около 80 потребителей, и модель обеспечит случайный перехват для всех 80 клиентов в ranef(m.ham), Когда я смотрю на стандартные ошибки этих случайных перехватов, они все одинаковы. См "se"объект. Почему стандартные ошибки будут одинаковыми для всех 80 клиентов. Глядя на несколько клиентов, использующих ham %>% filter(Consumer == 1) Вы можете видеть, что Informed.liking Значения для клиентов различны, поэтому разве это не должно приводить к разнице и стандартным ошибкам?

Смотрите строчку:

ham  %>% group_by(Consumer) %>%  summarise(COUNT = n(), consumerVariance = var(Informed.liking ) ) %>% as.data.frame()

который показывает дисперсию Informed.liking, отличается для каждой группы, и у всех групп есть 8 наблюдений.

library(lmerTest)
library(dplyr)
library(lme4)
ham  %>% group_by(Consumer) %>%  summarise(COUNT = n() ) %>% as.data.frame()

### look at a few cusomters
ham %>%  filter(Consumer == 1) 
ham %>%  filter(Consumer == 2)
ham %>%  filter(Consumer == 3)

m.ham <- lmer(Informed.liking ~ 1 + (1|Consumer), data = ham) #Product*Information*Gender + 
fixef(m.ham)
ranef(m.ham)
re = ranef(m.ham,condVar=TRUE)
qq = attr(ranef(m.ham,condVar=TRUE)[[1]],"postVar") #### VARIANCE OF RANDOM INTERCEPTS
cols <- 1:(dim(qq)[1])
se  <- unlist(lapply(cols, function(i) sqrt(qq[i, i, ]))) #standard errors OF RANDOM INTERCEPTS

0 ответов

Другие вопросы по тегам