如何使用R软件查找两个模型的AIC值？

Question

我正在研究生存分析。

我估计了Cox回归模型和Buckley＆James回归模型。

为了确定哪个模型更适合我的数据集，我使用了Akaike信息标准（AIC）。那么，如何使用R软件找到两个模型的AIC值？

Answer 1

如果您正在寻找AIC值，可以使用glm函数找到它们并将其保存为向量x。然后执行摘要（x），您将看到所有AIC，BIC等。以下是使用mtcars数据集的示例

> data(mtcars) #loads data
> head(mtcars) #summary view of data
                   mpg cyl disp  hp drat    wt  qsec vs am gear carb
Mazda RX4         21.0   6  160 110 3.90 2.620 16.46  0  1    4    4
Mazda RX4 Wag     21.0   6  160 110 3.90 2.875 17.02  0  1    4    4
Datsun 710        22.8   4  108  93 3.85 2.320 18.61  1  1    4    1
Hornet 4 Drive    21.4   6  258 110 3.08 3.215 19.44  1  0    3    1
Hornet Sportabout 18.7   8  360 175 3.15 3.440 17.02  0  0    3    2
Valiant           18.1   6  225 105 2.76 3.460 20.22  1  0    3    1

> x<-glm(mtcars$cyl~mtcars$mpg) #creates a regression model
> summary(x) #summary of regression model

Call:
glm(formula = mtcars$cyl ~ mtcars$mpg)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.8569  -0.6484   0.1205   0.5965   1.5876  

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) 11.26068    0.59304   18.99  < 2e-16 ***
mtcars$mpg  -0.25251    0.02831   -8.92 6.11e-10 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for gaussian family taken to be 0.9024651)

    Null deviance: 98.875  on 31  degrees of freedom
Residual deviance: 27.074  on 30  degrees of freedom
AIC: 91.463                                      #AIC value you are looking for

Number of Fisher Scoring iterations: 2