NoneLme4 - это пакет R для подбора и анализа линейных, нелинейных и обобщенных линейных смешанных моделей.
1 ответ

GLMER: ошибка: (maxstephalfit) разделение на шаги PIRLS не смогло уменьшить отклонение в pwrssUpdate

Я изучаю влияние различных характеристик на решение суда по конкретным преступлениям. Набор данных довольно большой (28928 наблюдений с 86 единицами уровня 2). Я смотрю на решение, заключать ли кого-либо в тюрьму или нет (= бинарная переменная резул…
24 июн '16 в 14:19
0 ответов

Как запретить модулю оптимизатора в lme4 изменять глобальные переменные?

Я использую lmer4 проверить процедуру начальной загрузки. Я надеялся воспользоваться его модульной структурой, изменив только некоторую часть выходных данных первого модуля (lFormula функция). Однако функция оптимизатора, похоже, изменяет вывод lFor…
19 ноя '17 в 13:30
2 ответа

Ищете специальную модель для Anova типа 2 (модель со смешанными эффектами)

У меня вопрос, как мне запустить post hoc на Anova типа 2 (модель со смешанными эффектами)? Пока я использую glmer() из пакета "lme4", Anova() из пакета "автомобиль", и пытается выполнить тест HSD из пакета "Agricolae". После поиска в течение некото…
19 июн '18 в 19:18
1 ответ

Используя as.formula() в glmer и углубляя его в MuMIn

Я пытаюсь использовать dredge функция от MuMIn пакет, как следует. Поскольку я хочу запустить его для разных наборов данных, я использую функцию, которая принимает вектор ковариат, и формула создается из этого вектора внутри функции. require(lme4);r…
16 ноя '13 в 22:22
1 ответ

anova.mer Mod(mod_null, mod_eins): модели не были все приспособлены к одному и тому же размеру набора данных

Я видел, что у некоторых ребят тоже была эта проблема, но я действительно не понимал данных ответов. Я сделал несколько линейных смешанных моделей, начиная с модели "только перехват". Впоследствии я хотел добавить больше переменных. Когда я пытался …
22 мар '17 в 14:51
0 ответов

Нелинейная смешанная модель на основе исходных данных

Я пытаюсь построить результирующую кривую из подгонки нелинейной смешанной модели. Это должно быть что-то вроде кривой нормального распределения, но с перекосом вправо. Я перешел по предыдущим ссылкам здесь и здесь, но когда я использую свои данные,…
16 окт '16 в 23:32
0 ответов

GLMM-скрипт с фиксированными и случайными переменными

Я пытаюсь создать модель в R, которая учитывает мой опрос. В исследовании участвуют многие объясняющие переменные, переменная ответа называется HUNTED и данные о 3 различных времени PERIODS, Мой набор данных слишком большой, я только что добавил его…
05 сен '18 в 19:45
0 ответов

Понимание структуры и синтаксиса случайных эффектов LME4

Вот проблема, которую я пытаюсь понять. Давайте предположим, что я провел относительно упрощенный эксперимент, в котором у меня есть несколько субъектов, и я записываю активность ЭЭГ на коже головы, которую я разделил на 9 областей. На самом базовом…
24 июн '17 в 16:40
1 ответ

lmertest и доверительные интервалы

Есть ли способ получить lmertest для получения доверительных интервалов? Я пытался использовать этот код из https://cran.r-project.org/web/packages/lmerTest/lmerTest.pdf contest(mymodel , L, rhs = 0, joint = TRUE, collect = TRUE, confint = TRUE, lev…
22 май '18 в 06:27
2 ответа

Вычисление p-значений для модели с нулевым случайным эффектом в lm4/lmerTest

Edit1: как сообщается @RolandASc ниже, похоже, ошибка в lmerTest , Я уже написал электронное письмо сопровождающему пакета с сообщением о проблеме. Edit2: Ответ сопровождающего: "мы работаем над обновлением, где, как мы надеемся, такие проблемы буду…
23 янв '18 в 09:18
0 ответов

Подгонка повторяющегося однофакторного эксперимента с gls, lme и gam

Базовый эксперимент: CRD с 4 обработками (N,L,M и ML) и 3 повторениями / трт. Итак, всего 12 прогонов / испытаний / предметов. Концентрацию каждого прогона определяли в день 0, 1,2,6,7, 9,10, 13 и 13, что в итоге составило 9 × 12 = 108 ак / точек да…
19 окт '18 в 20:30
0 ответов

Что означает это сообщение: "isS4" не является встроенной функцией?

Я установил lme4 пакет и попытался запустить example(glmer) но получил сообщение об ошибке: Error in function (..., deparse.level = 1) : "isS4" is not a BUILTIN function Тем не менее, когда я бегу builtins()[590], R возвращает "isS4"что заставляет м…
08 окт '13 в 01:42
1 ответ

rpy2: представляет NA в Python как аргумент функции R

Я пытаюсь передать NA для функции R, например, делать прогнозы с lme4 Смешанная модель, использующая только фиксированные эффекты (т.е. без случайных эффектов): import rpy2.rinterface as ri from rpy2.robjects.packages import importr rstats = importr…
24 фев '19 в 04:05
1 ответ

Ошибки при попытке запустить glmer с парными данными и повторными измерениями (lme4)

Я хочу проанализировать влияние двух обработок на изменение численности видов растений по временному градиенту. План эксперимента состоит из исключений (обработка = лось), в сочетании с контрольными участками (обработка = лось). Сайт состоит из 1 ис…
21 авг '17 в 18:06
0 ответов

Сообщение об ошибке при создании прогнозов для пользовательской функции связи для glmer

Я пытаюсь смоделировать выживание гнезда, используя функцию glmer в R. Чтобы включить дни воздействия в модель, я использовал отличную пользовательскую функцию связи Бена Болкера ( здесь). Однако, когда я пытаюсь сгенерировать прогнозы для каждого р…
17 авг '16 в 12:37
0 ответов

Извлечь остатки из цикла lme4

Я хочу извлечь остатки из линейного смешанного модельного цикла. results <- matrix(NA, nrow(meth),4, dimnames=list(rownames(meth), c("nobs", "coef", "coef.se", "Pval"))) Цикл создает матрицу с результатами для переменной X, но теперь я хочу извле…
18 июн '17 в 10:06
0 ответов

Странная картина не сходимости при сглаживании в смешанных эффектах логистической регрессии

Я провожу симуляционное исследование эффектов добавления дробных чисел успехов и неудач, которые я назову C, к логистическим регрессиям со смешанными эффектами. Я смоделировал 2000 наборов данных и смоделировал каждый с 5 логистическими регрессиями …
31 июл '15 в 01:01
0 ответов

Эффекты взаимодействия

У меня есть несколько вопросов, касающихся 3-стороннего эффекта взаимодействия в линейных смешанных модельных анализах в R с пакетом lme4. Моим результатом оценки является "MMSE", который является когнитивным тестом. В основном меня интересует сниже…
01 ноя '17 в 06:38
0 ответов

Модельные матрицы в R, для смешанных моделей

Я пишу свой собственный программный пакет для определенного типа смешанных моделей в R и хотел бы простой и эффективный способ формирования матриц проектирования. Фиксированная часть ("X") не проблема, я использую model.matrix, или же sparse.model.m…
11 авг '17 в 19:48
1 ответ

Почему функция Предиката () не принимает работу "новые данные" для моей модели lmer?

Когда я пытаюсь использовать Предикат (), или ПредикатИнтервал (), или bootMer(), предоставляя ему новые данные (даже если это те же данные, с которыми я построил модель), я получаю следующую ошибку: Error in [.data.frame(fr, vars) : undefined colum…
07 мар '17 в 14:46