Публикации по теме '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 просмотров
schedule
02.12.2022
Как рассчитать прогнозируемую вероятность модели отрицательной биномиальной регрессии?
Я использую функцию glm.nb() в пакете R MASS для оценки параметров модели отрицательной биномиальной регрессии. Как я могу рассчитать прогнозируемую вероятность (функцию массы вероятности) с учетом новых данных, какую функцию R я могу...
6000 просмотров
schedule
10.03.2023
Как подогнать авторегрессионную смешанную модель Пуассона (подсчет временных рядов) в R?
Моя задача состоит в том, чтобы оценить, как различные переменные среды влияют на годовые колебания численности населения. Для этого мне нужно подобрать авторегрессионную модель Пуассона для подсчета временных рядов:
Где N i,j — количество...
4986 просмотров
schedule
04.06.2023
Ошибка допустимых начальных значений для логистической регрессии GLM
Я относительно новичок в R и пытаюсь анализировать данные об успешном вложении, используя обновленную версию кода функции связи логистической экспозиции, созданную Shaffer 2004 и предоставленную в семействе файлов справки R.
Вот код, который я...
1551 просмотров
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