Вопросы по теме '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 просмотров
schedule
01.08.2023
Как загрузить 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