Я работал с некоторыми данными, которые имеют некоторые проблемы с повторными измерениями. При этом я заметил очень различное поведение lme()и lmer()использование моих тестовых данных и хочу знать почему. Поддельный набор данных, который я создал, содержит измерения роста и веса для 10 предметов,...
Является ли преобразование модели многомерной линейной регрессии в множественную линейную регрессию полностью эквивалентным? Я не имею в виду , просто запустив TTt отдельных регрессий. Я читал это в нескольких местах (Байесовский анализ данных - Гельман и др. И Многовариантная старая школа -...
Я использую R и анализирую свои данные с помощью GLM с биноминальной ссылкой. Я хочу знать, что означает перехват в выходной таблице. Перехват для одной из моих моделей существенно отличается, однако переменная - нет. Что это значит? Что такое перехват? Я не знаю, просто ли я запутываю себя, но,...
Я понимаю, что разница между ними связана с тем, оценивается ли группирующая переменная в модели как фиксированный или случайный эффект, но мне не ясно, почему они не одинаковы (если они не совпадают). Меня особенно интересует, как это работает при использовании оценки небольшой площади, если это...
У нас много хороших дискуссий по поводу идеального разделения в логистической регрессии. Например, логистическая регрессия в R привела к идеальному разделению (феномен Хаука-Доннера). Что теперь? и модель логистической регрессии не сходится . Лично я все еще чувствую, что это не интуитивно понятно,...
У меня есть некоторые данные в [0,1], которые я хотел бы проанализировать с помощью бета-регрессии. Конечно, что-то нужно сделать, чтобы приспособить значения 0,1. Мне не нравится изменять данные, чтобы соответствовать модели. Кроме того, я не верю, что нулевая и 1 инфляция - это хорошая идея,...
Мне нужны некоторые советы относительно двух основных дилемм в моем исследовании, которое представляет собой исследование трех крупных фармацевтических препаратов и инноваций. Количество патентов в год является зависимой переменной. Мои вопросы Каковы наиболее важные критерии для хорошей модели?...
Это кажется таким элементарным, но я всегда застреваю на этом этапе ... Большинство данных, с которыми я имею дело, являются ненормальными, и большинство анализов основано на структуре GLM. Для моего текущего анализа у меня есть переменная ответа, которая является "скоростью ходьбы" (метры в...
У меня есть график остаточных значений линейной модели в зависимости от подогнанных значений, где гетероскедастичность очень ясна. Однако я не уверен, как мне поступить сейчас, потому что, насколько я понимаю, эта гетероскедастичность делает мою линейную модель недействительной. (Это правильно?)...
Допустим, у меня есть набор данных «почечный катетер». Я пытаюсь смоделировать кривую выживания, используя модель Кокса. Если я рассматриваю модель Кокса: мне нужна оценка базовой опасности. Используя встроенную функцию пакета R , я легко могу сделать это так:ч (t,Z) =h0exp(b'Z) ,h(T,Z)знак...
Рассмотрим следующие данные из двустороннего дизайна предметов: df <- "http://personality-project.org/r/datasets/R.appendix4.data" df <- read.table(df,header=T) head(df) Observation Subject Task Valence Recall 1 1 Jim Free Neg 8 2 2 Jim Free Neu 9 3 3 Jim Free Pos 5 4 4 Jim Cued Neg 7 5 5 Jim...
Я пытаюсь использовать lmeиз nlmeпакета, чтобы повторить результаты aovдля повторных мер ANOVA. Я сделал это для однофакторного эксперимента с повторными измерениями и для двухфакторного эксперимента с одним фактором между субъектами и одним фактором внутри субъекта, но у меня возникли проблемы с...
У меня есть таблица с результатами команд. Первая команда с 10 очками побеждает. В каждой команде 2 игрока. Игроки все время играют с разными товарищами по команде, хотя они выбраны не случайно. Индивидуальные оценки не сохраняются. Таким образом, у нас Билл и Боб обыграли Энди и Алису со счетом...
Скажем, мы работаем над моделью случайных эффектов некоторых данных подсчета во времени, и мы хотим контролировать некоторые тенденции. Обычно вы делаете что-то вроде: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") включить квадратную форму для t. Можно ли использовать более...
В основном меня интересует, как применяются различные ковариационные структуры и как рассчитываются значения внутри этих матриц. Такие функции, как lme (), позволяют нам выбирать, какую структуру мы бы хотели, но я бы хотел знать, как они оцениваются. Рассмотрим модель линейных смешанных эффектов...
Скажем, мы должны GLMMs mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Эти модели не являются вложенными в обычном смысле: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) поэтому мы не можем поступить...
Если я фиксирую значения наблюдаемых узлов MRF, становится ли он
Я пытаюсь создать модель, для которой у меня есть переменная ответа, которая составляет пропорцию между 0 и 1, это включает довольно много 0 и 1, но также и много значений между ними. Я думаю о попытке бета-регрессии. Пакет, который я нашел для R (betareg), допускает только значения в диапазоне от...
Мне интересно узнать о том, как оценивать и оценивать отдельных лиц в группе, которые взаимодействуют / соревнуются только попарно (например, такие системы, как рейтинговая система ELO для шахмат). Есть ли какие-либо методы перехода или более точные и продвинутые методы там? Существуют ли пакеты R,...
Предположим, у вас есть единое сечение данных, где отдельные лица находятся в группах (например, учащиеся в школах), и вы хотите оценить модель формы, в Y_i = a + B*X_iкоторой Xесть вектор характеристик индивидуального уровня и aконстанта. В этом случае предположим, что ненаблюдаемая неоднородность...