本文章最後由 LouisBuLu 於 2014-6-25 01:08 編輯
看完數值分析部分後，我建議使用變異數分析(ANOVA)進行統計
建議最大值和最小值都不要刪掉，以增加樣本選取之個數
因為要進行多方比較，使用ANOVA會較佳
虛無假設: 六組火星塞輸出馬力無差異
對立假設: 六組火星塞輸出馬力有差異
在變異數分析後，F值若大於臨界值，
表示六組是有差異的
接下來可再進行後續檢定(例如幫佛諾利T檢定)，
以釐清哪幾組馬力輸出無差異、哪幾組馬力輸出是不同的
因為我電腦無法安裝馬力輸出的監看軟體，
以下僅能暫時以大魔王給的 「已刪除最大和最小值」的數據來分析
在原廠測試部分：
單因子變異數分析
摘要| 組 | 個數 | 總和 | 平均 | 變異數 | | NGK CR7E | 3 | 25.92 | 8.64 | 0.0016 | | NGK CR8EIX | 3 | 26.36 | 8.786667 | 3.33E-05 | | DENSO | 3 | 25.91 | 8.636667 | 0.000233 | | BOSCH | 3 | 26.44 | 8.813333 | 0.000633 | | VTA | 3 | 26.1 | 8.7 | 0.0001 | | NGK PT9 | 3 | 26.3 | 8.766667 | 0.003733 |
ANOVA
| 變源 | SS | 自由度 | MS | F | P-值 | 臨界值 | | 組間 | 0.086961 | 5 | 0.017392 | 16.47684 | 5.21E-05 | 3.105875 | | 組內 | 0.012667 | 12 | 0.001056 | | | | | | | | | | | | 總和 | 0.099628 | 17 | | | | |
F=16.48 (>3.11) 表示這六組的火星塞在原廠新競戰測試下是有差異的
但到底是差異在哪裡呢??
所以接下來以雙尾幫佛諾利檢定 (檢定過程略，太占版面了，因為要比較15次差異)
因為統計結果，無顯著差異居多，僅列出有顯著差異的部分
| t 檢定：兩個母體平均數差的檢定，假設變異數相等 | | | | | | NGK CR7E | NGK CR8EIX | | 平均數 | 8.64 | 8.786667 | | 變異數 | 0.0016 | 3.33E-05 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 0.000817 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | -6.28571 | | | P(T<=t) 單尾 | 0.001636 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.003272 | | | 臨界值：雙尾 | 6.254071 | |
| t 檢定：兩個母體平均數差的檢定，假設變異數相等 | | | | | | NGK CR7E | BOSCH | | 平均數 | 8.64 | 8.813333 | | 變異數 | 0.0016 | 0.000633 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 0.001117 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | -6.35281 | | | P(T<=t) 單尾 | 0.001573 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.003146 | | | 臨界值：雙尾 | 6.254071 | |
| t 檢定：兩個母體平均數差的檢定，假設變異數相等 | | | | | | NGK CR7E | BOSCH | | 平均數 | 8.64 | 8.813333 | | 變異數 | 0.0016 | 0.000633 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 0.001117 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | -6.35281 | | | P(T<=t) 單尾 | 0.001573 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.003146 | | | 臨界值：雙尾 | 6.254071 | |
| t 檢定：兩個母體平均數差的檢定，假設變異數相等 | | | | | | NGK CR8EIX | VTA | | 平均數 | 8.786667 | 8.7 | | 變異數 | 3.33E-05 | 0.0001 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 6.67E-05 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | 13 | | | P(T<=t) 單尾 | 0.000101 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.000202 | | | 臨界值：雙尾 | 6.254071 | |
| t 檢定：兩個母體平均數差的檢定，假設變異數相等 | | | | | | DENSO | BOSCH | | 平均數 | 8.636667 | 8.813333 | | 變異數 | 0.000233 | 0.000633 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 0.000433 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | -10.3942 | | | P(T<=t) 單尾 | 0.000242 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.000484 | | | 臨界值：雙尾 | 6.254071 | |
| t 檢定：兩個母體平均數差的檢定，假設變異數相等 | | | | | | DENSO | VTA | | 平均數 | 8.636667 | 8.7 | | 變異數 | 0.000233 | 0.0001 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 0.000167 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | -6.00833 | | | P(T<=t) 單尾 | 0.001931 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.003863 | | | 臨界值：雙尾 | 6.254071 | |
這部分比較有趣 雙尾檢定無差異，但單尾檢定是有差異的
| | BOSCH | VTA | | 平均數 | 8.813333 | 8.7 | | 變異數 | 0.000633 | 0.0001 | | 觀察值個數 | 3 | 3 | | Pooled 變異數 | 0.000367 | | | 假設的均數差 | 0 | | | 自由度 | 4 | | | t 統計 | 7.248824 | | | P(T<=t) 單尾 | 0.000961 | | | 臨界值：單尾 | 5.166801 | | | P(T<=t) 雙尾 | 0.001923 | | | 臨界值：雙尾 | 6.254071 | |
結論：
微觀來說差異不大，有可能是因為我沒有拿到原始數據
整體的話 BOSCH = NGK銥合金 = NGK賽車火星塞
但 BOSCH > VTA 也 > 原廠NGK 及 DENSO 火星塞
其餘的比較後 視為無差異，原因是因為測量變異(隨機誤差)所影響
粗糙的學理統計比較，以EXCEL進行統計
因為我電腦沒有安裝SPSS或SAS統計分析軟體，暫時以EXCEL分析
參考就好
|
