單因素分析-定性資料2013.3.26

流行病學(xué)系廖巍單因素分析-定性資料第一部分：背景知識(shí)復(fù)習(xí)idagesmokeracelowbwt123030275023611028773261112350419020370052501035006211312200………………………………研究假設(shè)：母親的吸煙史是否與嬰兒發(fā)生低出生體重有關(guān)聯(lián)？χ2檢驗(yàn)如果有關(guān)聯(lián)：吸煙母親生出低出生體重嬰兒的風(fēng)險(xiǎn)是不吸煙母親的多少倍？計(jì)算RR或者OR值及相應(yīng)的可信區(qū)間(95%CI)流行病學(xué)資料的四格表形式病例對(duì)照合計(jì)暴露aba+b非暴露cdc+d合計(jì)a+cb+dn=a+b+c+d暴露:研究對(duì)象接觸過(guò)某種物質(zhì)、具備某種待研究的特征或行為。擴(kuò)展->R×C列聯(lián)表雙向無(wú)序單向有序雙向有序母親的吸煙史是否與嬰兒發(fā)生低出生體重有關(guān)聯(lián)？將資料整理成四格表形式：H0：π1=π2(兩樣本總體的率相等)H1：π1≠π2

(兩樣本總體的率不等)χ2檢驗(yàn)的基本思想病例對(duì)照合計(jì)低出生體重率(%)吸煙30(23.10)a44(50.90)b74(a+b)40.54π1=a/(a+b)不吸煙29(35.90)c86(79.10)d115(c+d)25.22π2=c/(c+d)合計(jì)59(a+c)130(b+d)189(n=a+b+c+d)31.22π=(a+c)/n四格表χ2統(tǒng)計(jì)量的校正R×C列聯(lián)表χ2統(tǒng)計(jì)量的校正1≤T<5格子數(shù)不超過(guò)1/5且沒(méi)有格子T<1超過(guò)1/5格子數(shù)1≤T<5，或有格子T<1:增加樣本量；專業(yè)允許情況下，將理論頻數(shù)過(guò)小的行或列與性質(zhì)相近的鄰行或鄰列合并；刪除理論頻數(shù)過(guò)小的行/列；無(wú)序R×C表用確切概率法。關(guān)聯(lián)強(qiáng)度的計(jì)算隊(duì)列研究：相對(duì)危險(xiǎn)度(RelativeRisk,RR)暴露組發(fā)病率與非暴露組發(fā)病率之比病例-對(duì)照研究：比值比(OddsRatio,OR)病例組的暴露比值=對(duì)照組的暴露比值=病例組的暴露比值與對(duì)照組的暴露比值之比RR（OR）值的含義RR(OR)=1暴露與疾病無(wú)關(guān)聯(lián)

RR(OR)>1暴露與疾病正向關(guān)聯(lián)暴露是疾病的危險(xiǎn)因素RR(OR)<1暴露與疾病負(fù)向關(guān)聯(lián)

暴露是疾病的保護(hù)因素RR（OR）值的95%可信區(qū)間（CI）含義若RR（OR）的95%CI包含1該RR（OR）值在α檢驗(yàn)水平上無(wú)顯著性若RR（OR）的95%CI不包含1該RR（OR）值在α檢驗(yàn)水平上有顯著性第二部分：主要SAS程序介紹過(guò)程步介紹數(shù)據(jù)步介紹SAS過(guò)程步介紹FREQ過(guò)程:產(chǎn)生一維至多維頻數(shù)表計(jì)算表中格子的理論頻數(shù)、構(gòu)成比和各種率對(duì)分類變量資料作相應(yīng)的假設(shè)檢驗(yàn)FREQ過(guò)程使用的語(yǔ)句PROCFREQ<選擇項(xiàng)>;TABLES

請(qǐng)求式/<選擇項(xiàng)>；WEIGHT

<變量名>;/*分析原始庫(kù)不用寫，分析頻數(shù)表必須寫*/RUN;TABLES語(yǔ)句TABLES

請(qǐng)求式/<選擇項(xiàng)>;請(qǐng)求式：要求FREQ過(guò)程分析處理的一維或多維表的清單一維表：tablesa；

由一個(gè)變量名表示；二維表：tablesa*b；由“*”連接兩個(gè)變量名表示。a-行變量；b-列變量

三維表：tablesa*b*c；由“*”連接三個(gè)變量名表示。a-分層變量；b-行變量；c-列變量……以此類推每個(gè)FREQ過(guò)程可以有多個(gè)TABLES語(yǔ)句每個(gè)TABLES語(yǔ)句有多個(gè)請(qǐng)求式procfreqdata=bw;tables

low;run;procfreqdata=bw;tables

smokelow;run;TABLES請(qǐng)求式--以lbw.dbf數(shù)據(jù)庫(kù)為例一維表格一維表格-多個(gè)請(qǐng)求式procfreqdata=bw;tables

race*smoke*low;run;procfreqdata=bw;tables

smoke*low;run;二維表格三維表格-分層TABLES語(yǔ)句常用選擇項(xiàng)：TABLES請(qǐng)求式…/<選擇項(xiàng)>;CHISQ:Pearsonχ2檢驗(yàn)似然比χ2檢驗(yàn)Mantel-

Haenszel

χ2檢驗(yàn)2×2表計(jì)算Fisher精確概率EXACT:對(duì)大于2×2的表,計(jì)算Fisher精確概率檢驗(yàn)AGREE:配對(duì)資料MEASURES:2×2表,提供OR、RR和相應(yīng)的可信區(qū)間。R×C表，衡量變量間的關(guān)系，計(jì)算Pearson和Spearman相關(guān)系數(shù)等CMH:分層分析2×2

