SAS練習(xí)題及程序答案復(fù)習(xí)過程

1、Good is good, but better carries it.精益求精，善益求善。SAS練習(xí)題及程序答案-隨機(jī)取組有無重復(fù)試驗的兩種本題是無重復(fù)DATAPGM15G;DOA=1TO4;/*A為窩別*/DOB=1TO3;/*B為雌激素劑量*/INPUTX;/*X為子宮重量*/OUTPUT;END;END;CARDS;106116145426811570111133426387;RUN;odshtml;/*將結(jié)果輸出成網(wǎng)頁格式，SAS9.0以后版本可用*/PROCGLMDATA=PGM15G;CLASSAB;MODELX=AB/SS3;MEANSAB;/*給出因素A、B各水平下的均值和標(biāo)

2、準(zhǔn)差*/MEANSB/SNK;/*對因素B（即劑量）各水平下的均值進(jìn)行兩兩比較*/RUN;ODSHTMLCLOSE;2*3析因設(shè)計兩因素完全隨機(jī)統(tǒng)計方法2*3析因設(shè)計tiff=f的開方DATAaaa;DOzs=125,200;DOrepeat=1TO2;/*每種試驗條件下有2次獨(dú)立重復(fù)試驗*/dojs=0.015,0.030,0.045;INPUTcl;OUTPUT;END;END;END;CARDS;2.702.452.602.782.492.722.832.852.862.862.802.87;run;PROCGLM;CLASSzsjs;MODELcl=zsjszs*js/SS3;MEAN

3、Szs*js;LSMEANSzs*js/TDIFFPDIFF;/*對zs和js各水平組合而成的試驗條件進(jìn)行均數(shù)進(jìn)行兩兩比較*/RUN;ODSHTMLCLOSE;練習(xí)一：2*2橫斷面研究列鏈表方法：卡方矯正卡方FISHERDATAPGM19A;DOA=1TO2;DOB=1TO2;INPUTF;OUTPUT;END;END;CARDS;226821;run;PROCFREQ;WEIGHTF;TABLESA*B/CHISQ;RUN;樣本大小=57練習(xí)二：對裂列連表結(jié)果變量換和不換三部曲1橫斷面研究P0.05RDATAPGM19B;DOA=1TO2;DOB=1TO2;INPUTF;OUTPUT;END

4、;END;CARDS;403414119252;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;樣本大小=57練習(xí)三：病例對照2*2病例組中有何沒有那個基因是正常的3.8倍，則有可能導(dǎo)致癡呆要做前瞻性研究用對裂DATAPGM20;DOA=1TO2;DOB=1TO2;INPUTF;OUTPUT;END;END;CARDS;24060360340;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQcmh;RUN;odshtmlclose;總樣本大小=1000練習(xí)四：配對設(shè)計隱含金

5、標(biāo)準(zhǔn)2*2MC卡方檢驗34和0在總體上（B+C40用矯正卡方）是否相等則可得甲培養(yǎng)基優(yōu)于乙培養(yǎng)基一般都用矯正因卡方為近似計算DATAPGM19F;INPUTbc;chi=(ABS(b-c)-1)*2/(b+c);p=1-PROBCHI(chi,1);求概率1減掉從左側(cè)積分到卡方的值chi=ROUND(chi,0.001);IFp0.0001THENp=ROUND(p,0.0001);FILEPRINT;PUT（打印在輸出床口）#210Chisq30Pvalue（#表示行）#410chi30p;CARDS;340;run;odshtmlclose;練習(xí)五：雙向有序R*C列連表用KPAdataaa

6、a;doa=1to3;dob=1to3;inputf;output;end;end;cards;582314278917;run;odshtml;*簡單kappa檢驗;procfreqdata=aaa;weightf;（頻數(shù)）tablesa*b;testkappa;run;*加權(quán)kappa檢驗;procfreq;weightf;tablesa*b;testwtkap;run;odshtmlclose;SAS系統(tǒng)FREQ過程頻數(shù)百分比行百分比列百分比a*b表ab合計12315839.4692.0686.5721.363.173.7732.044.7611.116342.86210.682.001

7、.494228.5784.0079.2574.7614.0025.935034.01385.4423.5311.9496.1226.4716.981711.5650.0062.963423.13合計6745.585336.052718.37147100.00a*b表的統(tǒng)計量對稱性檢驗統(tǒng)計量(S)2.8561自由度3PrS0.4144對稱性檢驗指總體上主對角線的上三角數(shù)相加是否與下三角三個數(shù)相加對稱性檢驗與KPA檢驗是否一致是否一個可以代替另一個檢驗Pe理論觀察一致率獨(dú)立假設(shè)性基礎(chǔ)上計算的相互獨(dú)立簡單Kappa系數(shù)Kappa0.6809漸近標(biāo)準(zhǔn)誤差0.050095%置信下限0.583095%置信

