Вопросы с тегом «lm»

lm - это имя функции линейной модели (т. е. множественной регрессии) в статистическом пакете R. Для линейных моделей, как правило, используйте вместо него тег linear-model.

115
Как стандартные ошибки коэффициентов рассчитываются в регрессии?

Для моего собственного понимания я заинтересован в том, чтобы вручную повторить вычисление стандартных ошибок оценочных коэффициентов, поскольку, например, они поставляются с выходными данными lm()функции R, но не смогли ее определить. Какая формула / реализация...

37
Как найти подходящую для полусинусоидальной модели модель R?

Я хочу предположить, что температура поверхности моря в Балтийском море один и тот же год за годом, а затем описать это с помощью функции / линейной модели. У меня была идея просто ввести год в виде десятичного числа (или num_months / 12) и узнать, какой должна быть температура в это время. Бросив...

35
Что такое скорректированная формула R-квадрата в lm в R и как ее следует интерпретировать?

Какая точная формула используется в R lm() для Скорректированного R-квадрата? Как я могу интерпретировать это? Скорректированные R-квадрат формулы Кажется, существует несколько формул для расчета скорректированного R-квадрата. Формула...

22
Пошаговая линейная алгебраическая вычисление регрессии наименьших квадратов

В качестве приквела к вопросу о линейно-смешанных моделях в R и для использования в качестве справочного материала для любителей статистики начального и среднего уровня, я решил опубликовать в качестве независимого «Q & A-style» шаги, связанные с «ручным» вычислением Коэффициенты и прогнозные...

19
Лучший способ справиться с гетероскедастичностью?

У меня есть график остаточных значений линейной модели в зависимости от подогнанных значений, где гетероскедастичность очень ясна. Однако я не уверен, как мне поступить сейчас, потому что, насколько я понимаю, эта гетероскедастичность делает мою линейную модель недействительной. (Это правильно?)...

19
Является ли взвешенный

Я оценил надежную линейную модель Rс весами ММ, используя rlm()пакет MASS. `R`` не предоставляет значение для модели, но я хотел бы иметь его, если это значимое количество. Мне также интересно знать, есть ли смысл иметь значение которое взвешивает общую и остаточную дисперсию так же, как...

14
Использование R для GLM с гамма-распределением

В настоящее время у меня проблема с пониманием синтаксиса R для подгонки GLM с использованием гамма-распределения. У меня есть набор данных, где каждая строка содержит 3 ко-вариации ( ), переменную ответа ( ) и параметр формы ( K ). Я хочу смоделировать масштаб гамма-распределения как линейную...

13
R: проверить нормальность остатков линейной модели - какие остатки использовать

Я хотел бы сделать W-тест Шапиро Уилка и тест Колмогорова-Смирнова на невязках линейной модели, чтобы проверить на нормальность. Мне было просто интересно, какие остатки следует использовать для этого - необработанные остатки, остатки Пирсона, студентизированные остатки или стандартизированные...

13
Почему R's lm () возвращает оценки коэффициентов, отличные от моего учебника?

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

11
Оценка параметров динамической линейной модели

Я хочу реализовать (в R) следующую очень простую динамическую линейную модель, для которой у меня есть 2 неизвестных изменяющихся во времени параметра (дисперсия ошибки наблюдения и дисперсия ошибки состояния ). ε 2 тε1Tϵt1\epsilon^1_tε2Tϵt2\epsilon^2_t YTθт + 1знак равнознак равноθT+ ϵ1TθT+...

10
Робастный регрессионный вывод и сэндвич-оценки

Можете ли вы дать мне пример использования сэндвич-оценок для выполнения надежного регрессионного вывода? Я могу видеть пример ?sandwich, но я не совсем понимаю, как мы можем перейти от lm(a ~ b, data)( r- кодированного) к оценке и к значению p, полученному в результате регрессионной модели с...

10
Повторные измерения anova: lm vs lmer

Я пытаюсь воспроизвести несколько тестов взаимодействия между обоими lmи lmerна повторных измерениях (2x2x2). Причина, по которой я хочу сравнить оба метода, состоит в том, что GLM SPSS для повторных измерений дает те же самые результаты, что и lmподход, представленный здесь, поэтому в конце я хочу...

10
Пересчитать логарифмическое правдоподобие из простой модели R lm

Я просто пытаюсь пересчитать с помощью dnorm () логарифмическую вероятность, обеспечиваемую функцией logLik из модели lm (в R). Это работает (почти идеально) для большого количества данных (например, n = 1000): > n <- 1000 > x <- 1:n > set.seed(1) > y <- 10 + 2*x + rnorm(n, 0,...

10
Почему Anova () и drop1 () предоставили разные ответы для GLMM?

У меня есть GLMM формы: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Когда я использую drop1(model, test="Chi"), я получаю другие результаты, чем если бы я использовал Anova(model, type="III")из пакета автомобиля или summary(model). Последние...

9
Установка изменяющегося во времени коэффициента DLM

Я хочу приспособить DLM с изменяющимися во времени коэффициентами, то есть расширением к обычной линейной регрессии, .YT= θ1+ θ2Икс2YTзнак равноθ1+θ2Икс2y_t = \theta_1 + \theta_2x_2 У меня есть предиктор ( ) и переменная отклика ( y t ), ежегодный морской и внутренний вылов рыбы соответственно с...

9
Почему линейная регрессия не способна предсказать исход простой детерминированной последовательности?

Мой коллега прислал мне эту проблему, очевидно, делая обходы в Интернете: If $3 = 18, 4 = 32, 5 = 50, 6 = 72, 7 = 98$, Then, $10 =$ ? Ответ, кажется, 200. 3*6 4*8 5*10 6*12 7*14 8*16 9*18 10*20=200 Когда я делаю линейную регрессию в R: data <- data.frame(a=c(3,4,5,6,7), b=c(18,32,50,72,98)) lm1...

9
Подгонка гетероскедастической обобщенной линейной модели для биномиальных ответов

У меня есть данные из следующего экспериментального плана: мои наблюдения - это подсчет числа успехов ( K) из соответствующего числа испытаний ( N), измеренных для двух групп, каждая из которых состоит из Iиндивидуумов, из Tобработок, где в каждой такой комбинации факторов есть Rповторения , Таким...

9
Параметрический, полупараметрический и непараметрический бутстрап для смешанных моделей

Следующие прививки взяты из этой статьи . Я новичок в начальной загрузке и пытаюсь реализовать параметрическую, полупараметрическую и непараметрическую загрузку начальной загрузки для линейной смешанной модели с R bootпакетом. Код R Вот мой Rкод: library(SASmixed) library(lme4) library(boot)...