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

Извлечение RMSE из объекта lmList
У меня есть следующий код для извлечения коэффициентов из объекта lmList. library(lme4) library(plyr) lm.model <- lmList(Y ~ X | eventID, df) param <- ldply(lm.model, coef) Этот код работает нормально. Однако является ли это...
1656 просмотров
schedule 03.12.2022

Получение обобщенных средств наименьших квадратов для фиксированных эффектов в 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

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

Частотные веса в R (многоуровневый с использованием lme4)
Я строю многоуровневую модель с lme4 (два уровня: класс и уровень ученика с образовательными данными). Я использую данные ICCS 2009: http://www.iea.nl/iccs_2009.html ). ICCS 2009 рекомендует использовать гири. Я сделал «внутришкольный вес...
1380 просмотров
schedule 17.10.2023

Ошибка совместимости lme4 и languageR: "модель ввода не является объектом mer"
У меня есть набор данных dat2, на котором я хочу подогнать линейную модель смешанных эффектов. Раньше я использовал lmer () (пакет lme4) в дополнение к pvals.fnc для вычисления связанных p-значений. Однако, поскольку я переустановил версию R 3.0.2...
3272 просмотров
schedule 18.08.2022

неясное предупреждение lme4 с использованием lmer в optwrap
Используя lmer , я получаю следующее предупреждение: Warning messages: 1: In optwrap(optimizer, devfun, x@theta, lower = x@lower) : convergence code 3 from bobyqa: bobyqa -- a trust region step failed to reduce q Эта ошибка возникает после...
3278 просмотров
schedule 14.04.2024

Как использовать lme4.0 с lmerTest?
Мне интересно, сталкивался ли кто-нибудь с этой проблемой раньше. Я использую пакет lmerTest для запуска моделей смешанных эффектов в R, потому что он имеет удобный способ предоставления p-значений. Этот пакет по умолчанию загружает самую...
1173 просмотров
schedule 22.06.2022

модель glmer с начала 2013 года: предупреждающее сообщение о сходимости при повторном запуске
Более 1 года назад (февраль 2013 г.) я использовал lmer для запуска модели со смешанным эффектом, включающей биномиальный результат, с помощью следующей команды: nl3.lmer <- glmer( cul.bi ~ food.act + where + intlan + inter.cul + via.m +...
2108 просмотров
schedule 16.04.2022

функция drop1 для lmer
Я построил модель смешанных эффектов с тремя фиксированными эффектами и одним случайным эффектом. mdl1 <- lmer(yld.res ~ veg + rep + rip + (1|state), REML=FALSE,data=data2) Я хочу получить самую экономную модель из вышеописанной...
7058 просмотров
schedule 06.11.2023

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

Невозможно запустить модели glmer с na.action = na.fail, что необходимо для функции земснаряда MuMIn
Mac OS 10.9.5, R 3.2.3, MuMIn_1.15.6, lme4_1.1-10 Воспроизводимый пример кода с использованием пример данных руководство пользователя MuMIn рекомендует использовать na.action=na.fail , в противном случае функция земснаряда будет не...
3315 просмотров
schedule 19.01.2024

Извлечение p-значений для фиксированных эффектов из вывода nlme / lme4
Я пытаюсь извлечь отдельные элементы (в частности, p-значения) из таблицы фиксированных эффектов, содержащейся в объекте, созданном сводным вызовом модели смешанных эффектов. Данные игрушек: set.seed(1234) score <- c(rnorm(8, 20, 3),...
5189 просмотров
schedule 01.05.2022

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

Начальные значения смешанной модели для 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

Невозможно настроить средство для глмера с помощью lsmeans
У меня есть glm, который я хотел бы настроить для использования lsmeans. Следующий код создает модель (и, кажется, делает это правильно): library(lmerTest) data$group <- as.factor(data$grp) data$site <- as.factor(data$site) data$stimulus...
895 просмотров
schedule 16.11.2023

Ошибка при использовании пакета посредничества с моделью lme4: модель посредника еще не реализована - предыдущий ответ не работает для меня
Привет, я новичок в R и на этом форуме. У меня та же проблема, что и у предыдущий постер , но мой не нашел ответа. Я получаю сообщение: «Модель посредника еще не реализована», когда пытаюсь запустить многоуровневую модель посредничества с...
1596 просмотров
schedule 10.10.2022

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

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

линейная смешанная модель на продольных данных: время как непрерывный фактор против неупорядоченного фактора против упорядоченного фактора
Здесь я сделал некоторые продольные данные. Переменным результатом является уровень биомаркеров у пациентов. Переменные посещения представляют собой метки посещений. Переменное время означает дни от базового уровня t1. Есть два статуса ответа:...
264 просмотров
schedule 15.01.2023