Вопросы по теме 'mixed-models'

Получение обобщенных средств наименьших квадратов для фиксированных эффектов в nlme или lme4
Средние наименьших квадратов с их стандартными ошибками для объекта aov можно получить с помощью функции model.tables : npk.aov <- aov(yield ~ block + N*P*K, npk) model.tables(npk.aov, "means", se = TRUE) Интересно, как получить...
3651 просмотров
schedule 09.06.2024

Как подавить корреляционную таблицу в LME?
В стандартном примере функции lme() в пакете nlme R: fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1) summary(fm2) появляется корреляционная таблица: Correlation: (Intr) age age -0.813...
1843 просмотров
schedule 25.06.2022

Перезапустить оценку модели смешанного эффекта с ранее оцененными значениями
Я использую lmer() в пакете lme4 для оценки моделей со смешанными эффектами. Это работает хорошо, но теперь я хочу запустить процесс оценки для фиксированного количества итераций, а затем возобновить процесс, указав начальные значения,...
2894 просмотров
schedule 14.02.2022

Понимание кода модели смешанных эффектов с использованием Stata
От исследователя, использующего Stata, я получил данные и следующий код: . xtmelogit choicea network2 treatment2 treatment3 trustor2 if (trustor>0) & (period>0) || subgrid: || subjectidtr: , var Правильно ли, что ||...
294 просмотров
schedule 13.02.2024

R: mix() в пакете mixdist возвращает ошибку
Я установил пакет mixdist в R для объединения дистрибутивов. В частности, я использую функцию mix() . См. документацию. По сути, я получаю Error in nlm(mixlike, lmixdat = mixdat, lmixpar = fitpar, ldist = dist, : missing value in...
1264 просмотров
schedule 31.12.2023

Модель смешанных эффектов с отрицательной дисперсией
Я знаю, что это старый вопрос, но мне интересно, есть ли у него решение сейчас. Я обычно использую модель смешанных эффектов с lme4 пакетом с функцией lmer . Однако я знаю, что эта функция не позволяет мне включать в модель компоненты...
604 просмотров
schedule 17.10.2022

Как воспроизвести случайные эффекты в lme4 от SAS?
Я хочу запустить линейную смешанную модель для зависимой переменной DV , которая собирается под двумя разными Condition в трех разных Timepoint . Данные структурированы следующим образом: ## dput(head(RawData,5))...
205 просмотров
schedule 09.06.2024

Предупреждающие сообщения glmmLasso
Я пытаюсь запустить glmmLasso для оценки смешанной модели с помощью команды: glm1_final <- glmmLasso(Activity~Novelty + Valence + ROI, rnd = list(Subject=~1), data = KNov, lambda=lambda[opt],switch.NR=F,final.re=TRUE) Этот код в...
1763 просмотров
schedule 12.03.2023

Начальные значения смешанной модели для lme4
Я пытаюсь подогнать под смешанную модель с помощью функции lmer из пакета lme4 . Однако я не понимаю, что следует вводить в параметр start . Моя цель - использовать простую линейную регрессию, чтобы использовать оцененные там коэффициенты в...
2263 просмотров
schedule 28.10.2023

Функция lme4::glmer.nb выдает ошибку в family$family: оператор $ не определен для этого класса S4 в зависимости от порядка запуска моделей
library(lme4) dummy <- as.data.frame(cbind(speed = rpois(100, 10), pop = rep(1:4, each = 25), season = rep(1:2, each = 50), id = seq(1, 100, by = 1))) dummy2 <- as.data.frame(cbind(speed = c(rnbinom(50, 10, 0.6), rnbinom(50, 10, 0.1)), pop =...
565 просмотров
schedule 06.12.2022

Изменение CI для difflsmeans
Я запускаю линейную смешанную модель и использую функцию diffslmeans() из пакета lmerTest для вычисления попарных различий в моих фиксированных эффектах. Функция, кажется, производит 95% CI независимо от того, что я пытаюсь. Однако я хотел бы...
183 просмотров
schedule 08.12.2023

Общая предсказательная сила (например, R2) для байесовских линейных (смешанных) моделей
Следуя этому вопрос о журнальной отчетности, я хотел бы знать, есть ли какой-либо общий показатель размера эффекта для байесовских моделей, подобранных с использованием stan_lmer ? В частотной структуре есть, например, псевдо-R2 (вычисленный с...
186 просмотров
schedule 30.03.2022

Большие данные: обобщенные линейные модели со смешанными эффектами
Я ищу предложения по стратегии подбора обобщенных линейных моделей со смешанными эффектами для относительно большого набора данных. Допустим, у меня есть данные о 8 миллионах передач баскетбола в США примерно 300 командами за 10 лет. Данные...
1108 просмотров
schedule 06.05.2022

усредненные по модели коэффициенты линейных смешанных моделей в glmulti? Исправление больше не работает
Я использую пакет glmulti для выбора переменных фиксированных эффектов смешанной модели в lme4 . У меня была та же проблема с получением коэффициентов и доверительных интервалов, которую решил автор пакета в этой теме . А именно использование...
225 просмотров
schedule 20.01.2024

Случайные эффекты в моделях продольного многоуровневого импутации с использованием MICE
Я пытаюсь вменять данные в набор данных с продольным планом. Есть два предиктора (экспериментальная группа и время) и одна переменная результата (оценка). Переменная кластеризации - id. Вот данные игрушки set.seed(345) A0 <- rnorm(4,2,.5) B0...
1155 просмотров
schedule 14.06.2023

Предупреждение lme4: Модель не сходится с max | grad |
Мне нужно запустить lmer с преобразованной логарифмической переменной ответа, непрерывной переменной в качестве фиксированного эффекта и вложенным случайным эффектом: first<-lmer(logterrisize~spm + (1|studyarea/teriid), data =...
16454 просмотров
schedule 10.03.2022

Смешанная модель логита в Zelig (R), которая больше не работает, недоступна?
Я заинтересован в получении первых отличий от модели смешанного логита с использованием пакета Zelig. Однако я не могу запустить модель смешанного логита в Zelig. Я обновил пакет Zelig в соответствии с инструкциями веб-сайта Zelig. Я запустил...
319 просмотров

Как загрузить R-квадрат смешанной модели?
Я пытаюсь получить загрузочный R ^ 2 для модели смешанного эффекта. Поскольку уже существует только обходной путь для получения условного и предельного R^2, я попытался загрузить эту статистику на основе примера, приведенного методы статистики для...
223 просмотров
schedule 15.09.2023

Как написать формулу Лмера для модели смешанных эффектов с двумя фиксированными эффектами
Я новичок в линейных моделях смешанных эффектов и пытаюсь использовать их для проверки гипотез. В моих данных ( DF ) у меня есть две категориальные / факторные переменные: color (красный / синий / зеленый) и direction (вверх / вниз). Я хочу...
274 просмотров
schedule 23.03.2022

Прогнозирование с помощью lmer: ошибки и отсутствующие p-значения
В эксперименте самки рыб подвергались воздействию двух уровней фотопериода (окружающий и сжатый) и двух уровней температуры (4 и 7). Они находились в четырех резервуарах (два резервуара для каждого фотопериода, по одному резервуару для каждой...
192 просмотров
schedule 21.02.2022