Вопросы по теме 'lme4'

Есть ли способ создать таблицу LaTeX из подходящего объекта модели lme4 mer?
Кто-нибудь знает способ создать LaTeX-таблицу хорошего качества публикации из объекта lme4 mer ? Ни метод xtable (пакет xtable ), ни метод latex (пакет Hmisc ) не знают, как работать с mer объектами. Например, при такой посадке:...
7017 просмотров
schedule 16.10.2022

Ошибка в glmer: NA/NaN/Inf в вызове сторонней функции (arg 1)
Я пытаюсь соответствовать модели glmer(trans.dummy ~ pop + year + (year | munname), data=pool, family=binomial(link = "logit"), REML=T, verbose=T) но я продолжаю получать следующую ошибку: Error in glm.fit(fr$X, fr$Y, weights =...
11297 просмотров
schedule 28.06.2022

R: найти исходный код, который производит вывод для слота S4?
Добрый день Всем, Когда функция «lmer» в «lme4» запускается, она создает объект S4 с большим количеством слотов. Меня интересует один из этих слотов, а именно модель@X, и то, как создается вывод этого слота «X». Я хочу попытаться воспроизвести...
339 просмотров
schedule 25.05.2023

Оценка функции правдоподобия в линейных смешанных моделях (lme4)
В настоящее время я пишу сценарий для оценки (ограниченной) функции логарифмического правдоподобия для использования в линейных смешанных моделях. Мне это нужно для расчета вероятности модели с некоторыми параметрами, фиксированными на произвольные...
1382 просмотров
schedule 16.01.2024

Использование вывода объекта S4 в исходных пакетах
В последнее время я пытался организовать все свои функции в пакетах, чтобы не создавать сотни строк кода поверх каждого скрипта. Я написал ряд функций, которые использую после запуска набора многоуровневых моделей LME4. Эти функции предназначены для...
505 просмотров
schedule 20.05.2022

nlme и lme4 Игнорирование членов в квадрате
Я пытаюсь использовать стандартную функцию спроса на транслог, а именно: lnY = lnP + lnZ + lnY*lnZ + lnY^2 + lnZ^2 Где Y = спрос, P = цена и Z = доход. Однако, когда я включаю элементы в квадрате в nlme или lme4, они просто...
3300 просмотров
schedule 27.06.2022

Сообщение об ошибке в lme4::glmer: «что» должно быть строкой символов или функцией
У меня многоуровневая модель. Я использую следующие команды с validatedRS6 в качестве результата, random в качестве предиктора и clustno в качестве переменной случайных эффектов. new<-as.data.frame(read.delim("BABEX.dat", header=TRUE))...
1901 просмотров
schedule 15.06.2023

Как получить и подобрать прогнозы от glmer?
Я попробовал несколько предложений здесь и на других сайтах, но ни одно из них не работает. Моя модель выглядит так model = glmer(cbind(live, dead)~cont1*cont2+ (1|random), family = binomial, data = all) cont1 и cont2 —...
1076 просмотров
schedule 05.01.2023

сравнение коэффициентов между lmList и lmer
может ли кто-нибудь сказать мне, почему коэффициенты наклона отличаются между коэффициентами, извлеченными из модели lmer со случайным наклоном, и коэффициентами из модели lmList, подходящей для того же набора данных? Спасибо...
673 просмотров
schedule 05.03.2022

Как отменить масштабирование коэффициентов из lmer () - модели с масштабированным откликом
Я приспособил модель в R с lmer() -функцией из пакета lme4 . Я масштабировал зависимую переменную: mod <- lmer(scale(Y) ~ X + (X | Z), data = df, REML = FALSE) Я...
5965 просмотров
schedule 22.02.2022

R - Спецификация модели для повторных измерений GLMM (lme4)
У меня возникли проблемы с правильным указанием моей продольной модели в R. Мой анализ рассматривает гендерные различия в баллах, оцениваемых в трех временных точках. Фактически, я хочу увидеть, имеет ли тот или иной пол постоянно более высокие баллы...
5363 просмотров
schedule 19.04.2022

Сохранение постоянной дисперсии случайных эффектов в PROC MIXED по сравнению с lmer
Мне было интересно, можно ли поддерживать постоянную дисперсию случайных эффектов в функциях R's lme или lmer (или другой подпрограмме случайных эффектов в R ) или, по крайней мере, указать начальные значения. В SAS это возможно с помощью...
262 просмотров
schedule 11.05.2024

Прогнозирование с lme4 на новых уровнях
Я пытаюсь подогнать под модель смешанных эффектов, а затем использовать эту модель для создания оценок для нового набора данных, который может иметь разные уровни. Я ожидал, что оценки нового набора данных будут использовать среднее значение...
13716 просмотров

Как правильно учесть вложенные случайные эффекты с lme4?
У меня есть фрейм данных с переменными subject , wd и group и переменная ответа value . Каждый испытуемый разделен на одну группу, и каждый будний день выполняет 7 измерений. Поскольку каждый субъект полностью вложен в группу, я хочу...
1770 просмотров
schedule 28.07.2023

Несовместимость R stargazer, lme4 и lmerTest
R новичок здесь. У меня проблемы с работой lmerTest и stargazer . Я следовал инструкциям здесь, чтобы заставить stargazer работать с lme4 в R....
6378 просмотров
schedule 04.04.2022

Как я могу соответствовать регрессии Скеллама?
Есть ли простой способ подогнать многомерную регрессию в R, в которой зависимая переменная распределена в соответствии с распределением Скеллама (разница между двумя счетчиками с распределением Пуассона)? Что-то типа: myskellam <- glm(A ~ B...
985 просмотров
schedule 19.11.2022

Как указать разные случайные эффекты в nlme и lme4?
Я хочу указать различные случайные эффекты в модели, используя nlme::lme (данные внизу). Случайные эффекты: 1) intercept и position варьируются в пределах subject ; 2) intercept изменяется от comparison . Это просто использовать...
3855 просмотров
schedule 21.02.2022

Ошибка: объект «sigma» не экспортируется «пространством имен: lme4»
После обновления lme4 (1.1-12) я получаю эту ошибку при вызове пакета varComp . library(varComp) Error : object ‘sigma’ is not exported by 'namespace:lme4' Error: package or namespace load failed for ‘varComp’ Отредактировано Вот мой...
2754 просмотров
schedule 06.10.2022

Как получить ковариационную матрицу для случайных эффектов (BLUP/условные режимы) из lme4
Итак, я подогнал линейную смешанную модель с двумя случайными перехватами в R: Y = X beta + Z b + e_i, где b ~ MVN (0, Sigma) ; X и Z — матрицы модели фиксированных и случайных эффектов соответственно, а beta и b — параметры...
2856 просмотров
schedule 16.07.2023

Построение графика взаимодействия двух непрерывных переменных с данными lme4
Я пытаюсь изобразить взаимодействие между двумя непрерывными переменными в R. Однако мои данные многоуровневые (люди вложены в течение нескольких дней), поэтому мне нужно учитывать вложенную структуру моих данных, когда я строю ее график. Я...
1444 просмотров
schedule 25.01.2023