配對(duì)和兩樣本t檢驗(yàn)

配對(duì)和兩樣本t檢驗(yàn)第一頁(yè)，共二十九頁(yè)，編輯于2023年，星期五一、單樣本資料的

t檢驗(yàn)

目的:根據(jù)樣本均數(shù)推斷其總體均數(shù)是否與已知總體均數(shù)0相等——

與0的比較。應(yīng)用條件:獨(dú)立性正態(tài)性第二頁(yè)，共二十九頁(yè)，編輯于2023年，星期五例1某鎮(zhèn)痛新藥在藥廠投入量產(chǎn)，已知該藥的生產(chǎn)技術(shù)標(biāo)準(zhǔn)為：每片藥中平均有效藥物含量為40.4mg。為評(píng)價(jià)生產(chǎn)質(zhì)量是否達(dá)標(biāo)，現(xiàn)抽測(cè)10片藥物，結(jié)果為：41.3，40.8，38.7，42.7，43.0，38.6，42.7，39.5，42.1，42.1，問該廠是否達(dá)到生產(chǎn)要求？

第三頁(yè)，共二十九頁(yè)，編輯于2023年，星期五dataaa;inputa@@;x=a-40.4;cards;41.340.838.742.743.038.642.739.542.142.1;數(shù)據(jù)步第四頁(yè)，共二十九頁(yè)，編輯于2023年，星期五

過程步1調(diào)用means過程proc

means

mean

std

stderr

tprt;varx;run;第五頁(yè)，共二十九頁(yè)，編輯于2023年，星期五過程步2調(diào)用univariate過程proc

univariate

normal;varx;run;第六頁(yè)，共二十九頁(yè)，編輯于2023年，星期五

過程步3調(diào)用ttest過程proc

ttest

H0=40.4;/*H0=m：缺省情況下，H0=0；在單均數(shù)t檢驗(yàn)時(shí)，H0等于總體均數(shù)。*/vara;run;第七頁(yè)，共二十九頁(yè)，編輯于2023年，星期五方法1

方法2方法3第八頁(yè)，共二十九頁(yè)，編輯于2023年，星期五dataaa;inputa@@;x=a-40.4;cards;41.340.838.742.743.038.642.739.542.142.1；proc

means

mean

std

stderr

tprt;varx;run;proc

univariate

normal;varx;run;proc

ttest

h0=40.4;vara;run;第九頁(yè)，共二十九頁(yè)，編輯于2023年，星期五結(jié)果分析本例資料服從正態(tài)分布（W=0.876128，P=0.1178），故選t檢驗(yàn)。本例t=1.41，雙側(cè)檢驗(yàn)p=0.1917，按a=0.05水準(zhǔn)，不拒絕H0

，差別無統(tǒng)計(jì)學(xué)意義，尚不認(rèn)為該廠生產(chǎn)質(zhì)量不達(dá)標(biāo)。第十頁(yè)，共二十九頁(yè)，編輯于2023年，星期五二、配對(duì)資料t檢驗(yàn)資料類型：兩個(gè)同質(zhì)對(duì)象接受不同處理；同一受試對(duì)象分別接受不同的處理，同一受試對(duì)象處理前后。條件：差值d服從正態(tài)分布第十一頁(yè)，共二十九頁(yè)，編輯于2023年，星期五例2第十二頁(yè)，共二十九頁(yè)，編輯于2023年，星期五數(shù)據(jù)步dataaa;

inputx1x2@@;

d=x1-x2;

cards;

0.94 0.921.02 1.011.14 1.11

1.23 1.221.31 1.321.41 1.42

1.53 1.511.61 1.611.72 1.72

1.81 1.821.93 1.932.02 2.04

;第十三頁(yè)，共二十九頁(yè)，編輯于2023年，星期五

procmeanstprt；vard；run；procunivariatenormal；vard；run；procttest;pairedx1*x2;run;方法1方法2方法3第十四頁(yè)，共二十九頁(yè)，編輯于2023年，星期五方法一

