В линейных моделях и особенно в ANOVA контраст представляет собой линейную комбинацию параметров с коэффициентами, суммируемыми до нуля. Он используется для проверки соответствующей нулевой гипотезы. Контрасты особенно часто используются с категориальными предикторами (факторами) для сравнения групп (категорий). [См. Также тег "категориальное кодирование"]
Пожалуйста, предоставьте R код, который позволяет проводить ANOVA между субъектами с -3, -1, 1, 3 контрастами. Я понимаю, что есть спор относительно соответствующего типа суммы квадратов (SS) для такого анализа. Однако тип SS по умолчанию, используемый в SAS и SPSS (тип III), считается стандартом в...
Я выполнил повторные измерения ANOVA в R следующим образом: aov_velocity = aov(Velocity ~ Material + Error(Subject/(Material)), data=scrd) summary(aov_velocity) Какой синтаксис в R можно использовать для выполнения пост-специального теста после ANOVA с повторными измерениями? Подойдет ли тест Тьюки...
Я не могу понять использование полиномиальных контрастов в регрессионном подборе. В частности, я имею в виду кодировку, используемую Rдля выражения интервальной переменной (порядковой переменной с одинаково разнесенными уровнями), описанную на этой странице . В примере на этой странице , если я...
Допустим, у меня есть простой факторный эксперимент 2x2, на котором я хочу провести ANOVA. Вот так, например: d <- data.frame(a=factor(sample(c('a1','a2'), 100, rep=T)), b=factor(sample(c('b1','b2'), 100, rep=T))); d$y <- as.numeric(d$a)*rnorm(100, mean=.75, sd=1) + as.numeric(d$b)*rnorm(100,...
Если вы подходите к линейной или смешанной модели, существуют различные типы кодировок, доступных для преобразования категориальной или номинальной вариабельной переменной в ряд переменных, для которых оцениваются параметры, такие как фиктивная кондукция (по умолчанию R) и кодирование эффектов. Я...
Я пытаюсь выполнить повторные измерения Anova в R с последующими некоторыми конкретными контрастами в этом наборе данных. Я думаю, что правильный подход будет использовать Anova()из пакета автомобилей. Позволяет проиллюстрировать мой вопрос на примере взят из ?Anovaиспользуя OBrienKaiserданные...
Я отправил это в mathoverflow, и никто не отвечает: Метод Шеффе для выявления статистически значимых контрастов широко известен. Контраст среди средств , из популяций является линейной комбинацией , в котором , и скалярное кратное контрастности - это, по сути, один и тот же контраст, поэтому можно...
Начинающий пользователь R здесь борется с повторными измерениями ANOVA. У меня есть набор данных, который состоит из одного фактора между субъектами с 4 уровнями (закодированного в единую переменную, называемую «группы»), и одного из фактора субъектов с 4 уровнями (кодированного в четырех отдельных...
Пожалуйста, дайте мне идею, как эффективно перекодировать категориальную переменную (фактор) в набор ортогональных полиномиальных контрастных переменных. Для многих типов переменных контраста (например, отклонения, простые, Хелмерта и т. Д.) Проход: Составьте матрицу коэффициентов контрастности,...
Я делаю односторонний ANOVA (для каждого вида) с индивидуальными контрастами. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 где я сравниваю интенсивность 0,5 против 5, 5 против 12,5 и так далее. Это данные, над которыми я работаю со следующими результатами...
Используя R, я установил линейную модель для одной переменной ответа из смеси непрерывных и дискретных предикторов. Это базовый уровень, но мне трудно понять, как работает коэффициент для дискретного фактора. Концепция: Очевидно, что коэффициент непрерывной переменной «х» применяется в форме, y =...
У меня есть модель регрессии, которая выглядит следующим образом:Y=β0+β1X1+β2X2+β3X3+β12X1X2+β13X1X3+β123X1X2X3Y=β0+β1X1+β2X2+β3X3+β12X1X2+β13X1X3+β123X1X2X3Y = \beta_0+\beta_1X_1 + \beta_2X_2 + \beta_3X_3 +\beta_{12}X_1X_2+\beta_{13}X_1X_3+\beta_{123}X_1X_2X_3 ... или в обозначении R: y ~ x1 + x2...
Я использую lmer в R, чтобы проверить влияние условия ( cond) на некоторый результат. Вот некоторые составные данные, где s - идентификатор субъекта a, bи c- условия. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30, -2, 1), b = rnorm(30, -3,...