8、上限0.7788H0檢驗:Kappa=0總體的H0下的漸近標(biāo)準(zhǔn)誤差0.0597Z11.4112單側(cè)PrZ|Z|S0.4144簡單Kappa系數(shù)Kappa0.6809漸近標(biāo)準(zhǔn)誤差0.050095%置信下限0.583095%置信上限0.7788加權(quán)的Kappa系數(shù)加權(quán)的Kappa0.6614漸近標(biāo)準(zhǔn)誤差0.056095%置信下限0.551695%置信上限0.7711H0檢驗:加權(quán)的Kappa=0H0下的漸近標(biāo)準(zhǔn)誤差0.0646Z10.2406單側(cè)PrZ|Z|.0001對加權(quán)的KPA檢驗與簡單的（利用對角線上的數(shù)據(jù)分析）加權(quán)還要利用對角線以外的數(shù)據(jù)分析樣本大小=147練習(xí)六：雙向無序R*C列連表用卡

9、方理論頻數(shù)小于5沒有超過五分之一，一般用卡方實在不行用FISHER檢驗超過用KPA兩種血型都是按小中大排列相互不影響?yīng)毩⒌慕邮蹾0不一致行與列變量相互不影響DATAPGM20A;DOA=1TO4;DOB=1TO3;INPUTF;OUTPUT;END;END;CARDS;431490902388410800495587950137179325;run;odshtml;PROCFREQ;WEIGHTF;TABLESA*B/CHISQ;*exact;RUN;odshtmlclose;樣本大小=6094練習(xí)七：單向有序R*C秩和檢驗*方法1;（單因素非參數(shù)HO三個藥物療效相同H1不完全相等）DATAP

10、GM20C;DOA=1TO4;DOB=1TO3;INPUTF;OUTPUT;END;END;CARDS;15414991531504552224;run;odshtml;PROCNPAR1WAYWILCOXON;FREQF;CLASSB;VARA;RUN;*方法2;（FIQCHIM）procfreqdata=PGM20C;weightf;tablesb*a/cmhscores=rank;run;odshtmlclose;總樣本大小=270練習(xí)八：雙向有序?qū)傩圆煌琑*C4種目的4種方法SPEARMAN秩相關(guān)分析DATAPGM20E;DOA=1TO3;DOB=1TO3;INPUTF;OUTPUT;

11、END;END;CARDS;215131148671011284463132;run;odshtml;PROCCORRSPEARMAN;VARAB;FREQF;RUN;odshtmlclose;統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題一、定量資料上機(jī)實習(xí)題要求：先判斷定量資料所對應(yīng)的實驗設(shè)計類型；假定資料滿足參數(shù)檢驗的前提條件，請選用相應(yīng)設(shè)計的定量資料的方差分析，并用SAS軟件實現(xiàn)統(tǒng)計計算；摘錄主要計算結(jié)果并合理解釋，給出統(tǒng)計學(xué)結(jié)論和專業(yè)結(jié)論?！揪毩?xí)1】取4窩不同種系未成年的大白鼠，每窩3只，隨機(jī)分配到三個實驗組中，分別注射不同劑量雌激素，經(jīng)過一定時間后處死大白鼠測子宮重量，資料見表1。問劑量和窩

12、別的各自水平下子宮重量之間的差別有無統(tǒng)計學(xué)意義？若劑量間差別有統(tǒng)計學(xué)意義，請作兩兩比較。表1未成年大白鼠注射不同劑量雌激素后的子宮重量窩別子宮重量(g)劑量(g/100g):0.20.40.8合計1106116145367242681152253701111333144426387192合計2603584801098【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次課實驗課2.1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/G:中科院研究生統(tǒng)計學(xué)講課2011年有關(guān)資

13、料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生程序1.sas程序1【練習(xí)2】一位工程師研究由鉆頭壓力產(chǎn)生的沖力。考察了A（鉆孔速度）和B（進(jìn)料速度），兩因素分別取2與3水平，各水平組合下均做了兩次獨(dú)立重復(fù)實驗，資料見表2。假定資料滿足參數(shù)檢驗的前提條件，且兩因素對觀測結(jié)果的影響地位平等，已知沖力越小越好，試作分析，盡可能給出較為明確的統(tǒng)計和專業(yè)結(jié)論。表2在鉆孔速度和進(jìn)料速度取不同水平的條件下沖力的測定結(jié)果鉆孔沖力（單位）速度進(jìn)料速度：0.0150.0300.0451252.702.452.602.782.492.722002.832.852.862.862.802.

14、87【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次課實驗課2.1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/G:中科院研究生統(tǒng)計學(xué)講課2011年有關(guān)資料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生程序2.sas程序2.二、定性資料上機(jī)實習(xí)題要求：（1）若題目中未給出表格，請列出標(biāo)準(zhǔn)的列聯(lián)表，并對其命名；（2）若題目中已列出不規(guī)范的表格，先修改，然后對其命名；（3）根據(jù)分析目的或自己提出分析目的、資料的前提條件選用相應(yīng)的統(tǒng)計分析方法，

