1. ttest與ttest2的區(qū)別:
- 單樣本梦皮,配對(duì)樣本用 ttest
- 不配對(duì)樣本(獨(dú)立樣本)用ttest2
二者的公式也不一樣淘衙,但我不知道怎么打公式就不打了垃帅。
2. 導(dǎo)入數(shù)據(jù)
df = xlsread("C:\Users\Bai\Desktop\ex_corr.xlsx")
x=df(:1)
y=df(:2)
3.
[h,p,ci,stats] = ttest(x) ##數(shù)據(jù)均值是否等于0
[h,p,ci,stats] = ttest(x,m) ##數(shù)據(jù)均值是否等于m
[h,p,ci,stats] = ttest(x,y) ##兩個(gè)配對(duì)樣本
[h,p,ci,stats] = ttest2(x,y) ##兩個(gè)不配對(duì)樣本
h 0 or 1 接受/拒絕原假設(shè)
p pvalue
ci 置信區(qū)間
stats 統(tǒng)計(jì)量
tail: right 均值是否大于m ; left 小于m; both 等于m
4. vartype
help ttest2后铝阐,有一個(gè)參數(shù)值得注意:
vartype: equal or unequal
equal: 兩組的方差一致
unequal: 兩組的方差不一致
- 如何判斷兩組方查是否一致呢耗美,palisade給出很好的解釋腌乡,搬運(yùn)過來:
直接說結(jié)論:選擇unequal
原因:很多書和計(jì)算器用“pooling”這個(gè)術(shù)語盟劫,如果兩組的方差一致,你就可以“pool the dataset”与纽,假定他們是從同一個(gè)數(shù)據(jù)集中出來的侣签。
檢驗(yàn)兩組不同數(shù)據(jù)的方差是否相等塘装,需要用F檢驗(yàn),但如果F檢驗(yàn)出得出很大的P值影所,你只能說沒法成功拒絕原假設(shè)蹦肴,你并不能證明他。而且猴娩,F(xiàn)檢驗(yàn)需要兩組數(shù)據(jù)符合正態(tài)分布阴幌,不像T檢驗(yàn)大概是正態(tài)分布就可以了。你永遠(yuǎn)不知道什么時(shí)候是正態(tài)分布的卷中。所以保守的來講矛双,要選擇unequal。即使你選擇unequal蟆豫,pooling or not 只有微小的區(qū)別议忽。用unequal不需要你作出“我確定事實(shí)就是如此”的假設(shè),也不會(huì)對(duì)你的結(jié)果造成很大的影響十减。
原文:http://kb.palisade.com/index.php?pg=kb.page&id=1708
最后栈幸,學(xué)matblab隨手記感想:為什么實(shí)驗(yàn)室都用matlab我就也要從頭學(xué),這些東西明明用R已經(jīng)可以做的很好了帮辟。并且實(shí)驗(yàn)室不分享腳本速址,我就更找不到學(xué)他的理由了……