如果我只有模型拟合调用和摘要，有没有办法重用glm拟合？

Question

我正在处理一个由前程序员构建的项目。程序员在其文档中对

glm()

的调用是

glm(formula = AVAL ~ AUC, family = binomial(), data = logreg.dat)

我有该代码的输出，但没有用于构建模型的数据。

Deviance Residuals:
Min 1Q Median 3Q Max
-0.5991 -0.4784 -0.4220 -0.3028 2.7425

Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -1.62668 0.34738 -4.683 2.83e-06 ***
AUC -0.14730 0.07448 -1.978 0.048 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

Null deviance: 119.87 on 193 degrees of freedom
Residual deviance: 113.23 on 192 degrees of freedom
AIC: 117.23

Number of Fisher Scoring iterations: 6

有办法重建这种契合吗？我想使用拟合来模拟一些基于新 AUC 向量的新 AVAL 数据。

Answer 1

如果我理解正确，您想使用参数拟合的估计来模拟模型。你可以这样做：

# take the estimates from the output
intercept <- -1.62668
slope <- -0.14730 

# new AUC values
AUC <- c(1, 2, 3, 4, 5)

# simulate AVAL values
logit_probs <- intercept + slope * AUC
probs <- 1 / (1 + exp(-logit_probs))
n <- length(AUC)
rbinom(n, 1, probs)