Вопросы с тегом «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), data=d) для …


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), так как он, кажется, приспосабливает эту настройку напрямую. Мои данные …

1
Остаточная диагностика в регрессионных моделях на основе MCMC
Недавно я приступил к подгонке регрессионно-смешанных моделей в байесовской структуре, используя алгоритм 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.