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

Змішані (також багаторівневі або ієрархічні) моделі - це лінійні моделі, що включають як фіксовані ефекти, так і випадкові ефекти. Вони використовуються для моделювання поздовжніх або вкладених даних.

2
Чому lme та aov повертають різні результати для повторних заходів ANOVA в R?
Я намагаюся перейти від використання ezпакету до lmeповторних заходів ANOVA (оскільки, я сподіваюся, я зможу використовувати власні контрасти на с lme). Дотримуючись поради з цього допису в блозі, я зміг налаштувати ту саму модель, використовуючи як aov(так і ezза запитом) та lme. Однак, хоча у прикладі, наведеному в цій публікації, …

2
Чи має сенс фіксований ефект вкладатись у випадковий, або як кодувати повторні заходи в R (aov та lmer)?
Я переглядав цей огляд формул lm / lmer R від @conjugateprior і збентежився наступним записом: Тепер припустимо, що A є випадковим, але B є фіксованим, а B вкладений у межах A. aov(Y ~ B + Error(A/B), data=d) Нижче подано аналогічну змішану формулу моделі lmer(Y ~ B + (1 | A:B), …


2
Як слід порівняти та затвердити моделі змішаних ефектів?
Яким чином (лінійні) моделі змішаних ефектів зазвичай порівнюються одна з одною? Я знаю, що можна використовувати тести на коефіцієнт ймовірності, але це не працює, якщо одна модель не є "підмножиною" іншої правильної? Чи завжди оцінка моделей df прямолінійна? Кількість фіксованих ефектів + ​​кількість дисперсійних компонентів, що оцінюються? Чи ігноруємо ми …

2
Чому я отримую нульову дисперсію випадкового ефекту в моїй змішаній моделі, незважаючи на певні зміни в даних?
Ми виконали логістичну регресію зі змішаними ефектами, використовуючи наступний синтаксис; # fit model fm0 <- glmer(GoalEncoding ~ 1 + Group + (1|Subject) + (1|Item), exp0, family = binomial(link="logit")) # model output summary(fm0) Тема та предмет - випадкові ефекти. Ми отримуємо непарний результат, який є коефіцієнтом і стандартним відхиленням для предмета, …

1
Показано, що 100 вимірювань для 5 суб'єктів дають набагато менше інформації, ніж 5 вимірювань для 100 предметів
На конференції я почув таке твердження: 100 вимірювань для 5 суб'єктів дають набагато менше інформації, ніж 5 вимірювань для 100 предметів. Це начебто очевидно, що це правда, але мені було цікаво, як можна це довести математично ... Я думаю, що можна використовувати лінійну змішану модель. Однак я не знаю багато …

2
Як застосувати біноміальний GLMM (glmer) до відсотків, а не так-ні?
У мене є експеримент повторного вимірювання, де залежна змінна є відсотком, і я маю кілька факторів як незалежних змінних. Я хотів би використати glmerз пакету R lme4трактувати це як проблему логістичної регресії (шляхом уточнення family=binomial), оскільки, здається, це вміщення безпосередньо відповідає. Мої дані виглядають так: > head(data.xvsy) foldnum featureset noisered …

1
Залишкова діагностика в регресійних моделях на основі МСМС
Нещодавно я взявся за пристосування регресійних змішаних моделей у байєсівській основі, використовуючи алгоритм MCMC (фактично функція MCMCglmm в R). Я вважаю, що я зрозумів, як діагностувати конвергенцію процесу оцінки (слід, графік гевеке, автокореляція, задній розподіл ...). Одне з того, що мене вражає в байєсівських рамках, - це те, що, здається, …

3
Як інтерпретувати основні ефекти, коли ефект взаємодії не суттєвий?
Я запустив узагальнену лінійну змішану модель в R і включив ефект взаємодії між двома прогнозами. Взаємодія була не суттєвою, але основні ефекти (обидва прогнози) були. Зараз багато прикладів підручників говорять про те, що якщо є значний вплив взаємодії, то основні ефекти не можна інтерпретувати. Але що робити, якщо ваша взаємодія …

4
Як спроектувати новий вектор на простір PCA?
Після проведення аналізу основних компонентів (PCA) я хочу спроектувати новий вектор на простір PCA (тобто знайти його координати в системі координат PCA). Я розрахував PCA мовою R за допомогою prcomp. Тепер я повинен мати можливість помножити свій вектор на матрицю обертання PCA. Чи повинні головні компоненти в цій матриці розташовуватися …
21 r  pca  r  variance  heteroscedasticity  misspecification  distributions  time-series  data-visualization  modeling  histogram  kolmogorov-smirnov  negative-binomial  likelihood-ratio  econometrics  panel-data  categorical-data  scales  survey  distributions  pdf  histogram  correlation  algorithms  r  gpu  parallel-computing  approximation  mean  median  references  sample-size  normality-assumption  central-limit-theorem  rule-of-thumb  confidence-interval  estimation  mixed-model  psychometrics  random-effects-model  hypothesis-testing  sample-size  dataset  large-data  regression  standard-deviation  variance  approximation  hypothesis-testing  variance  central-limit-theorem  kernel-trick  kernel-smoothing  error  sampling  hypothesis-testing  normality-assumption  philosophical  confidence-interval  modeling  model-selection  experiment-design  hypothesis-testing  statistical-significance  power  asymptotics  information-retrieval  anova  multiple-comparisons  ancova  classification  clustering  factor-analysis  psychometrics  r  sampling  expectation-maximization  markov-process  r  data-visualization  correlation  regression  statistical-significance  degrees-of-freedom  experiment-design  r  regression  curve-fitting  change-point  loess  machine-learning  classification  self-study  monte-carlo  markov-process  references  mathematical-statistics  data-visualization  python  cart  boosting  regression  classification  robust  cart  survey  binomial  psychometrics  likert  psychology  asymptotics  multinomial 

2
Попередження "Модель не вдалося збільшити" в lmer ()
За допомогою наступного набору даних я хотів побачити, чи змінюється реакція (ефект) щодо сайтів, сезону, тривалості та їх взаємодії. Деякі інтернет-форуми зі статистикою запропонували мені продовжувати використовувати лінійні моделі змішаних ефектів, але проблема полягає в тому, що оскільки репліки рандомізовані в межах кожної станції, у мене є мало шансів зібрати …

1
lme () та lmer (), що дають суперечливі результати
Я працював з деякими даними, які мають певні проблеми з повторними вимірюваннями. Роблячи це, я помітив дуже різну поведінку між lme()і lmer()використовуючи свої тестові дані, і хочу знати, чому. Створений мною підроблений набір даних містить вимірювання висоти та ваги для 10 предметів, зроблених двічі кожен. Я встановив дані, щоб між …

1
Чому орієнтовні значення найкращого лінійного неупередженого прогноктора (BLUP) відрізняються від найкращого лінійного неупередженого оцінювача (BLUE)?
Я розумію, що різниця між ними пов’язана з тим, чи змінна групування в моделі оцінюється як фіксований або випадковий ефект, але мені незрозуміло, чому вони не однакові (якщо вони не однакові). Мене конкретно цікавить, як це працює при використанні оцінки невеликих площ, якщо це актуально, але я підозрюю, що питання …

4
Як працює розподіл Пуассона при моделюванні постійних даних і чи призводить це до втрати інформації?
Співробітниця аналізує деякі біологічні дані для своєї дисертації з деякою неприємною гетероседастичністю (малюнок нижче). Вона аналізує це за змішаною моделлю, але все ще має проблеми з залишками. Перетворення змінних відповідей змін журналу очищує речі, і на основі відгуків на це питання, здається, є відповідним підходом. Однак спочатку ми думали, що …

2
Як я можу об'єднати задні засоби та достовірні інтервали після багаторазової імпутації?
Я використовував багаторазову імпутацію, щоб отримати ряд завершених наборів даних. Я використовував байєсівські методи для кожного із завершених наборів даних для отримання заднього розподілу параметра (випадковий ефект). Як можна об'єднати / об'єднати результати для цього параметра? Більше контексту: Моя модель є ієрархічною у розумінні окремих учнів (одне спостереження на одного …

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