Запитання з тегом «lme4-nlme»

lme4 і nlme - пакети R, які використовуються для встановлення лінійних, узагальнених лінійних та нелінійних моделей змішаних ефектів. Для загальних питань щодо змішаних моделей використовуйте тег [змішана модель].

1
Як інтерпретувати коефіцієнти багатоваріантної змішаної моделі в lme4 без загального перехоплення?
Я намагаюся встановити багатоваріантну (тобто множину реакцій) змішану модель R. Окрім пакунків ASReml-rта SabreRпакунків (для яких потрібне зовнішнє програмне забезпечення), здається, це можливо лише в MCMCglmm. У статті, що додається до MCMCglmmпакету (с. 6), Джаррод Хадфілд описує процес встановлення такої моделі, як переформатування декількох змінних відповідей в одну змінну довгого …

2
Ієрархічні моделі для численних порівнянь - контекст множинного результату
Я щойно читав (пере) читав "Гельман", чому нам (як правило) не потрібно турбуватися про багаторазове порівняння . Зокрема, у розділі "Кілька результатів та інші виклики" згадується використання ієрархічної моделі для ситуацій, коли в одній особі / підрозділі є різні пов'язані заходи в різні періоди / умови. Схоже, він має ряд …

1
Запитання щодо визначення лінійних змішаних моделей у R для даних повторних вимірювань з додатковою структурою вкладення
Структура даних > 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 …

1
Чому я отримую різко різні результати для poly (raw = T) проти poly ()?
Я хочу моделювати дві різні змінні часу, деякі з яких сильно колінеарні в моїх даних (вік + когорта = період). Роблячи це, у мене виникли проблеми з lmerвзаємодією та взаємодією poly(), але, мабуть, це не обмежується lmer, я отримав однакові результати з nlmeIIRC. Очевидно, що мого розуміння того, що робить …

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

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

1
Що повідомляє lsmeans для узагальненої лінійної моделі, такої як змішана модель Пуассона (підходить до glmer)?
Я аналізую дані відстеження очей із розробленого експерименту. Спрощена версія моїх даних виглядає приблизно так (дані dput () можна отримати тут ), head(lookDATA) participant fixationImage fixationCount 1 9 Automobile 81 2 9 Bird 63 3 9 Chair 82 4 9 Dog 64 5 9 Face 90 6 9 Plant 75 …

2
P значення для терміну взаємодії в моделях зі змішаними ефектами з використанням lme4
Я аналізую деякі поведінкові дані, використовуючи, lme4в Rосновному дотримуючись чудових навчальних посібників Bodo Winter , але я не розумію, чи правильно обробляю взаємодію. Гірше, що ніхто більше не бере участь у цьому дослідженні, не використовує змішаних моделей, тому я трохи похитнувся, коли потрібно переконатися, що все правильно. Замість того, щоб …

2
Варіаційно-коваріаційна структура для випадкових ефектів у lme4
Яка структура дисперсії-коваріації за замовчуванням для випадкових ефектів у пакеті glmerчи lmerв ньому lme4? Як можна вказати іншу структуру коваріації дисперсії для випадкових ефектів у коді? Я не зміг знайти жодної інформації щодо цього в lme4документації.

2
REML проти ML stepAIC
Я відчуваю себе переповненим після спроби зануритися в літературу про те, як запустити мішаний аналіз моделі після його використання з використанням AIC для вибору найкращої моделі чи моделі. Я не думаю, що мої дані є настільки складними, але я шукаю підтвердження того, що те, що я зробив, є правильним, а …

3
Бактерії, зібрані на пальцях після декількох поверхневих контактів: ненормальні дані, повторні заходи, схрещені учасники
Вступ У мене є учасники, які неодноразово торкаються забруднених поверхонь кишковою паличкою за двох умов ( А = носіння рукавичок, В = рукавички відсутні). Хочу знати, чи є різниця між кількістю бактерій на кінчиках пальців із рукавичками та без них, а також між кількістю контактів. Обидва фактори є всередині учасника. …

3
Що робити зі співвідношенням випадкових ефектів, що дорівнює 1 або -1?
Не настільки рідкісне явище при роботі зі складними максимально змішаними моделями (оцінка всіх можливих випадкових ефектів для даних та моделі) є ідеальним (+1 або -1) або майже ідеальним співвідношенням між деякими випадковими ефектами. З метою обговорення, дотримуйтесь наступної моделі та резюме моделі Model: Y ~ X*Cond + (X*Cond|subj) # Y …

1
Допоможіть інтерпретувати дані підрахунку GLMM за допомогою lme4 glmer та glmer.nb - Від'ємний біном на проти Пуассона
У мене є деякі питання щодо специфікації та інтерпретації GLMM. 3 питання, безумовно, статистичні, а 2 - конкретніше про Р. Я публікую тут, оскільки, зрештою, я думаю, що це інтерпретація результатів ГЛМ. На даний момент я намагаюся вписатись в GLMM. Я використовую дані перепису США з Бази даних про поздовжні …

1
Модель змішаних ефектів із сплайнами
Мені підходить модель змішаних ефектів із терміном сплайна в додатку, де тенденція з часом є криволінійною. Однак те, що я хотів би оцінити, - чи виникає крива-лінійна тенденція через індивідуальне відхилення від лінійності, чи це вплив на рівні групи, який робить груповий рівень придатним, здається криволінійним. Я наводжу відтворювальний приклад, …
9 r  splines  lme4-nlme 

1
Як встановити власні контрасти з lmer в R
Я використовую 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, 1), c = rnorm(30, …

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.