Средние наименьших квадратов с их стандартными ошибками для объекта aov можно получить с помощью функции model.tables : npk.aov <- aov(yield ~ block + N*P*K, npk) model.tables(npk.aov, "means", se = TRUE) Интересно, как получить...

Я немного погуглил и не могу найти документацию по этому поводу. Я пытаюсь оценить возможную модель обобщенных наименьших квадратов (FGLS) для данных временных рядов поперечного сечения в R. Например: library(nlme) foo <- gls(Y ~...

Я пытаюсь запустить линейную модель со смешанными эффектами (пакет nlme), но постоянно получаю сообщение об ошибке: объект типа «закрытие» не является подмножеством. > apoeht <- read.csv("apoeht.csv") > library(nlme) > model.a...

У меня есть следующие линейные модели library(nlme) fm2 <- lme(distance ~ age + Sex, data = Orthodont, random = ~ 1) fm2.lm <- lm(distance ~ age + Sex,data = Orthodont) Как получить стандартную ошибку расстояния с учетом возраста и...

В настоящее время я пытаюсь восстановить пример кода для lmeSplines (страницы 4 и 5 на официальная документация ). Первая строка кода (в синтаксисе R): # smoothing spline curve fit data(smSplineEx1) Как я могу просмотреть эти данные...

Я пытаюсь подогнать кривую логистического роста к некоторым данным с помощью функции gnls. Данные: structure(list(Nest = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 16L, 10L,...

В следующем примере я использую функцию-оболочку для размещения объекта gls. Я могу успешно вернуть прогноз модели для исходных данных, но не при использовании newdata , который возвращает ошибку « Error in eval(expr, envir, enclos) : 'nthcdr' needs...

У меня есть 100 групп по 40 наблюдений в каждой. Я знаю, что могут быть и другие подходящие модели, но в настоящее время меня интересует только следующее. reg<-lmList(Y ~ Intercept + a + b + c + d | grp,data=data, pool=F) Чтобы сохранить...

Я пытаюсь извлечь отдельные элементы (в частности, p-значения) из таблицы фиксированных эффектов, содержащейся в объекте, созданном сводным вызовом модели смешанных эффектов. Данные игрушек: set.seed(1234) score <- c(rnorm(8, 20, 3),...

Я хотел бы сделать модель 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)...

Я работал с Моделями смешанных эффектов в S и S-Plus Пиньеро и Бейтса и обнаружил много проблем, заставляющих работать модели из книги. Последним является использование функции постоянной плюс мощности для моделирования гетероскедастической...

Я использую R для анализа линейной модели смешанных эффектов (nlme::lme). Я хочу отобразить значения ACF пунктирной линией на альфа-уровне (как это делает plot.ACF). После того, как я require(nlme) , я могу запустить ACF, но plot.ACF по-прежнему...

Мне нужно ОЦЕНИТЬ ПРОДОЛЬНОЕ ИЗМЕНЕНИЕ ЖИРА МЕЖДУ 2 ПОСЕЩЕНИЯМИ с помощью линейной модели смешанных эффектов. У меня есть некоторые, которые будут меняться от визита 1 к визиту 2, так как это статус гипертонии, статус диабета, ИМТ, окружность...

Я пытаюсь построить график остаточных эффектов модели смешанных эффектов, используя функцию ggplot2. Однако после выполнения поиска я нашел некоторые доступные функции, но мне кажется, что для функции nlme они не работают. Графики, которые я...