Вопросы по теме '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