Публикации по теме 'glm'
Обобщенная линейная модель (GLM)
Обобщенная линейная модель (GLM) — это статистическая модель, которая используется для прогнозирования переменной отклика на основе одной или нескольких переменных-предикторов. Это гибкий и широко используемый инструмент для моделирования данных, которые подчиняются определенному распределению, такому как биномиальное, пуассоновское или нормальное распределение.
В GLM переменная отклика связана с переменными-предикторами через линейный предиктор, который представляет собой линейную..
Основы машинного обучения (часть 5)
В этой части давайте немного обсудим множественную регрессию (ознакомьтесь с основами линейной модели в части 2), а затем сразу перейдем к общей линейной модели.
Множественная регрессия
Вычисление p-значения точно такое же, как и раньше, только на этот раз мы вносим изменения с помощью «pmean».
Мы также можем сравнить простую линейную регрессию с множественной регрессией и посмотреть, стоит ли вписываться в множественную регрессию с p-значением, как в следующем уравнении:
полученная выше ошибка исправлена следующим кодом
model <- glm(income ~. ,family=binomial(logit),data = train)
###Error### ---- Error in eval(family$initialize) : y values must be 0 <= y <= 1
полученная выше ошибка исправлена следующим кодом
Эта ошибка отслеживания может быть решена с помощью as.factor
model <- glm(as.factor(income) ~. ,family=binomial(logit),data = train)
###Error### ---- Error in eval(family$initialize) : y values must be 0 <= y <= 1
принимая в качестве фактора дохода ошибка очищается..
Пример биномиальной классификации в Scala и GLM с H2O
Вот пример задачи биномиальной классификации с использованием алгоритма H2O GLM с использованием набора данных кредитной карты на языке Scala.
Следующий пример предназначен для задачи полиномиальной классификации. Этот образец создан с использованием Spark 2.1.0 с Sparkling Water 2.1.4 .
import org.apache.spark.h2o._
import water.support.SparkContextSupport.addFiles
import org.apache.spark.SparkFiles
import java.io.File
import water.support.{H2OFrameSupport, SparkContextSupport,..
Вопросы по теме 'glm'
Как кратко написать формулу с множеством переменных из фрейма данных?
Предположим, у меня есть переменная ответа и данные, содержащие три ковариаты (в качестве игрушечного примера):
y = c(1,4,6)
d = data.frame(x1 = c(4,-1,3), x2 = c(3,9,8), x3 = c(4,-4,-2))
Я хочу подогнать к данным линейную регрессию:
fit...
135019 просмотров
schedule
30.05.2022
Почему каретный поезд занимает так много памяти?
Когда я тренируюсь только с использованием glm , все работает, и я даже близко не истощаю память. Но когда я запускаю train(..., method='glm') , мне не хватает памяти.
Это потому, что train хранит много данных для каждой итерации перекрестной...
9125 просмотров
schedule
20.07.2023
Как извлечь более фиксированные эффекты путем наблюдения?
У меня есть объект lme, созданный на основе данных о потреблении питательных веществ с повторными измерениями (два 24-часовых периода потребления на RespondentID):
Male.lme2 <- lmer(BoxCoxXY ~ -1 + AgeFactor + IntakeDay + (1|RespondentID),...
12119 просмотров
schedule
16.07.2022
Медленный расчет glm при использовании rpy2
Я хочу рассчитать параметры логистической регрессии с помощью пакета R glm. Я работаю с python и использую для этого rpy2. По какой-то причине, когда я запускаю функцию glm с помощью R, я получаю гораздо более быстрые результаты, чем при...
938 просмотров
schedule
07.08.2022
Векторизация glm и извлечение всей информации
У меня есть набор данных "ключевые слова" с несколькими группами. Я хочу применить glm к каждой группе отдельно, чтобы создать список подгонок glm с одной подгонкой для каждой группы.
Я мог бы сделать это с помощью цикла for , но это не в духе R....
396 просмотров
schedule
06.02.2023
Переключитесь на другую функцию при возникновении ошибки и сбоях программы в R
Возможное дублирование: Обработка исключений в R
Я использую функцию glm.nb для соответствия модели NB, но иногда программа выдает ошибку:
Error in while ((it <- it + 1) < limit && abs(del) > eps) { :
missing...
343 просмотров
schedule
12.01.2023
Как я могу выбрать модель с помощью AIC с помощью Gamma GLM в R?
Как объясняется в документации для glm (), компонент aic значения, возвращаемого glm (), не является допустимым AIC:
Для гауссовых, гамма- и обратных гауссовых семейств дисперсия оценивается по остаточному отклонению, а количество параметров...
2036 просмотров
schedule
28.02.2022
ошибка данных в bic.glm пакета BMA, но не в glm
Я оцениваю модель Пуассона по набору коэффициентов взаимодействия, а bic.glm пакета BMA помогает ориентироваться в пространстве модели. Я использую его в течение многих лет, но когда я вчера вечером обновил R с 2.10.x до 2.14.2, он перестал работать....
1553 просмотров
schedule
25.09.2022
Приступаем к работе с объектом rms lrm
Я не могу заставить логистическую регрессию из пакета rms работать с confint(), вот пример:
library(rms)
data(mtcars)
dd <- datadist(mtcars)
options(datadist = "dd")
fit <- lrm(am ~ gear + mpg, data=mtcars)
confint(fit)
Это дает ошибку:...
2889 просмотров
schedule
29.07.2023
glm и относительный риск - реплицировать код Stata в R
Может ли кто-нибудь помочь мне воспроизвести эти расчеты относительного риска (и его доверительный интервал) в R?
Аналогичная процедура, используемая в Stata, описана здесь . Может ли кто-нибудь показать мне, как это сделать в R (мои данные имеют...
2847 просмотров
schedule
15.09.2023
Подбор модели: glm vs glmmPQL
Я подбираю модель, касающуюся данных об отсутствии-присутствии, и хотел бы проверить, является ли случайный фактор значимым или нет. Для этого нужно сравнить glmm с glm и проверить с помощью LR-теста, какой из них наиболее значимый, если я правильно...
3385 просмотров
schedule
20.09.2022
Выбор статистически значимых переменных в модели R glm
У меня есть переменная результата, скажем Y, и список из 100 измерений, которые могут повлиять на Y (скажем, X1 ... X100).
После запуска моего glm и просмотра сводки моей модели я вижу статистически значимые переменные. Я хотел бы иметь...
22515 просмотров
schedule
24.06.2023
На сутки вперед с использованием модели GLM в R
У меня есть следующий код, чтобы получить прогноз потребления нагрузки на сутки вперед с 15-минутным интервалом с использованием температуры наружного воздуха и TOD (96 категориальных переменных, время дня). Когда я запускаю приведенный ниже код, я...
2546 просмотров
schedule
04.09.2022
прогнозировать и мультипликативные переменные / условия взаимодействия в пробит-регрессиях
Я хочу определить предельные эффекты каждой зависимой переменной в пробит-регрессии следующим образом:
предсказать (базовую) вероятность со средним значением каждой переменной
для каждой переменной предсказать изменение вероятности по сравнению...
3148 просмотров
schedule
16.09.2022
отрицательная биномиальная функция в r
Я сравниваю заболеваемость в четырех разных регионах, чтобы увидеть, равны они или нет. Сначала я подогнал модель Пуассона, но я не уверен, как я могу проверить предположение о дисперсии = среднему, когда я использую смещение. Следовательно, я хотел...
1271 просмотров
schedule
03.07.2022
Ошибка R, которая говорит, что не все модели были приспособлены к одному и тому же размеру набора данных.
Я создал две обобщенные линейные модели следующим образом:
glm1 <-glm(Y ~ X1 + X2 + X3, family=binomial(link=logit))
glm2 <-glm(Y ~ X1 + X2, family=binomial(link=logit))
Затем я использую функцию anova :
anova(glm2,glm1)
но...
34772 просмотров
schedule
15.06.2022
Как увеличить количество перестановок в функции перекрестной проверки GLM `cv.glm`
Мне интересно, есть ли у кого-нибудь опыт указания количества перестановок в функции перекрестной проверки GLM cv.glm (пакет: загрузка)? Я использую настройки K=2 для разделения наборов для обучения и проверки примерно на две равные группы. Мои...
379 просмотров
schedule
14.06.2022
R Пакет сегментирован со смещением
Я использую пакет сегментированный в R. Сначала я установил пуассоновский glm с функцией ссылки журнала, использующей функцию glm , которая включает термин смещения, например:
M1=glm(Y~X1+X2+X3+offset(log(X)),data=dat.1,family=poisson)
M1...
379 просмотров
schedule
13.06.2023
Извлечь остаточное отклонение из анова (glm) в R
Я установил модель glm в R и взял таблицу Anova. Мне нужно извлечь столбец «Остаточное отклонение». Но выдает ошибку. Вот коды:
Создание данных:
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)...
3676 просмотров
schedule
08.07.2023
Почему функция R `glm` не работает правильно, когда строковый аргумент передается для формулы одновременно с указанием подмножества?
Первые два вызова glm работают правильно, а последний выдает ошибку.
Почему это так, и есть ли какие-либо известные обходные пути, если я хочу использовать строку для формулы, а также указать подмножество? Вариант использования заключается в...
2568 просмотров
schedule
30.07.2023
