Добавить столбец в вывод регрессии - PullRequest
1 голос
/ 13 апреля 2020

Я хотел бы добавить столбец с X² (хи-квадрат), а также столбец с Exp (B) в мой вывод регрессии. Есть ли идея, как это сделать? Спасибо заранее. На данный момент я рассчитал это вручную для каждой модели и переменной, что довольно много времени.

model_simple <- as.formula("completion_yesno ~  ac + ov + UCRate + FirstWeek + LastWeek + DayofWeekSu + DayofWeekMo + DayofWeekTu + DayofWeekWe + DayofWeekTh + DayofWeekFr + MonthofYearJan + MonthofYearFeb + MonthofYearMar + MonthofYearApr +MonthofYearMay+ MonthofYearJun + MonthofYearJul + MonthofYearAug + MonthofYearSep + MonthofYearOct + MonthofYearNov")
clog_simple1 = glm(model_simple,data=cllw,family = binomial(link = cloglog))
summary(clog_simple1)

1 Ответ

1 голос
/ 13 апреля 2020

Может быть, вы можете уточнить, что вы подразумеваете под хи-квадрат и exp (B). Вы можете сделать следующее:

da <- MASS::Pima.tr
model <- glm(type ~ .,data=da,family = binomial(link = cloglog))
results <- data.frame(coefficients(summary(model)),check.names=FALSE)
# some random values
results$chisq = rchisq(nrow(results),1)
results$expB = exp(results$Estimate)

Или вы можете использовать приборку из метлы:

library(broom)
results = tidy(model)
results$expB = exp(results$Estimate)
...