Вопросы по теме '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 просмотров
schedule
11.03.2023
Стандартная ошибка компонента дисперсии на выходе 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