Публикации по теме 'glm'


Установка GLM вручную (в R)
Команда glm в R (или функция GLM из библиотеки statsmodel.api в Python) может соответствовать модели GLM без необходимости знать все сложные уравнения и методы оптимизации, используемые «за кулисами». Но умение что-то кодировать вручную обычно является признаком глубокого понимания. Итак, в этой статье мы собираемся установить GLM вручную в R и сравнить его со встроенным glm . Уравнения вероятности Уравнения правдоподобия относятся к первой производной логарифмического..

Вопросы по теме 'glm'

Предсказание / вменение недостающих значений регрессии Пуассона GLM в R?
Я пытаюсь изучить способы вменения недостающих значений в набор данных. Мой набор данных содержит количество случаев (Неестественное, Естественное и Сумма Всего) для Года (2001-2009), Месяца (1-12), Пола (М / Ж) и Возрастной группы (4 группы)....
3359 просмотров
schedule 07.10.2022

Интерпретация имен коэффициентов в glmnet в R
Я использую glmnet для прогнозирования вероятностей на основе набора из 5 функций, используя следующий код. Мне нужна фактическая формула, потому что мне нужно использовать ее в другой (не R) программе. deg = 3 glmnet.fit <-...
8719 просмотров
schedule 24.11.2023

Разница между svyglm и взвешенным glm
Я хотел бы знать, в чем разница между использованием svyglm или взвешенного glm . Например: M1 = glm(formula = yy ~ age + gender + country , family = binomial(link = "probit"), data = P2013, subset =...
3908 просмотров
schedule 02.03.2023

Доступ к fit.values ​​при использовании ddply
Я использую ddply для выполнения glm на подмножествах моих данных. У меня возникают трудности с доступом к оценочным значениям Y. Я могу получить оценки параметров модели, используя приведенный ниже код, но все варианты, которые я пытался...
560 просмотров
schedule 17.08.2022

бенчмаркинг логистической регрессии с использованием glm.fit, bigglm, speedglm, glmnet, LiblineaR
Я имитирую данные и сравниваю glm.fit , bigglm, speedglm, glmnet, LiblineaR для бинарной логит-модели. testGLMResults_and_speed <- function(N, p, chunk=NULL, Ifsample=FALSE, size=NULL, reps=5){ library(LiblineaR) library(speedglm)...
3630 просмотров
schedule 15.06.2023

Ошибка при запуске glm с family=gamma
Я обычно использую SAS, поэтому я не слишком хорошо знаком с R, поэтому извините, если это основной вопрос. Я запустил модель, и она выдает следующую ошибку Error in family() : 0 arguments passed to 'gamma' which requires 1 Кто-нибудь знает,...
4498 просмотров
r glm
schedule 02.12.2022

Как рассчитать прогнозируемую вероятность модели отрицательной биномиальной регрессии?
Я использую функцию glm.nb() в пакете R MASS для оценки параметров модели отрицательной биномиальной регрессии. Как я могу рассчитать прогнозируемую вероятность (функцию массы вероятности) с учетом новых данных, какую функцию R я могу...
6000 просмотров
r glm
schedule 10.03.2023

Как подогнать авторегрессионную смешанную модель Пуассона (подсчет временных рядов) в R?
Моя задача состоит в том, чтобы оценить, как различные переменные среды влияют на годовые колебания численности населения. Для этого мне нужно подобрать авторегрессионную модель Пуассона для подсчета временных рядов: Где N i,j — количество...
4986 просмотров
schedule 04.06.2023

Ошибка допустимых начальных значений для логистической регрессии GLM
Я относительно новичок в R и пытаюсь анализировать данные об успешном вложении, используя обновленную версию кода функции связи логистической экспозиции, созданную Shaffer 2004 и предоставленную в семействе файлов справки R. Вот код, который я...
1551 просмотров
r glm
schedule 24.09.2022

Решающий порог для модели логистической регрессии glm в R
У меня есть некоторые данные с предикторами и двоичной целью. Например: df <- data.frame(a=sort(sample(1:100,30)), b= sort(sample(1:100,30)), target=c(rep(0,11),rep(1,4),rep(0,4),rep(1,11))) Я обучил модель логистической...
31661 просмотров
schedule 22.02.2024

Пользовательская стандартная ошибка начальной загрузки: числовой аргумент envir не равен единице
Я пишу собственный сценарий для загрузки стандартных ошибок в GLM в R и получаю следующую ошибку: Ошибка в eval (predvars, data, env): числовой аргумент envir не равен единице Может кто-нибудь объяснить, что я делаю не так? Мой код:...
2268 просмотров
schedule 26.05.2022

Байес предсказывает, нижний индекс выходит за пределы
У меня возникли проблемы с функцией прогнозирования при использовании байесовской программы. Я читал некоторые сообщения, в которых говорится, что эта проблема может возникнуть, когда данные вне выборки имеют больше уровней, чем данные в выборке, но...
3233 просмотров
schedule 28.04.2022

Модель работает с glm, но не с bigglm
Я пытался запустить логистическую регрессию на 320 000 строк данных (6 переменных). Пошаговый выбор модели на выборке данных (10000) дает довольно сложную модель с 5 элементами взаимодействия: Y~X1+ X2*X3+ X2*X4+ X2*X5+ X3*X6+ X4*X5 . Функция...
2959 просмотров
schedule 19.01.2023

Как рассчитать срок взаимодействия как отношение шансов в логистической регрессии?
Я использую логистическую регрессию с фиксированными эффектами в R, используя функцию glm . Я прочитал кое-что об интерпретации терминов взаимодействия в обобщенных линейных моделях. При использовании логарифма шансов модель является линейной, и...
5086 просмотров
schedule 28.06.2023

Полиномиальная модель MCMCglmm в R
Я пытаюсь создать модель с помощью пакета MCMCglmm в R. Данные структурированы следующим образом, где диада, фокус, прочее - все случайные эффекты, предикторы 1-2 - переменные-предикторы, а ответы 1-5 - переменные результата, которые фиксируют...
4511 просмотров
schedule 22.03.2022

Как я могу создать предельные эффекты для логит-модели при использовании весов обследования?
Обычно я генерирую предельные эффекты модели logit, используя пакет mfx и функцию logitmfx. Однако текущий опрос, который я использую, имеет веса (которые имеют большое влияние на долю DV в выборке из-за передискретизации в некоторых популяциях), и...
2329 просмотров
schedule 23.03.2022

Ошибка при выполнении пошаговой регрессии
Я пытаюсь запустить модель пошаговой регрессии. Я продолжаю получать это сообщение: #Error in step(cdc.fit, direction = "backward") : # number of rows in use has changed: remove missing values? #In addition: There were 50 or more warnings (use...
842 просмотров
schedule 16.06.2023

Определите точки на остаточных графиках GLM в R
Как я могу идентифицировать точки на остаточных графиках, которые я получаю, запуская plot(glm.object) в R? Я понимаю, как работает identify() , но где я могу получить значения x и y, необходимые для этих графиков? Я также знаю о...
1338 просмотров
schedule 19.12.2022

Расчет предельных эффектов для взвешенной логит-модели
Я рассчитал взвешенную логит-модель с помощью пакета «Survey». library(survey) Mod_design <- svrepdesign(variables = Data, weights = weights, repweights = REP_WGT, type = "JKn",...
420 просмотров
schedule 17.07.2022

Вывод обобщенной линейной модели через texreg
Я могу использовать texreg , чтобы получить красивый вывод glm, который будет использоваться для knitr . Иногда нам нужно преобразовать вывод glm обратно в ответ, используя обратную ссылку. Интересно, как получить обратную ссылку с texreg ....
915 просмотров
schedule 03.01.2022