Я пытаюсь выполнить повторные измерения Anova в R с последующими некоторыми конкретными контрастами в этом наборе данных. Я думаю, что правильный подход будет использовать Anova()из пакета автомобилей. Позволяет проиллюстрировать мой вопрос на примере взят из ?Anovaиспользуя OBrienKaiserданные...
У меня очень большой набор данных и около 5% случайных значений отсутствуют. Эти переменные связаны друг с другом. В следующем примере набор данных R - просто игрушечный пример с фиктивными коррелированными данными. set.seed(123) # matrix of X variable xmat <- matrix(sample(-1:1, 2000000,...
Как правило, когда встречаются непрерывные, но искаженные показатели результата в продольном дизайне (скажем, с одним эффектом между субъектами), общий подход заключается в преобразовании результата в нормальность. Если ситуация экстремальная, например, с усеченными наблюдениями, можно подумать и...
Начинающий пользователь R здесь борется с повторными измерениями ANOVA. У меня есть набор данных, который состоит из одного фактора между субъектами с 4 уровнями (закодированного в единую переменную, называемую «группы»), и одного из фактора субъектов с 4 уровнями (кодированного в четырех отдельных...
mgcvПакет Rимеет две функции для установки взаимодействия Тензор продукта: te()и ti(). Я понимаю основное разделение труда между ними (подгонка нелинейного взаимодействия против разложения этого взаимодействия на основные эффекты и взаимодействие). Чего я не понимаю, так это почему te(x1, x2)и...
Мой вопрос основан на этом ответе, который показал, какая lme4::lmerмодель соответствует двусторонним повторным измерениям ANOVA: require(lme4) set.seed(1234) d <- data.frame( y = rnorm(96), subject = factor(rep(1:12, 4)), a = factor(rep(1:2, each=24)), b = factor(rep(rep(1:2, each=12))), c =...
Таким образом, предполагая, что есть смысл проверить допущение нормальности для ановы (см. 1 и 2 ) Как это можно проверить в R? Я ожидал бы сделать что-то вроде: ## From Venables and Ripley (2002) p.165. utils::data(npk, package="MASS") npk.aovE <- aov(yield ~ N*P*K + Error(block), npk)...
Я хотел бы оценить несколько различных моделей, которые обеспечивают предсказания поведения на ежемесячном уровне. Данные сбалансированы, и 100 000 и T = 12. В результате посещение концерта происходит в определенном месяце, поэтому оно равно нулю для ~ 80% людей в любом месяце, но есть длинный...
У меня проблемы с определением условий ошибки для двухсторонних повторяющихся измерений ANOVA в R. Мои данные состоят из оценок плотности древесины для трех радиальных положений (внутреннего, среднего и внешнего) вдоль керна, извлеченного из дерева. Всего существует 20 видов деревьев, 6 особей...
Я анализирую результаты эксперимента времени реакции в R. Я провел повторные измерения ANOVA (1 внутри-субъектный фактор с 2 уровнями и 1 между субъектный фактор с 2 уровнями). Я запустил аналогичную линейную смешанную модель, и я хотел обобщить результаты lmer в виде таблицы ANOVA с использованием...
Ради простого примера предположим, что есть две модели линейной регрессии Модель 1 имеет три предсказатели, x1a, x2b, иx2c Модель 2 имеет три предиктора из модели 1 и два дополнительных предиктора x2aиx2b Существует уравнение регрессии населения, где объясняется дисперсия населения для Модели 1 и...
Я работаю над набором данных. После использования некоторых методов идентификации моделей я разработал модель ARIMA (0,2,1). Я использовал detectIOфункцию в пакете TSAв R, чтобы обнаружить инновационный выброс (IO) на 48-м наблюдении за моим исходным набором данных. Как включить этот выброс в мою...
Мой вопрос очень тесно связан с предыдущим постом указания срока Error () в повторных измерениях ANOVA в R . Тем не менее, я хотел бы получить более полное представление о том, как определить термин ошибки. Предположим, у меня есть двухсторонний повторный ANOVA. Фактором межгруппового эффекта...
У меня есть набор данных, который содержит ~ 7500 анализов крови от ~ 2500 человек. Я пытаюсь выяснить, увеличивается ли вариабельность анализов крови с течением времени между двумя тестами. Например - я беру вашу кровь для базового теста, а затем сразу же беру второй образец. Через полгода я...
Я пытаюсь вписать модель с дискретным временем в R, но я не уверен, как это сделать. Я читал, что вы можете организовать зависимую переменную в разных строках, по одной для каждого временного наблюдения, и использовать glmфункцию со ссылкой logit или cloglog. В этом смысле, у меня есть три колонки:...
Это всего лишь пример, с которым я сталкивался несколько раз, поэтому у меня нет примеров данных. Запуск модели линейной регрессии в R: a.lm = lm(Y ~ x1 + x2) x1является непрерывной переменной x2является категориальным и имеет три значения, например, «Низкий», «Средний» и «Высокий». Однако вывод,...
Проще говоря, у меня есть модель множественной регрессии или ANOVA, но переменная отклика для каждого человека является криволинейной функцией времени. Как я могу сказать, какие из переменных правой стороны ответственны за существенные различия в формах или вертикальных смещениях кривых? Это...
Почему повторные измерения ANOVA предполагают сферичность? Под сферичностью я подразумеваю предположение, что дисперсия всех парных различий между группами должна быть одинаковой. В частности, я не понимаю, почему это должно быть предположение, а не то, что отклонения наблюдаемых групповых оценок...
Структура данных > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1 1 2 1 2...
Я выполнил трехсторонние повторные измерения ANOVA; какие специальные анализы действительны? Это полностью сбалансированный дизайн (2x2x2) с одним из факторов, имеющих повторное измерение внутри субъекта. Мне известны многовариантные подходы к повторным измерениям ANOVA в R, но мой первый инстинкт...