По вопросам, связанным с обобщенными линейными моделями. Для библиотеки математики GLM, см. [Tag: glm-math].

Подробнее про glm...

Я пытаюсь выполнить анализ логистической регрессии для каждой строки от R1 до R6 по строке C1 кадра данных "dat" Input data frame dat <- data.frame(list(X1 = c(1, 1, 0, 1, 1, 1, 0), X2 = c(1, 1, 0, 1, 1, 1, 0), X3 = c(1, 1, 1, 0, 1, 1, 0), X4 = c(1, 1, 0, 1, 1, 1, 0), X5 = c(0, 1, 0, 0, 1, 1, 0), X....
9 Окт 2021 в 17:48
Я пытаюсь выяснить, как коэффициенты логистической регрессии с полиномиальным членом связаны с прогнозами. В частности, меня интересует место на оси x, где прогноз наиболее высок. Пример ниже: set.seed(42) # Setup some dummy data x <- 1:200 y <- rep(0, length(x)) y[51:150] <- rbinom(100, 1, 0.5) ....
3 Авг 2021 в 16:43
Я пытаюсь создать определяемую пользователем функцию, которая выводит отношения шансов (95% ДИ) из моделей логистической регрессии. Я в основном застрял на Строке 2, где я предсказываю отношение шансов с помощью функции oddsratio::or_glm(). Мне нужно передать ввод предиктора в именованный список в ....
Следующий код df <- data.frame(place = c("South","South","North","East"), temperature = c(30,30,20,12), outlookfine=c(TRUE,TRUE,FALSE,FALSE) ) glm.fit <- glm(outlookfine ~ .,df , family= binomial) coef.glm <-coef(summary(glm.fit)) coef.glm Выходы ....
21 Май 2021 в 01:38
В следующем примере df <- data.frame(place = c("South","South","North"), temperature = c(30,30,20), outlookfine=c(TRUE,TRUE,FALSE) ) glm.fit <- glm(outlookfine ~ .,df, family=binomial) glm.fit На выходе Call: glm(formula = outlookfine ~ ., family....
r glm
21 Май 2021 в 01:24
Я пытаюсь настроить обобщенную линейную модель, определенную ниже: Данные: https://drive.google.com/file/d/ 1iJdZbEVwf0_zFnGUrNmVBCkmpzgI_fr1 / view? Usp = sharing Следует отметить, что переменная ответа Var1, а также переменная регрессора Var2 имеют нулевые значения, для которых была добавлена ​​к....
12 Май 2021 в 16:00
Я пытаюсь использовать edgeR для анализа дифференциальной экспрессии набора данных биологического подсчета. Мои образцы разделены на случайные и контрольные, и я хотел бы знать, какие гены регулируются в большей или меньшей степени в образцах клинических случаев (т. Е. Тех, которые имеют заболевани....
9 Май 2021 в 19:36
Поэтому я хочу найти оценочный параметр с помощью GLM и сравнить его с пакетом mle2. Вот мой код для GLM d <- read.delim("http://dnett.github.io/S510/Disease.txt") d$disease=factor(d$disease) d$ses=factor(d$ses) d$sector=factor(d$sector) str(d) glm2 <- glm(disease~ses+sector, family=binomial(link=....
26 Апр 2021 в 10:23
Я использую взвешенные обобщенные линейные модели (statsmodels) для классификации: import statsmodels.api as sm model= sm.GLM(y, x_with_intercept, max_iter=500, random_state=42, family=sm.families.Binomial(),freq_weights=weights) Одна из переменных в x_with_intercept - двоичная. Я думал, что не во....
24 Апр 2021 в 09:46
Я новичок в инструменте R, и у меня возникли проблемы с функцией glm(). У меня есть некоторые данные, которые я показал ниже. Когда линейный предиктор равен x, функция glm() работает нормально, но как только я меняю линейный предиктор на x + x ^ 2, она начинает давать мне те же результаты, что и д....
r glm
23 Апр 2021 в 01:12
Рассмотрим список аннотаций: set.seed(42) bin_var <- sample(0:1, 125, T) indep_1 <- rnorm(125) indep_2 <- rexp(125) indep_3 <- runif(125) summary_list <- list(summary(glm(bin_var~indep_1)), summary(glm(bin_var~indep_2)), summary(glm(bin_var~indep_3))) В и....
16 Апр 2021 в 11:48
Кросс размещен на CrossValidated . Я пытаюсь загрузить свои результаты для варианта подхода 2SLS (2SRI) на основе эту ссылку. По какой-то причине бутстрап не дает никаких результатов. library(sure) # for residual function and sample data sets library(MASS) # for polr function rm(df1) df1 <- df1 df....
9 Апр 2021 в 09:59
Чтобы проверить, существует ли связь между группами болезней (категориальная_переменная) и заболеванием (исход; дихотомия), я провожу логистическую регрессию. Чтобы проверить наличие смешения со стороны других переменных, я запускаю 3 модели с различным количеством ковариат. Для отображения OR и CI....
7 Апр 2021 в 17:22
Моя цель - найти параметры оценки с помощью пакета optim () в R. И я сравниваю свой результат с моделью GLM в R. d <- read.delim("http://dnett.github.io/S510/Disease.txt") d$disease=factor(d$disease) d$ses=factor(d$ses) d$sector=factor(d$sector) str(d) oreduced <- glm(disease~age+sector, family=bin....
4 Апр 2021 в 12:53
Я хочу попытаться найти оценку максимального правдоподобия с помощью пакета optim () на R. Я начинаю с моделирования данных с помощью модели GLM, чтобы сравнить оценку с optim (). Вот мой код d....
3 Апр 2021 в 13:35
Это моя первая публикация здесь, поэтому заранее извиняюсь, если есть какие-либо проблемы с моим вопросом или воспроизводимым примером. Я проверил похожие вопросы и ответы на них, но не смог .......
31 Мар 2021 в 19:08
Я работаю с набором данных orings в пакете faraway в R. Я написал следующую сгруппированную биномиальную модель: orings_model....
31 Мар 2021 в 17:08
Я не уверен, что то, что я пытаюсь сделать, возможно, поэтому я подумал, что посоветуюсь с более опытными пользователями SAS. Я хочу запускать Proc GLM много-много раз, используя оператор «by», а затем иметь выходные параметры для каждого запуска в одном файле. Выходные параметры - это что-то вроде....
30 Мар 2021 в 22:29
Похоже, что эталонный уровень выбран как первый уникальный элемент категориального значения. Однако в моем случае это P (а не A). ols_all= lm( @formula( value ~ Treatment ), s_treat) Дает value ~ 1 + Treatment Coefficients: ──────────────────────────────────────────────────────────────────────── ....
23 Мар 2021 в 17:02
Оценить точность модели с помощью результатов Бернулли довольно просто, но я не уверен, как сгенерировать значимые прогнозы на основе агрегированной биномиальной регрессии. Возьмем этот пример. Мы хотим моделировать .......
18 Мар 2021 в 23:17
Я пытаюсь подогнать модель GLM в Python 3.1. Я использую функцию связи Гаусса. Поэтому я написал import statsmodels.api как sm glm = sm (y_train, X_train, family = sm.families.Gaussian ()) Но я получаю сообщение об ошибке .......
11 Мар 2021 в 17:03
Я получаю постоянную ошибку «Ошибка: переменные 'x1', 'x2', 'x3' были указаны с разными типами, отличными от соответствия» при попытке предсказать результат новых данных, используя прогноз в R. Я запустил это код с успехом на других моделях, но по какой-то причине я не могу понять, что не так с это....
8 Мар 2021 в 20:44
Я пытаюсь вставить в график значения отношения шансов. Я также пытаюсь придать фигурам другой цвет со значительным отношением шансов. Значимое отношение шансов - это те, доверительные интервалы которых .......
25 Фев 2021 в 14:58
Я пытаюсь получить информацию из результатов OLS после регрессии цикла. Например, depvars = ['y1', 'y2', 'y3', ...] models = [ "~ x1 + x2", "~ x1 + x2 + x3", ...] results = [] for depvar in depvars: for model in models: results.append(smf.glm(formula = depvar + model, data= data).fit()) ....
25 Фев 2021 в 14:15
Я запустил регресс примерно так, import statsmodels.api as sm import statsmodels.formula.api as smf from statsmodels.genmod.generalized_linear_model import GLMResults result = smf.glm(formula = 'y ~ x1 + x2', data = data).fit() Я мог получить оценки, p-значения и количество наблюдений result.para....
25 Фев 2021 в 11:17