Вопросы с тегом «mixed-model»

Предположим, у меня есть некоторая переменная ответа которая была измерена от го брата в м семействе. Кроме того, некоторые поведенческие данные были собраны в то же время от каждого субъекта. Я пытаюсь проанализировать ситуацию с помощью следующей линейной модели смешанных эффектов: j i x i jyijyijy_{ij}jjjiiixijxijx_{ij} yij=α0+α1xij+δ1ixij+εijyij=α0+α1xij+δ1ixij+εijy_{ij} = \alpha_0 + …

14 mixed-model  lme4-nlme  covariance-matrix  non-independent 

Может кто-нибудь сказать мне, как заставить R оценить точку разрыва в кусочно-линейной модели (как фиксированный или случайный параметр), когда мне также нужно оценить другие случайные эффекты? Ниже я привел пример с игрушкой, который соответствует регрессии хоккейной клюшки / сломанной палки со случайными отклонениями наклона и случайной дисперсией y-точки пересечения для …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

Дуглас Бейтс утверждает, что следующие модели эквивалентны «если матрица дисперсии-ковариации для векторно-значных случайных эффектов имеет особую форму, называемую составной симметрией» ( слайд 91 в этой презентации ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor + (1|group) + (1|group:factor), data) В частности, Бейтс …

13 r  anova  mixed-model  repeated-measures  lme4-nlme 

Я вижу упомянутое в разных местах, что ANOVA делает оценку, используя метод моментов. Меня смущает это утверждение, потому что, хотя я не знаком с методом моментов, я понимаю, что это нечто отличное от метода максимальной вероятности и не эквивалентное ему; с другой стороны, ANOVA можно рассматривать как линейную регрессию с …

13 anova  mixed-model  maximum-likelihood  method-of-moments 

Из элементарной статистики я узнал, что с общей линейной моделью, чтобы выводы были достоверными, наблюдения должны быть независимыми. Когда происходит кластеризация, независимость может больше не сохраняться, приводя к неверному выводу, если это не учитывается. Одним из способов учета такой кластеризации является использование смешанных моделей. Я хотел бы найти примерный набор …

13 r  mixed-model  inference  independence 

Давайте рассмотрим этот гипотетический набор данных: set.seed(12345) num.subjects <- 10 dose <- rep(c(1,10,50,100), num.subjects) subject <- rep(1:num.subjects, each=4) group <- rep(1:2, each=num.subjects/2*4) response <- dose*dose/10 * group + rnorm(length(dose), 50, 30) df <- data.frame(dose=dose, response=response, subject=subject, group=group) мы можем использовать lmeдля моделирования ответа с моделью случайного эффекта: require(nlme) model <- …

13 mixed-model  linear-model  prediction 

Предположим, я измерил некоторую переменную у братьев и сестер, которые вложены в семьи. Структура данных выглядит следующим образом: семейный брат ------ ------- ----- 1 1 год_11 1 2 год_12 2 1 год_21 2 2 y_22 2 3 y_23 ... ... ... Я хочу знать корреляцию между измерениями, сделанными на братьях …

13 mixed-model  intraclass-correlation 

Может кто-нибудь сказать мне разницу между использованием aov()и lme()для анализа продольных данных и как интерпретировать результаты этих двух методов? Ниже я анализирую один и тот же набор данных aov()и получаю lme()2 разных результата. Таким образом aov(), я получил значительный результат во времени при взаимодействии с лечением, но при использовании линейной …

13 r  mixed-model  repeated-measures  panel-data 

Я часто слышал, что модели LME более надежны при анализе данных о точности (т. Е. В психологических экспериментах), поскольку они могут работать с биномиальными и другими ненормальными распределениями, которые не могут традиционные методы (например, ANOVA). Какова математическая основа моделей LME, которая позволяет им включать эти другие распределения, и что некоторые …

13 mixed-model 

Я хотел бы извлечь наклоны для каждого человека в модели смешанного эффекта, как обрисовано в общих чертах в следующем параграфе Модели смешанных эффектов использовались для характеристики отдельных путей изменения показателей когнитивной сводки, включая термины для возраста, пола и года обучения в качестве фиксированных эффектов (Laird and Ware, 1982; Wilson et …

13 r  mixed-model 

Я читал о расчете значений в смешанных моделях и после прочтения FAQ по R-sig, других постов на этом форуме (я бы связал несколько, но мне не хватает репутации) и нескольких других ссылок, которые я понимаю, используя Значения в контексте смешанных моделей сложны.R 2р2R2R^2р2R2R^2 Однако недавно я наткнулся на эти две …

13 r  mixed-model  r-squared  lme4-nlme 

Я анализировал некоторые данные с использованием линейного моделирования смешанных эффектов в R. Я планирую сделать плакат с результатами, и мне было просто интересно, может ли кто-нибудь, имеющий опыт работы с моделями смешанных эффектов, предложить, какие графики использовать для иллюстрации результатов модель. Я думал об остаточных графиках, графике подгоночных значений против …

13 r  data-visualization  mixed-model 

Я хотел бы получить 95% доверительные интервалы на предсказаниях нелинейной смешанной nlmeмодели. Поскольку для этого не предусмотрено ничего стандартного nlme, мне хотелось бы знать, правильно ли использовать метод «интервалов прогнозирования населения», как описано в главе книги Бена Болкера в контексте моделей, подходящих с максимальной вероятностью , основанных на идее пересобрать …

12 r  mixed-model  confidence-interval  lme4-nlme 

В многоуровневом моделировании внутриклассовая корреляция часто рассчитывается по ANOVA со случайными эффектами yij=γ00+uj+eijyij=γ00+uj+eij y_{ij} = \gamma_{00} + u_j + e_{ij} где - остатки уровня 2, а - остатки уровня 1. Затем мы получаем оценки и для дисперсии и соответственно и включаем их в следующее уравнение:е я J σ 2 U …

12 correlation  mixed-model  intraclass-correlation 

TL; DR: lme4оптимизация кажется линейной по количеству параметров модели по умолчанию и намного медленнее, чем эквивалентная glmмодель с фиктивными переменными для групп. Что я могу сделать, чтобы ускорить это? Я пытаюсь соответствовать довольно большой иерархической модели логита (~ 50 тыс. Строк, 100 столбцов, 50 групп). Подгонка нормальной модели логита к …

12 r  mixed-model  optimization  lme4-nlme  numerics