Средние наименьших квадратов с их стандартными ошибками для объекта aov можно получить с помощью функции model.tables : npk.aov <- aov(yield ~ block + N*P*K, npk) model.tables(npk.aov, "means", se = TRUE) Интересно, как получить...

У меня проблемы с пониманием поведения метода наименьших квадратов. Ниже приведен игрушечный пример, в котором используется случайный набор данных для демонстрации моей проблемы. Сценарий таков: есть 10 станций, отобранных для показателя,...

У меня есть glm, который я хотел бы настроить для использования lsmeans. Следующий код создает модель (и, кажется, делает это правильно): library(lmerTest) data$group <- as.factor(data$grp) data$site <- as.factor(data$site) data$stimulus...

Я пытаюсь использовать lsmeans и его контраст для F-теста взаимодействия. По сути, я хотел бы воспроизвести то, что делает Stata с помощью команды contrast . Я хотел бы сделать это по двум причинам: в регрессионной модели, которая...

Я определил обобщенную линейную модель следующим образом: glm(formula = ParticleCount ~ ParticlePresent + AlgaePresent + ParticleTypeSize + ParticlePresent:ParticleTypeSize + AlgaePresent:ParticleTypeSize, family = poisson(link = "log"), data =...

Я хотел бы получить прогноз вероятности события из модели Коксма для заданных значений моих независимых переменных и на заданном временном шаге. Мои данные структурированы так: # Generate data set.seed(123) mydata <- data.frame(Site =...

Чтобы создать фрейм данных: num <- sample(1:25, 20) x <- data.frame("Day_eclosion" = num, "Developmental" = c("AP", "MA", "JU", "L"), "Replicate" = 1:5) model <- glmer(Day_eclosion ~ Developmental + (1 | Replicate), family =...