R語言開發之協方差分析瞭解下
阿新 • • 發佈:2018-12-10
我們通常使用迴歸分析來建立描述預測變數變數對響應變數的影響的模型,有時,如果我們有類似於是/否或男/女等值的分類變數,簡單迴歸分析為分類變數的每個值提供多個結果。在這種情況下,我們可以通過使用分類變數和預測變數來研究分類變數的影響,並比較分類變數的每個級別的迴歸線。 這樣的分析被稱為協方差分析,也稱為ANCOVA。
我們接下來使用資料集mtcars
建立一個包含欄位“mpg”
,“hp”
和“am”
的資料框。 這裡我們將“mpg”
作為響應變數,將“hp”
作為預測變數,將“am”
作為分類變數,來看資料展示:
input <- mtcars[,c("am","mpg","hp")] print(head(input))
輸出結果如下:
am mpg hp
Mazda RX4 1 21.0 110
Mazda RX4 Wag 1 21.0 110
Datsun 710 1 22.8 93
Hornet 4 Drive 0 21.4 110
Hornet Sportabout 0 18.7 175
Valiant 0 18.1 105
我們再來嘗試建立一個迴歸模型,將“hp”
作為預測變數,將“mpg”
作為響應變數,我們需要考慮到“am”
和“hp”
之間的相互作用,來看模型與分類變數和預測變數之間的互動,如下:
# Get the dataset.
input <- mtcars
# Create the regression model.
result <- aov(mpg~hp*am,data = input)
print(summary(result))
輸出結果為:
Df Sum Sq Mean Sq F value Pr(>F) hp 1 678.4 678.4 77.391 1.50e-09 *** am 1 202.2 202.2 23.072 4.75e-05 *** hp:am 1 0.0 0.0 0.001 0.981 Residuals 28 245.4 8.8 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
該結果表明,馬力和變速箱型均具有對每加侖英里數的顯著影響,因為這兩種情況下的p
值均小於0.05
。但是,這兩個變數之間的相互作用並不重要,因為p
值大於0.05。
再來看分類變數與預測變數之間沒有互動的模型,如下:
# Get the dataset.
input <- mtcars
# Create the regression model.
result <- aov(mpg~hp+am,data = input)
print(summary(result))
輸出結果為:
Df Sum Sq Mean Sq F value Pr(>F)
hp 1 678.4 678.4 80.15 7.63e-10 ***
am 1 202.2 202.2 23.89 3.46e-05 ***
Residuals 29 245.4 8.5
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
該結果表明,馬力和變速箱型均具有對每加侖英里數的顯著影響,因為這兩種情況下的p
值均小於0.05。
最後我們就可以使用anova()
函式比較這兩個模型來確定變數的相互作用是否真的有意義的,如下:
# Get the dataset.
input <- mtcars
# Create the regression models.
result1 <- aov(mpg~hp*am,data = input)
result2 <- aov(mpg~hp+am,data = input)
# Compare the two models.
print(anova(result1,result2))
結果如下:
Model 1: mpg ~ hp * am
Model 2: mpg ~ hp + am
Res.Df RSS Df Sum of Sq F Pr(>F)
1 28 245.43
2 29 245.44 -1 -0.0052515 6e-04 0.9806
由於p
值大於0.05
,我們可以得出結論,馬力與傳播型別之間的相互作用不明顯。所以每加侖里程將取決於自動和手動變速模式的馬力。
好啦,本次記錄就到這裡了。
如果感覺不錯的話,請多多點贊支援哦。。。