Запитання з тегом «mixed-model»

Припустимо, у мене є деяка змінна відповіді яка вимірювалася від ї побратимів у й сім'ї. Крім того, деякі дані поведінки були зібрані одночасно від кожного суб'єкта. Я намагаюся проаналізувати ситуацію за допомогою наступної лінійної моделі змішаних ефектів:yijyijy_{ij}jjjiiixijxijx_{ij} yij=α0+α1xij+δ1ixij+εijyij=α0+α1xij+δ1ixij+εijy_{ij} = \alpha_0 + \alpha_1 x_{ij} + \delta_{1i} x_{ij} + \varepsilon_{ij} де і …

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

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

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

Чи може хто-небудь сказати мені різницю між використанням aov()та lme()аналізом поздовжніх даних та як інтерпретувати результати цих двох методів? Нижче я аналізую той же набір даних , використовуючи aov()і lme()та отримав 2 різні результати. З aov(), я отримав вагомий результат за час взаємодії з лікуванням, але, якщо підходити до лінійної …

14 r  mixed-model  repeated-measures  panel-data 

Дуглас Бейтс зазначає, що наступні моделі є еквівалентом "якщо матриця коваріації дисперсії для випадкових ефектів, що оцінюються за вектором, має спеціальну форму, що називається складовою симетрією" ( слайд 91 у цій презентації ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor + (1|group) …

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 y_11 1 2 y_12 2 1 y_21 2 2 y_22 2 3 y_23 ... ... ... Я хочу знати співвідношення між вимірюваннями, проведеними на …

13 mixed-model  intraclass-correlation 

Я часто чув, що моделі LME є більш надійними при аналізі даних про точність (тобто, в експериментах з психології), оскільки вони можуть працювати з біноміальними та іншими ненормальними розподілами, які традиційні підходи (наприклад, ANOVA) не можуть. Що є математичною основою моделей LME, що дозволяють їм включати ці інші дистрибутиви, і …

13 mixed-model 

Я хотів би витягнути нахили для кожної людини в моделі змішаного ефекту, як це викладено в наступному параграфі Моделі змішаних ефектів використовувались для характеристики окремих шляхів зміни когнітивних підсумкових заходів, включаючи терміни віку, статі та років навчання як фіксованих ефектів (Laird and Ware, 1982; Wilson et al., 2000, 2002c) ... …

13 r  mixed-model 

Я читав про обчислення значень R2R2R^2 у змішаних моделях і після прочитання відповідей на R-sig, інших публікацій на цьому форумі (я б зв’язав декілька, але мені не вистачає репутації) та кілька інших посилань, я розумію, що використовуючи Значення R2R2R^2 у контексті змішаних моделей є складним. Однак я нещодавно натрапив на …

13 r  mixed-model  r-squared  lme4-nlme 

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

13 r  data-visualization  mixed-model 

Я хотів би отримати 95% довірчі інтервали для прогнозів нелінійної змішаної nlmeмоделі. Оскільки для цього не передбачено нічого стандартного nlme, мені було цікаво, чи правильно використовувати метод "інтервалів прогнозування населення", як це викладено в книзі книги Бен Болкер в контексті моделей, що відповідають максимальній вірогідності , заснованих на ідеї перекомпонування …

12 r  mixed-model  confidence-interval  lme4-nlme 

При багаторівневому моделюванні внутрішньокласова кореляція часто розраховується на основі ANOVA випадкових ефектів уi j= γ00+ уj+ еi jyij=γ00+uj+eij y_{ij} = \gamma_{00} + u_j + e_{ij} де - рівня 2, а - залишки рівня 1. Тоді отримуємо оцінки, та для дисперсії та відповідно, і підключаємо їх до наступного рівняння:е я J …

12 correlation  mixed-model  intraclass-correlation 

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

12 r  mixed-model  optimization  lme4-nlme  numerics