方法二方法三第十五頁(yè)，共二十九頁(yè)，編輯于2023年，星期五結(jié)果解釋本例t＝0.771，P=0.4569，在α=0.05水平上不能拒絕H0，差別沒有統(tǒng)計(jì)學(xué)意義

，尚不能認(rèn)為兩法測(cè)定結(jié)果不同。

第十六頁(yè)，共二十九頁(yè)，編輯于2023年，星期五三、兩個(gè)獨(dú)立樣本的t檢驗(yàn)

應(yīng)用條件：（1）觀察值之間是獨(dú)立的；（2）每組觀察值來自正態(tài)分布的總體；（3）兩個(gè)獨(dú)立組的方差相等。第十七頁(yè)，共二十九頁(yè)，編輯于2023年，星期五例3：第十八頁(yè)，共二十九頁(yè)，編輯于2023年，星期五dataaa;inputxgroup@@;cards;241361251141261341231201151191142182202152222242212252272232;proc

univariate

normal;classgroup;varx;

run;proc

ttest;classgroup;varx;run;第十九頁(yè)，共二十九頁(yè)，編輯于2023年，星期五第二十頁(yè)，共二十九頁(yè)，編輯于2023年，星期五F=2.92,P=0.1258,提示兩總體方差相等。第二十一頁(yè)，共二十九頁(yè)，編輯于2023年，星期五正態(tài)性檢驗(yàn)結(jié)果：

w1=0.93988，P=0.551；

w2=0.96219，P=0.8106；提示兩組資料均服從正態(tài)分布。方差齊性檢驗(yàn)結(jié)果：

F=2.92,P=0.1258,提示兩總體方差相等。t檢驗(yàn)結(jié)果：

t=1.02，P=0.3215，不拒絕H0

，還不能認(rèn)為兩組藥物的療效有差別。結(jié)果解釋：第二十二頁(yè)，共二十九頁(yè)，編輯于2023年，星期五作業(yè)：1、課本page105，第1題2、課本page106，第4題，再補(bǔ)充一個(gè)問題：試問B組放療前后的血清sil-2R水平（U/ml）有無差異？

第二十三頁(yè)，共二十九頁(yè)，編輯于2023年，星期五作業(yè)要求：1.計(jì)算均數(shù)、標(biāo)準(zhǔn)差、標(biāo)準(zhǔn)誤;2.寫出完整的假設(shè)檢驗(yàn)過程，其中結(jié)論要包括統(tǒng)計(jì)學(xué)和專業(yè)上的結(jié)論；3.寫出程序；4.寫出程序的結(jié)果，并對(duì)結(jié)果作結(jié)論。筆算機(jī)算第二十四頁(yè)，共二十九頁(yè)，編輯于2023年，星期五第二十五頁(yè)，共二十九頁(yè)，編輯于2023年，星期五單純放療放療+平消第二十六頁(yè)，共二十九頁(yè)，編輯于2023年，星期五ttest過程的格式PROCTTEST

選項(xiàng)串

CLASS

變量表

VAR

變量表

BY

變量表第二十七頁(yè)，共二十九頁(yè)，編輯于2023年，星期五語(yǔ)句說明

(一)PROCTTEST用于啟動(dòng)TTEST過程，是必須語(yǔ)句，先擇串可以是：1、DATA＝輸入文件名指明對(duì)哪個(gè)文件進(jìn)行t檢驗(yàn)。若省略，則SAS會(huì)自動(dòng)找出此程序之前最后形成的SAS數(shù)據(jù)文件，即當(dāng)前文件，對(duì)它進(jìn)行t檢驗(yàn)。2、ALPHA=p：確定檢驗(yàn)水準(zhǔn)α，如ALPHA=0.05，即檢驗(yàn)水準(zhǔn)為0.05。3、H0=m：缺省情況下，H0=0；在單均數(shù)t檢驗(yàn)時(shí)，H0等于總體均數(shù)。

第二十八頁(yè)，共二十九頁(yè)，編輯于2023年，星期五

(二)CLASS語(yǔ)句該語(yǔ)句的變量是分類變量，該變量只能取二值，分別代表所要比較的兩組資料，以分類變量分組，是必選語(yǔ)句