表,提供OR、RR和相應(yīng)的可信區(qū)間。R×C列聯(lián)表，單向有序R×C列聯(lián)表，雙向有序ALL:=CHISQ+MEASURES+CMHEXPECTED：在獨(dú)立性假設(shè)下，單元格頻數(shù)的期望值（T）低出生體重正常體重合計(jì)吸煙304474不吸煙2986115合計(jì)59130189四格表資料dataa1;inputsmokelowcount@@;cards;1130104401290086;run;dataa2;dosmoke=1to2;dolow=1to2;inputcount@@;output;end;end;cards;30442986;run;smokelow001011000011000001…………數(shù)據(jù)庫(kù)資料行變量（R）、列變量（C）、頻數(shù)變量（F）或DO/END循環(huán)語(yǔ)句;里層優(yōu)先循環(huán);循環(huán)控制變量大于終值時(shí)終止output；SAS數(shù)據(jù)步介紹數(shù)據(jù)集形式根據(jù)四格表建立數(shù)據(jù)集需要幾個(gè)變量？數(shù)據(jù)集形式dataa2;dosmoke=1to2;dolow=1to2;inputcount@@;output;end;end;cards;30442986;run;procfreq;weightcount;tablessmoke*low;run;低出生體重正常體重合計(jì)吸煙304474不吸煙2986115合計(jì)59130189未使用output語(yǔ)句使用output語(yǔ)句Xdataa1;inputsmokelowcount@@;cards;1130104401290086;run;procfreq;weight

count;tablessmoke*low;run;procfreq;tablessmoke*low;run;正確輸出的表格未使用weight語(yǔ)句建立的數(shù)據(jù)集要在過(guò)程步中通過(guò)WEIGHT

<變量名>;給頻數(shù)變量加權(quán)!!!WEIGHT語(yǔ)句dataa1;inputsmokelowcount@@;cards;1130104401290086;run;procfreq;weightcount;tablessmoke*low;run;procfreqorder=data;weightcount;tablessmoke*low;run;低出生體重正常體重合計(jì)吸煙304474不吸煙2986115合計(jì)59130189第三部分：SAS分析實(shí)例常見單因素定性資料類型2×2列聯(lián)表配對(duì)資料精確概率法分層分析R×C列聯(lián)表某個(gè)病例對(duì)照研究擬探討“母親首次分娩的年齡與母親乳腺癌發(fā)病是否存在關(guān)聯(lián)”。研究中將母親首次分娩的年齡大于或等于30歲設(shè)定為暴露組，低于30歲設(shè)定為非暴露組。結(jié)果如下：探討首次分娩時(shí)年齡與乳腺癌是否存在關(guān)聯(lián)？2×2列聯(lián)表建立SAS數(shù)據(jù)集dataaa1;inputagecasecount;cards;11683

121498

212537

228747;run;dataaa2;doage=1to2;do

case=1to2;inputcount@@;output;end;end;cards;683149825378747;run;SAS過(guò)程步procfreqdata=aa1;tablesage*case/chisq;weightcount;run;procfreqdata=aa1;tablesage*case/expected;weightcount;run;procfreqdata=aa1;tablesage*case/measures;weightcount;run;procfreqdata=aa1;tablesage*case/all;weightcount;run;procfreqdata=aa1;tablesage*case;weightcount;run;某研究者隨機(jī)選取565名調(diào)查對(duì)象,其中血壓偏高者80名,血壓正常者485名,經(jīng)多年追蹤得到下表資料。試分析當(dāng)初血壓情況對(duì)調(diào)查對(duì)象后來(lái)否患冠心病可能造成什么樣的影響？2×2列聯(lián)表-練習(xí)1冠心病非冠心病血壓偏高1961血壓正常20465隊(duì)列研究配對(duì)資料公式及應(yīng)用條件-McNemar檢驗(yàn)b+c>40b+c<40乙法甲法+-+ab-cd配對(duì)資料例2：某醫(yī)院對(duì)168例痢疾標(biāo)本同時(shí)用SPAIF法和培養(yǎng)法（金標(biāo)準(zhǔn)）進(jìn)行檢驗(yàn)，用以分析兩種方法的檢驗(yàn)結(jié)果之間的差別有無(wú)顯著性意義，結(jié)果如下表。問(wèn)：兩種方法在檢驗(yàn)培養(yǎng)痢疾菌結(jié)果方面是否有差異？SAS過(guò)程步建立SAS數(shù)據(jù)集databb;doa=1to2;dob=1to2;inputcount@@;output;end;end;cards;64402044;run;procfreqdata=bb;tablesa*b/agree;weightcount;run;用兩種不同方法對(duì)53例肺癌患者進(jìn)行診斷，結(jié)果見下表。問(wèn)：兩種方法的檢測(cè)結(jié)果有無(wú)差別？配對(duì)資料-練習(xí)2甲法乙法合計(jì)+-+25227-111526合計(jì)361753SAS過(guò)程步建立SAS數(shù)據(jù)集datalx2;doa=1to2;dob=1to2;inputcount@@;output;end;end;cards;2521115;run;/*b+c<40校正公式*/procfreqdata=lx2;tablesa*b;exactMCNEM;weightcount;run;procfreqdata=lx2;tablesa*b/agree;weightcount;run;應(yīng)用條件：理論頻數(shù)T<1

或者資料總例數(shù)n<40

例3：研究三種藥物治療某病的結(jié)果，由于是少見疾病，所以例數(shù)很少（見下表）問(wèn)三種藥物治療效果有無(wú)差別？精確概率法(R×C表