使用R软体进行双样本t检定,.PDF

基本統計 15 使用 R 軟體進行雙樣本 t 檢定， 比較男女服務員的小費差異 t 檢定常常被用來比較兩組樣本。我們使用 tips(小費)資料來說明此檢定， 比較男性和女性服務員所收到的小費是否有差別。在進行 t 檢定之前，我們 必須檢測兩組樣本的變異數。傳統的 t 檢定要求兩組的變異數必須一樣，而 Welch 雙樣本 t 檢定則可以處理變異數不同的樣本群。我們從數值上和圖表上 來探討這些檢定，如圖 15.5 顯示。 圖 15.5 以性別來做分類的小費額直方圖。可以注意到兩個分佈都沒有呈現常態。 data (tips, package = reshape2] # 首先對各組計算變異數; # 使用 formula 介面 # 計算每個性別小費的變異數 aggregate(tip ~ sex, data=tips, var) sex tip 1 Female 1.3444282 2 Male 2.217424 # 現在檢測小費的分佈是否為常態 shapiro.test(tips$tip) Shapiro-Wilk normality test data: tips$tip W = 0.8978, p-value = 8.2e-12 shapiro.test(tips$tip[tips$sex == Female]) Shapiro-Wilk normality test data: tips$tip[tips$sex == Female] W = 0.9568, p-value = 0.005448 shapiro.test(tips$tip[tips$sex == Male]) Shapiro-Wilk normality test data: tips$tip[tips$sex == Male] W = 0.8759, p-value = 3.708e-10 # 用目測可以判斷所有檢測都不通過 ggplot(tips, aes(x=tip, fill=sex)) + + geom histogram(binwidth=.5, alpha=1/2) 由於資料都不呈現常態分佈，標準的 F 檢定(通過 var.test 函數)和 Bartlett 檢定(通過 bartlett.test 函數)都不能被使用。因此，我們用無母數 Ansari-Bradley 檢定來檢測變異數是否相等。 基本統計 15 ansari.test(tip ~ sex, tips) Ansari-Bradley test data: tip by sex AB = 5582.5, p-value = 0.376 alternative hypothesis: true ratio of scales is not equal to 1 檢定結果顯示兩組樣本的變異數相等，因此我們可以用最標準的雙樣本 t 檢定。 # 設定 var.equal=TRUE 將進行標準的雙樣本 t 檢定 # 設定 var.equal=FALSE(預設)則進行 Welch 檢定 t.test(tip ~ sex, data = tips, var.equal = TRUE) Two Sample t-test data: tip by sex t = -1.3879, df = 242, p-value = 0.1665 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval: -0.6197558 0.1074167 sample estimates: mean in group Female mean in group Male 2.833448 3.089618 根據此檢定，結果為不顯著的，因此結論為男性和女性服務員所收到的小費 額是差不多一樣的。雖然用精確的統計方法進行檢