作者celestialgod (攸藍)
看板R_Language
標題Re: [問題] t檢定之t值相同p值不同問題請益
時間Mon Jun 1 18:12:29 2015
※ 引述《lapavoni (lapavoni)》之銘言：
: 文章分類提示:
: - 問題: 當你想要問問題時，請使用這個類別
: [問題類型]:
: 經驗諮詢(我覺得我的t檢定是不是出問題了...)
: [軟體熟悉度]:
: 使用者(已經有用R 做過不少作品)
: [問題敘述]:
: 目前在處理一份問卷分析，受試者目前依性別分為兩個Group(兩組人數不同)，
: 題目作答類型為五點計分，其中有兩題的t檢定跑出來的結果t值都是-2.0，但
: 是兩個p值差很多，一個有顯著，另一個則沒有，還有它df的算法我也不太懂...
: [程式範例]:
: ======================================================================
: 第八題的兩組作答資料為：
: > gb8
: [1] 4 3 2 2 4 3 4 3 3 2 3 4 2 1 4 2 3 4 2 3 1 2 3 3 4 3 2 2 2 3 2 2 3 3 2
2
: [70] 3 3 2 3 3 3 4 2 3 4 3 4 3 3 3 3 3 3 2 3 3 4 3 3 3 4 2 3 2 4 2 3 3 3 4
4
: [139] 4 2 2 3 2 1 2 3 2 3 2 4 2 2 3 3 2 3 3 3 4 2 4 1 3 3 3 1 3 3 3 3 2 3 2
3
: [208] 3 3 2 2 2 1 1 4 2 3 4 2 3 2 3 4 3 2 3 3 4 3 2 3 3 3 4 2 2 3 2 3 3 3 4
4
: [277] 3 2 3 3 3 3 3 3 3 3 4 2 2 3 3 4 3 3 3 3 4 3 3 2 4 3 2 2 2 3 3 4 3 4 3
2
: [346] 3 3 3 3 2 4 4 3 3 3 4 3 3 3 3 3 2 3 3 1 3 3 3 3 3 3 3 3 2 2 3 4 3 2 3
4
: [415] 3 4 4 4 3 3 3 3 3 2 4 3 4 3 3 3 4 4 4 3 3 3 2 3 3 3 2
: > sb8
: [1] 3 2 2 3 3 3 3 3 3 3 3 3 2 3 4 3 4 2 3 3 4 4 2 3 4 3 3 3 2 2 2 3 3 4 2 4
3
: [70] 3 3 2 3 3 3 2 3 3 3 4 3 3 3
: t檢定的結果為(兩組變異數有顯著差異)
: > t.test(gb8,sb8,var.equal = FALSE)
: Welch Two Sample t-test
: data: gb8 and sb8
: t = -2, df = 100, p-value = 0.1
: alternative hypothesis: true difference in means is not equal to 0
: 95 percent confidence interval:
: -0.2691 0.0259
: sample estimates:
: mean of x mean of y
: 2.90 3.02
我不知道你哪裏出問題，我把你資料下載下來之後，結果是這下面這樣
一個是-1.6，另外一個是-2.3....
嚴重懷疑變成四捨五入了...
gb8 = as.numeric(strsplit(
"4 3 2 2 4 3 4 3 3 2 3 4 2 1 4 2 3 4 2 3 1 2 3 3 4 3 2 2 2 3 2 2 3 3 2 2
1 3 4 4 3 3 4 3 4 3 3 2 3 4 3 3 2 4 4 4 3 3 3 4 2 4 3 2 2 2 3 2 3 3 3 2
3 3 3 4 2 3 4 3 4 3 3 3 3 3 3 2 3 3 4 3 3 3 4 2 3 2 4 2 3 3 3 4 4
4 3 2 3 3 3 3 3 2 3 3 3 3 3 3 3 3 2 3 3 3 1 4 2 3 3 3 3 3 3 2 2 3
4 2 2 3 2 1 2 3 2 3 2 4 2 2 3 3 2 3 3 3 4 2 4 1 3 3 3 1 3 3 3 3 2 3 2 3
3 4 3 3 4 3 3 3 3 4 3 2 3 3 3 3 2 4 3 3 3 3 3 3 2 3 3 3 3 2 2 3 4
3 3 2 2 2 1 1 4 2 3 4 2 3 2 3 4 3 2 3 3 4 3 2 3 3 3 4 2 2 3 2 3 3 3 4 4
3 3 3 3 1 2 3 3 3 3 4 4 3 3 3 3 3 3 3 4 2 3 3 3 4 3 2 3 3 3 3 3 3
3 2 3 3 3 3 3 3 3 3 4 2 2 3 3 4 3 3 3 3 4 3 3 2 4 3 2 2 2 3 3 4 3 4 3 2
1 1 3 3 3 1 3 3 3 3 2 3 3 3 2 3 4 4 4 4 3 3 4 4 3 2 1 3 3 4 3 3 3
3 3 3 3 2 4 4 3 3 3 4 3 3 3 3 3 2 3 3 1 3 3 3 3 3 3 3 3 2 2 3 4 3 2 3 4
3 3 3 2 3 3 3 3 3 3 3 3 2 3 4 3 4 2 4 4 3 2 3 3 3 3 2 3 2 2 4 2 3
3 4 4 4 3 3 3 3 3 2 4 3 4 3 3 3 4 4 4 3 3 3 2 3 3 3 2", "\\s+")[[1]])
sb8 = as.numeric(strsplit(
"3 2 2 3 3 3 3 3 3 3 3 3 2 3 4 3 4 2 3 3 4 4 2 3 4 3 3 3 2 2 2 3 3 4 2 4 3
3 3 3 3 2 3 4 3 3 4 4 4 3 3 3 2 3 4 3 3 4 3 4 3 4 3 2 3 3 3 3 3
3 3 2 3 3 3 2 3 3 3 4 3 3 3", "\\s+")[[1]])
t.test(gb8,sb8,var.equal = FALSE)
Welch Two Sample t-test
data: gb8 and sb8
t = -1.6309, df = 129.19, p-value = 0.1053
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-0.26912095 0.02591684
sample estimates:
mean of x mean of y
2.902494 3.024096
: =========================================================================
: 第9題的兩組資料為
: > gb9
: [1] 4 3 2 2 3 3 4 3 3 2 3 4 2 1 3 2 4 4 2 3 2 3 3 3 3 2 3 2 2 3 2 3 3 3 2
1
: [70] 3 3 2 3 2 3 4 2 3 3 3 4 3 3 3 3 3 3 2 3 3 4 3 3 3 4 2 3 2 3 3 3 3 2 4
4
: [139] 4 2 2 3 2 2 2 2 3 4 3 4 2 2 3 3 2 3 3 3 4 2 4 1 3 2 2 2 3 3 3 2 2 3 2
3
: [208] 3 3 2 2 2 3 2 3 2 3 3 2 3 2 3 3 3 3 3 3 4 3 2 3 3 3 4 3 3 2 2 3 3 3 4
4
: [277] 3 2 2 3 3 3 3 3 3 3 3 3 2 3 3 4 3 3 2 3 4 3 2 3 4 3 2 3 2 3 3 4 3 3 3
2
: [346] 3 3 3 3 2 3 4 3 3 4 2 3 3 3 3 3 2 3 3 1 2 3 2 3 3 3 3 3 2 2 3 4 3 3 3
3
: [415] 3 4 3 4 3 3 3 3 3 3 3 3 4 3 3 3 4 4 3 3 3 3 2 3 3 3 2
: > sb9
: [1] 3 2 3 2 3 3 3 3 4 3 2 3 3 3 4 3 4 2 3 3 4 4 2 3 3 3 4 3 2 3 3 3 3 4 3 4
3
: [70] 3 3 3 3 3 3 2 3 3 3 4 3 3 3
: t檢定結果為
: > t.test(gb9,sb9,var.equal = FALSE)
: Welch Two Sample t-test
: data: gb9 and sb9
: t = -2, df = 100, p-value = 0.02
: alternative hypothesis: true difference in means is not equal to 0
: 95 percent confidence interval:
: -0.3043 -0.0265
: sample estimates:
: mean of x mean of y
: 2.87 3.04
gb9 = as.numeric(strsplit(
"4 3 2 2 3 3 4 3 3 2 3 4 2 1 3 2 4 4 2 3 2 3 3 3 3 2 3 2 2 3 2 3 3 3 2 1
1 3 3 4 3 3 3 3 4 3 3 2 3 3 3 3 3 4 3 4 3 3 3 2 2 4 3 2 2 2 2 4 2
3 3 2 3 2 3 4 2 3 3 3 4 3 3 3 3 3 3 2 3 3 4 3 3 3 4 2 3 2 3 3 3 3 2 4 4
4 2 3 3 3 3 3 3 2 3 3 3 3 3 3 3 3 3 2 3 3 1 4 2 3 3 3 3 3 3 2 2 3
4 2 2 3 2 2 2 2 3 4 3 4 2 2 3 3 2 3 3 3 4 2 4 1 3 2 2 2 3 3 3 2 2 3 2 3
3 4 3 3 4 3 3 3 3 4 4 2 3 3 3 3 2 3 3 3 2 3 3 3 2 3 4 3 3 2 2 3 3
3 3 2 2 2 3 2 3 2 3 3 2 3 2 3 3 3 3 3 3 4 3 2 3 3 3 4 3 3 2 2 3 3 3 4 4
3 3 3 3 2 3 3 3 2 3 4 4 4 2 4 2 3 3 2 4 3 3 3 2 4 3 3 4 3 3 3 2 2
3 2 2 3 3 3 3 3 3 3 3 3 2 3 3 4 3 3 2 3 4 3 2 3 4 3 2 3 2 3 3 4 3 3 3 2
1 2 3 3 3 1 3 3 3 2 2 3 3 3 2 3 4 4 4 4 3 3 4 3 3 2 3 3 3 4 3 3 4
3 3 3 3 2 3 4 3 3 4 2 3 3 3 3 3 2 3 3 1 2 3 2 3 3 3 3 3 2 2 3 4 3 3 3 3
3 3 3 2 3 3 3 3 3 3 3 3 3 3 3 3 4 2 3 4 2 2 3 3 2 3 2 3 2 2 4 3 3
3 4 3 4 3 3 3 3 3 3 3 3 4 3 3 3 4 4 3 3 3 3 2 3 3 3 2", "\\s+")[[1]])
sb9 = as.numeric(strsplit(
"3 2 3 2 3 3 3 3 4 3 2 3 3 3 4 3 4 2 3 3 4 4 2 3 3 3 4 3 2 3 3 3 3 4 3 4
3 3 3 3 3 3 3 4 3 3 4 3 2 3 3 3 2 2 4 3 3 4 3 4 3 4 2 2 3 3 3 3 3
3 3 3 3 3 3 2 3 3 3 4 3 3 3 ", "\\s+")[[1]])
t.test(gb9,sb9,var.equal = FALSE)
Welch Two Sample t-test
data: gb9 and sb9
t = -2.356, df = 126.14, p-value = 0.02002
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-0.3043242 -0.0264684
sample estimates:
mean of x mean of y
2.870748 3.036145
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.205.27.107
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1433153553.A.0E5.html
※ 編輯: celestialgod (123.205.27.107), 06/01/2015 18:16:10
1F：推 lapavoni: 謝謝!!我再檢查一次後，真的是四捨五入問題XD，我錯了..06/01 23:18
2F：→ lapavoni: 還讓您費心跑一次，深感不安...XD06/01 23:27
有找到問題比較重要
※ 編輯: celestialgod (123.205.27.107), 06/02/2015 19:55:21