[R]{stats}influence.measures() : 外れ値削除法による回帰診断を実行する関数

説明

influence.measures()は、線形回帰分析モデルおよび一般化線形回帰分析モデルに対して、外れ値削除法による回帰診断を実行する関数である。Belsley, Kuh and Welsch (1980), Cook and Weisberg (1982)などの書籍で詳細は論じられている。

使用法

influence.measures(model)

引数

model ：１つのRオブジェクト、典型的にはlm()やglm()の返り値オブジェクト。

実行例

Rヘルプのサンプルスクリプトに実行結果をつけたものを掲げる。

> lm.SR <- lm(sr ~ pop15 + pop75 + dpi + ddpi, data = LifeCycleSavings)
> inflm.SR <- influence.measures(lm.SR)
# モデルを歪ませている可能性のある外れ値のみを挙げる
> which(apply(inflm.SR$is.inf, 1, any))
        Chile United States        Zambia         Libya 
            7            44            46            49 
> summary(inflm.SR) # only these。summary()を介さないと全データが出力される
Potentially influential observations of
	 lm(formula = sr ~ pop15 + pop75 + dpi + ddpi, data = LifeCycleSavings) :

              dfb.1_ dfb.pp15 dfb.pp75 dfb.dpi dfb.ddpi dffit   cov.r   cook.d hat    
Chile         -0.20   0.13     0.22    -0.02    0.12    -0.46    0.65_*  0.04   0.04  
United States  0.07  -0.07     0.04    -0.23   -0.03    -0.25    1.66_*  0.01   0.33_*
Zambia         0.16  -0.08    -0.34     0.09    0.23     0.75    0.51_*  0.10   0.06  
Libya          0.55  -0.48    -0.38    -0.02   -1.02_*  -1.16_*  2.09_*  0.27   0.53_*
>