15、并用SAS軟件實現(xiàn)計算；（4）將主要計算結(jié)果摘錄出來，給出統(tǒng)計學(xué)和專業(yè)結(jié)論。【練習(xí)1】某衛(wèi)生防疫站對屠宰場及肉食零售點(diǎn)的豬肉,檢查其表層沙門氏菌帶菌情況，結(jié)果如下表。試比較屠宰場與肉食零售點(diǎn)豬肉表層沙門氏菌的帶菌率之間差別有無統(tǒng)計學(xué)意義?表1屠宰場及肉食零售點(diǎn)豬肉表層沙門氏菌抽檢結(jié)果采樣地點(diǎn)檢查數(shù)陽性數(shù)值帶菌率(%)屠宰場2827.14零售點(diǎn)29827.59【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次課實驗課2.1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/

16、G:中科院研究生統(tǒng)計學(xué)講課2011年有關(guān)資料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生練習(xí)1.sas練習(xí)1【練習(xí)2】有人對某部門22707名雇員中，普查了HBsAg,其中3454名陽性，19253名為陰性。從1975年起，追蹤了3年，發(fā)現(xiàn)在陽性組有40名患了肝癌，陰性組僅一名患肝癌。試選用合適的方法對資料進(jìn)行全面分析?！維AS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次課實驗課2.1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/G:中科院

17、研究生統(tǒng)計學(xué)講課2011年有關(guān)資料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生練習(xí)2.sas練習(xí)2【練習(xí)3】APOE-4等位基因與老年癡呆性的關(guān)聯(lián)研究：以600名晚發(fā)及散發(fā)老年癡呆患者和400名正常對照為研究對象，分析APOE-4等位基因與老年癡呆性的關(guān)系。表4APOE-4等位基因與老年癡呆性病例對照關(guān)聯(lián)研究AGT等位基因例數(shù)合計病例組對照組APOE-424060300非APOE-4360340700合計6004001000【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5次課實驗課2.

18、1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/G:中科院研究生統(tǒng)計學(xué)講課2011年有關(guān)資料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生練習(xí)3.sas練習(xí)3【練習(xí)4】請分析下表資料。已從專業(yè)上認(rèn)定培養(yǎng)的陽性結(jié)果就是“真陽性”，而不會出現(xiàn)假陽性。表3兩種培養(yǎng)基對同一批痰液標(biāo)本同時培養(yǎng)的結(jié)果甲培養(yǎng)基培養(yǎng)結(jié)果例數(shù)乙培養(yǎng)基結(jié)果：+-合計+363470-0135135合計36169205配對設(shè)計22列聯(lián)表資料總體率差異性檢驗統(tǒng)計量的計算公式若b+c40時若b+c40時【SAS程序】：HYPERLINKmhtml:file:/C:

19、Documents%20and%20Settingsdell桌面第5次課實驗課2.1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/G:中科院研究生統(tǒng)計學(xué)講課2011年有關(guān)資料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生練習(xí)4.sas練習(xí)4【練習(xí)5】請分析下表資料。表6兩法檢查室壁收縮運(yùn)動的符合情況對比法測冠心病人數(shù)定的結(jié)果核素法正常減弱異常正常5823減弱1427異常8917合計675327【SAS程序】：HYPERLINKmhtml:file:/C:Documents%20and%20Settingsdell桌面第5

20、次課實驗課2.1轉(zhuǎn)學(xué)生20121215中科院研究生院統(tǒng)計分析與SAS實現(xiàn)第1次上機(jī)實習(xí)題.mht!file:/G:中科院研究生統(tǒng)計學(xué)講課2011年有關(guān)資料2011年春季統(tǒng)計分析與SAS實現(xiàn)教學(xué)材料第5次課實驗實驗課2.1轉(zhuǎn)學(xué)生練習(xí)5.sas練習(xí)5簡單kappa檢驗和加權(quán)kappa檢驗這兩種方法都是用來檢驗兩種評價方法是否具有一致性的方法。其主要的區(qū)別是兩種方法計算的公式不一樣，更具體地說是對兩個變量的打分不一樣，簡單kappa檢驗主要是利用對角線上的信息，加權(quán)kappa檢驗除了利用對角線上的數(shù)據(jù)外，還將對角線外的數(shù)據(jù)進(jìn)行加權(quán)打分，將對角線外的信息也充分利用。所以在選擇方法時應(yīng)根據(jù)專業(yè)知識，如果兩個變量取值的界線比較明確，如“+”“+”“+”等，這時可以選用簡單的kappa檢驗，如果兩個變量的取值不十分明確，人為的因素較多時，就可選用加權(quán)kappa檢驗?！揪毩?xí)6】請分析下表資料。