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

Ошибка при попытке запустить lmer() в R
Итак, вот моя проблема. У меня есть набор данных в R, на котором мне нужно запустить модель смешанных эффектов. Вот код: data <- read.csv("D:/blahblah.csv") analysis.data <- lmer(intdiff ~ stress_limit * word_position * follows +...
11600 просмотров
schedule 12.04.2023

вывод anova пуст при запуске на модели lmer
Я пытаюсь запустить модель анова, используя библиотеку lme4. У меня нет реальных данных, поэтому я просто взял данные Cars93 из MASS. Хотя я мог заставить lmer(...) работать, я не могу распечатать таблицу anova. Вот что произошло: >...
903 просмотров
schedule 11.03.2023

r: package arm: ошибка при загрузке: обнаружены функции при экспорте методов из пространства имен «arm», которые не являются общими для S4: «fixef», «ranef»
> library(arm) Loading required package: lme4 Loading required package: RcppEigen Loading required package: R2WinBUGS Loading required package: foreign Error : Functions found when exporting methods from the namespace ‘arm’ which are not S4...
2667 просмотров
schedule 26.03.2024

Как построить отдельные траектории модели lme
У меня есть пример данных и модель x<-rep(seq(0,100,by=1),10) y<-15+2*rnorm(1010,10,4)*x+rnorm(1010,20,100) id<-NULL for (i in 1:10){ id<-c(id, rep(i,101))} dtfr<-data.frame(x=x,y=y, id=id) library(nlme) with (dtfr,...
2578 просмотров
schedule 22.01.2023

glmer - прогнозирование с использованием биномиальных данных (данные подсчета cbind)
Я пытаюсь предсказать значения с течением времени (дни по оси x) для модели glmer, которая была запущена на моих биномиальных данных. Total Alive и Total Dead - это данные подсчета. Это моя модель и соответствующие шаги ниже....
18493 просмотров
schedule 15.12.2023

Невозможно запустить lmer из функции
Я столкнулся с проблемой, пытаясь встроить lmer в функцию. Вот воспроизводимый пример с использованием данных из lexdec . Если я запускаю lmer напрямую во фрейме данных, проблем не возникает. Например, предположим, что я хочу узнать,...
841 просмотров
schedule 14.07.2022

R: Странные ошибки при использовании lmList от lme4
Я читаю эту книгу о данных. анализ и графика в R, и я думаю, что есть некоторые проблемы между текущей версией lme4 и той, которую они использовали (у меня есть версии 2007 и 2010). Моя текущая проблема касается страницы 340 книги, где авторы...
441 просмотров
schedule 26.04.2023

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

прогнозировать lme4 для ответа типа (двоичный выбор) Доверительные интервалы
Я провел логистическую регрессию со смешанными эффектами с lme4 (type = "response"). Теперь я использовал функцию прогнозирования и хотел также определить доверительные интервалы. Я нашел этот код ( http://glmm.wikidot.com/faq ) для прогнозов, он...
463 просмотров
schedule 20.12.2022

Просмотр стандартных ошибок и оценок параметров в lme4
Я только начал работать с пакетом lme4 для линейных смешанных моделей в R. Я начал следовать книге Дугласа Бейтса, чтобы понять процесс. В книге вывод модели, подходящей с использованием набора данных «sleepstudy» в пакете, дается как: Linear...
1103 просмотров
schedule 18.11.2022

подбор линейной смешанной модели к очень большому набору данных
Я хочу запустить смешанную модель (используя lme4::lmer ) для 60 млн наблюдений следующего формата; все предикторы/зависимые переменные являются категориальными (факторами), за исключением непрерывной зависимой переменной tc ; patient — это...
4132 просмотров

Стандартная ошибка компонента дисперсии на выходе lmer
Мне нужно извлечь компонент дисперсии standard error из вывода lmer . library(lme4) model <- lmer(Reaction ~ Days + (1|Subject), sleepstudy) Следующее дает оценки компонента дисперсии: s2 <- VarCorr(model)$Subject[1] Это НЕ...
7965 просмотров
schedule 20.02.2024

Другая версия lme4 дает другой результат
Я столкнулся с проблемой в пакете lme4 , которая кажется мне очень странной. R версии 3.0.2 (2013-09-25) и lme4_1.1-7 ведут себя по-разному с R версией 3.2.1 (2015-06-18) и lme4_1.1-8 . Например, в R версии 3.0.2 при запуске кода...
174 просмотров
schedule 23.07.2023

Первые линейные контрасты: Кенвард Роджерс или Саттертуэйт DF и SE
В R я ищу способ оценить доверительные интервалы для линейных контрастов для lmer моделей, которые используют либо степени свободы Кенварда-Роджера, либо Саттертуэйта и SE. Например, я могу вычислить CI для параметра с фиксированным эффектом в...
1654 просмотров
schedule 10.07.2022

Извлечение оценок уникальных коэффициентов из glm
У меня есть 57 особей (с ~ 60 000 повторений), помеченных «ID_Year», которые используются для создания модели выбора среды обитания GLME. Как извлечь уникальные оценки коэффициентов для каждого "ID_Year"? Вот модель: Cand.mod[[25]]...
58 просмотров
schedule 06.07.2023

Ошибка в qlogis(p) в lme4 с семейством logit
Я пытаюсь подогнать GLMM к данным о размерах (длине тела) личинок рыб из разных регионов при разных температурах. Каждая комбинация регион-температура имеет измерения как минимум для 10 и до 30 человек. Рассматриваемые данные были собраны в течение...
981 просмотров
schedule 24.11.2022

Мультиколлинеарность в линейной смешанной модели
Я подгоняю следующие модели lmer в R: lmer(rt~deadline*cond+age+(1+deadline|task/pp), REML=FALSE) ##Model 1 lmer(rt~deadline+cond+age+(1+deadline|task/pp), REML=FALSE) ##Model 2 Где rt - время реакции, deadline имеет 2 уровня, короткий...
727 просмотров
schedule 15.12.2022

Ошибка в r.squaredGLMM()
Я создаю GLMM (используя glmer() из пакета R "lme4"), и иногда я получаю сообщение об ошибке при оценке значений R2 (используя r.squaredGLMM() из пакета "MuMIn"). Модель, которую я пытаюсь подогнать, похожа на эту: library(lme4) lmA <-...
622 просмотров
schedule 29.07.2023

Как мне преодолеть эту ошибку при выполнении функции lmer: Ошибка в if (REML) p else 0L: аргумент не интерпретируется как логический?
Я проверил защитные методы, как обсуждалось в это сообщение , чтобы предотвратить эту ошибку, но она все еще не исчезла. model<-lmer(Proportion~Plot+Treatment+(1|Plot/Treatment),binomial,data=data) Ошибка в if (REML) p else 0L:...
1730 просмотров
schedule 01.06.2023

Построение наклона с фиксированными эффектами из модели lmer
У меня есть следующий набор данных: > dput(df) structure(list(Subject = c(1L, 2L, 3L, 5L, 6L, 6L, 6L, 7L, 7L, 7L, 8L, 8L, 8L, 9L, 9L, 9L, 10L, 10L, 11L, 11L, 11L, 12L, 12L, 13L, 13L, 14L, 14L, 15L, 15L, 16L, 16L, 16L, 17L, 17L, 17L,...
602 просмотров
schedule 01.03.2023