Пакет R для подбора статистических моделей нелинейных смешанных эффектов.
1 ответ

Пропущенные значения в объекте в R lme

Я работаю над простой моделью в R: fit<-lme(x~y, data, random=~1|subject) но я продолжаю получать эту ошибку: Ошибка в na.fail.default (список (h1 = c(103L, 20L, 34L, 85L, 47L, 136L, 76L,: пропущенные значения в объекте (h1 - одно из названий мое…
11 апр '16 в 21:29
0 ответов

Альтернативные решения для многоуровневого моделирования из-за того, что на уровнях 2 и 3 недостаточно кластеров

Я использовал несколько многоуровневых моделей для своих данных, но теперь мне сказали, что мои кластеры слишком малы, и мне нужно рассмотреть возможность вынести два уровня и добавить их в качестве ковариат (фиктивные переменные). Я сделал это, и э…
0 ответов

Перекодирование SAS PROC MIXED MODEL в R

Я пытаюсь перекодировать смешанную модель из SAS в R, и мне трудно, так как SAS не является языком кодирования, с которым я когда-либо работал. Я пытался использовать три документа, чтобы помочь мне: Резюме смешанных моделей в SAS и R lmer для смеша…
07 янв '18 в 16:42
1 ответ

ggplot2 выводит два набора данных в одну картинку

Это должно быть FAQ, но я не могу найти точно такой же пример в других ответах (не стесняйтесь закрыть это, если вы можете указать аналогичные вопросы и ответы). Я все еще новичок в ggplot2 и, похоже, не могу так легко обернуться вокруг него. У меня…
07 мар '16 в 16:15
0 ответов

Укажите модель ANOVA внутри субъекта и между субъектами, используя lme или lmer, в качестве фиксированных эффектов

Есть ли способ указать двусторонний ANOVA, с одним предиктором внутри субъекта и одним предиктором между субъектами, используя lme (из nlme) или lmer (из lme4)? Возможно, это перекрестный вопрос, но я не думаю, что меня интересуют случайные эффекты,…
04 авг '17 в 03:36
0 ответов

Возможная ошибка: R Проверка CMD завершается неудачно в Ubuntu 18.04 LTS при использовании nlme

С помощью devtools::check() выдает следующее сообщение об ошибке: R CMD check results 1 error | 0 warnings | 0 notes checking whether package ‘abc’ can be installed ... ERROR Installation failed. See ‘/tmp/RtmpSpjP2M/abc.Rcheck/00install.out’ for de…
20 авг '18 в 08:29
0 ответов

nlme() не запускается при использовании в пользовательской функции в R

Я пытаюсь создать пользовательскую функцию, используя nlme(), и я не могу понять, почему моя функция не запускается... Я предполагаю, что nlme не может получить доступ к объектам в локальной среде моей функции. Чтобы дать вам некоторое представление…
04 сен '17 в 14:48
0 ответов

Различия между p-значениями nlm и доверительными интервалами

У меня есть набор данных, в котором я оценил размер домашнего диапазона, используя 10 оценщиков из 41 человека. Я хотел проверить, насколько эти оценки существенно отличаются друг от друга, поэтому я подгоняю линейную модель смешанных эффектов к nlm…
25 июн '14 в 18:18
1 ответ

Сбой glmmPQL при включении объекта corSpatial

Ссылка на данные (1170 объектов, 9 переменных, файл.Rd) Просто прочитайте это, используя readRDS(file), Я пытаюсь настроить GLMM, используя glmmPQL функция от MASS пакет, включающий часть случайных эффектов и учет пространственной автокорреляции. Т…
21 июл '16 в 16:03
1 ответ

R эквивалентный случайный остаток по субъекту в SAS

Я могу закодировать эту проблему в SAS с остаточным значением как случайный эффект (я считаю, что это случайный перехват r-стороны рыбой) proc glimmix data=one method=mmpl ; class fish; model increment =age growth_year age*growth_year; random residu…
23 мар '18 в 20:57
0 ответов

Модель смешанных эффектов в R с использованием lme

Я анализирую данные для наблюдательного исследования роста форели. Таким образом, в датафрейме много NA. Кроме того, разные методы лечения имеют разное количество наблюдений. Так я думаю это можно назвать несбалансированным дизайном? Year Site FishI…
23 апр '13 в 13:41
1 ответ

Построение прогнозируемых значений по модели Ime (с полиномами) в R

Я использую линейную модель смешанного эффекта (запустить с lme() функция в пакете nlme в R), которая имеет один фиксированный эффект и один случайный член перехвата (для учета разных групп). Модель представляет собой кубическую полиномиальную модел…
31 янв '15 в 17:26
1 ответ

Будут ли переменные случайного эффекта автоматически приниматься как факторы в lmer (или lme) в R?

Я понимаю, что иметь непрерывную или числовую переменную в качестве случайного эффекта в модели смешанных эффектов не имеет особого смысла (например, см. Здесь). Но что мне интересно, если lme4::lmer или же nlme::lme в R целенаправленно мешают вам с…
10 окт '16 в 16:30
1 ответ

Почему линейная смешанная модель работает в SAS и nlme, но не в lme4?

Мои данные состоят из 20 предметов в контрольной группе и 20 в экспериментальной группе. Представляющий интерес DV представляет собой оценку изменения пиковой мощности, измеренную для каждого участника. Существует также фиктивная переменная xVarExp …
30 июн '16 в 02:26
0 ответов

Нелинейная регрессия с группами по nlsList в R

Возьмите пример mtcars в R: я хочу сделать нелинейную регрессию для отношения между mpg и disp для каждого цилиндра и каждой группы против. В моем случае у меня есть кадр данных, намного больший, чем mtcars. Я знаю, что недостаточно данных. Но чтобы…
19 дек '16 в 14:01
0 ответов

Укажите модель нелинейных смешанных эффектов для степенного закона в nlme

Я хочу использовать nlme пакет, чтобы соответствовать степенному закону в смешанных эффектах для моих данных. Я был в состоянии использовать brms чтобы соответствовать байесовской модели, см. ниже, но я хочу соответствовать "частой" версии той же мо…
06 июн '18 в 14:08
0 ответов

Проблема с моделью GLS, ошибка сходимости

Я хотел бы сделать модель GLS, используя эти векторы: v1 <- c( 3.27985464, -0.40586830, -0.40586830, -0.09872472, -0.09872472, -0.40586830, -0.40586830, 0.82270601, -0.40586830, -0.40586830, -0.25229651, -0.40586830, -0.40586830, -0.40586830) yea…
22 авг '17 в 13:42
0 ответов

Указанная ковариационная матрица в lmer в R для штрафных сплайнов

Я пытался вписать регрессионную модель штрафованных сплайнов в R, используя связь между штрафованными сплайнами и линейными смешанными моделями. Хотя я довольно знаком с функцией R lme, используя своего конкурента lmer представляет некоторые труднос…
26 апр '18 в 15:04
1 ответ

Я с переменным для повторных мер

Я столкнулся со странными результатами при использовании Ime с повторными измерениями и VarIdent. Любая помощь с этим будет очень признательна! Я проверяю, отличается ли сигнал 13C листьев по временному ряду между двумя видами (A и B). Меня в основн…
23 сен '16 в 12:04
0 ответов

Ошибка документации пакета NLME: суммы квадратов anova.gls() и anova.lme()

Документация по функциям nlme anova.lme() а также anova.gls() четко говорится: "Когда присутствует только один подобранный объект модели, возвращается фрейм данных с суммами квадратов...". При выполнении приведенного ниже примера сумма квадратов не …
12 окт '16 в 09:46