統(tǒng)計SAS實例及答案

1、統(tǒng)計SAS實例及答案（專業(yè)結(jié)論待完善）目錄 TOC o 1-3 h z u HYPERLINK l _Toc422168267 第一部分 計量資料 PAGEREF _Toc422168267 h 2 HYPERLINK l _Toc422168268 （一）t檢驗 PAGEREF _Toc422168268 h 2 HYPERLINK l _Toc422168269 1、單樣本實例3-5 PAGEREF _Toc422168269 h 2 HYPERLINK l _Toc422168270 2、兩組配對實例3-6 PAGEREF _Toc422168270 h 2 HYPERLINK l _T

2、oc422168271 3、兩組獨立實例3-7 PAGEREF _Toc422168271 h 2 HYPERLINK l _Toc422168272 （二）方差分析 PAGEREF _Toc422168272 h 2 HYPERLINK l _Toc422168273 1、完全隨機實例4-2 PAGEREF _Toc422168273 h 2 HYPERLINK l _Toc422168274 2、隨機區(qū)組實例4-4 PAGEREF _Toc422168274 h 2 HYPERLINK l _Toc422168275 3、兩階段交叉實例4-6 PAGEREF _Toc422168275 h

3、 2 HYPERLINK l _Toc422168276 4、析因設(shè)計 PAGEREF _Toc422168276 h 2 HYPERLINK l _Toc422168277 22補充1個 PAGEREF _Toc422168277 h 2 HYPERLINK l _Toc422168278 22實例11-1 PAGEREF _Toc422168278 h 2 HYPERLINK l _Toc422168279 33實例11-2 PAGEREF _Toc422168279 h 2 HYPERLINK l _Toc422168280 522實例11-3 PAGEREF _Toc422168280

4、 h 2 HYPERLINK l _Toc422168281 5、正交設(shè)計 PAGEREF _Toc422168281 h 2 HYPERLINK l _Toc422168282 有空白列，無重復(fù)實例11-4 PAGEREF _Toc422168282 h 2 HYPERLINK l _Toc422168283 無空白列，有重復(fù)補充打印1 PAGEREF _Toc422168283 h 2 HYPERLINK l _Toc422168284 有空白列，有重復(fù)補充打印2 PAGEREF _Toc422168284 h 2 HYPERLINK l _Toc422168285 （三）非參 PAGER

5、EF _Toc422168285 h 2 HYPERLINK l _Toc422168286 1、單樣本類似3-5實例8-2 PAGEREF _Toc422168286 h 2 HYPERLINK l _Toc422168287 2、兩組配對類似3-6實例8-1 PAGEREF _Toc422168287 h 2 HYPERLINK l _Toc422168288 3、兩組獨立類似3-7實例8-3 PAGEREF _Toc422168288 h 2 HYPERLINK l _Toc422168289 4、多組非參實例8-5 PAGEREF _Toc422168289 h 2 HYPERLINK

6、 l _Toc422168290 第二部分 計數(shù)資料 PAGEREF _Toc422168290 h 2 HYPERLINK l _Toc422168291 （一）四格表（22） PAGEREF _Toc422168291 h 2 HYPERLINK l _Toc422168292 1、獨立四格表 PAGEREF _Toc422168292 h 2 HYPERLINK l _Toc422168293 基本卡方2實例7-1 PAGEREF _Toc422168293 h 2 HYPERLINK l _Toc422168294 連續(xù)校正c2或Fisher精確概率實例7-2 PAGEREF _Toc

7、422168294 h 2 HYPERLINK l _Toc422168295 Fisher精確概率實例7-4 PAGEREF _Toc422168295 h 2 HYPERLINK l _Toc422168296 2、配對四格表 PAGEREF _Toc422168296 h 2 HYPERLINK l _Toc422168297 Kappa值實例7-3 PAGEREF _Toc422168297 h 2 HYPERLINK l _Toc422168298 （二）行列表（rc） PAGEREF _Toc422168298 h 2 HYPERLINK l _Toc422168299 1、教材三

8、個 PAGEREF _Toc422168299 h 2 HYPERLINK l _Toc422168300 32率實例7-6 PAGEREF _Toc422168300 h 2 HYPERLINK l _Toc422168301 23構(gòu)成比實例7-7 PAGEREF _Toc422168301 h 2 HYPERLINK l _Toc422168302 43列聯(lián)系數(shù)c實例7-8 PAGEREF _Toc422168302 h 2 HYPERLINK l _Toc422168303 2、補充三個 PAGEREF _Toc422168303 h 2 HYPERLINK l _Toc42216830

9、4 23列有序-CMH2補充1個 PAGEREF _Toc422168304 h 2 HYPERLINK l _Toc422168305 33雙向有序-CMH2補充1個 PAGEREF _Toc422168305 h 2 HYPERLINK l _Toc422168306 43分層卡方-CMH2補充1個 PAGEREF _Toc422168306 h 2 HYPERLINK l _Toc422168307 （三）非參 PAGEREF _Toc422168307 h 2 HYPERLINK l _Toc422168308 25列有序-非參實例8-4 PAGEREF _Toc422168308 h

10、 2 HYPERLINK l _Toc422168309 44列有序-非參實例8-7 PAGEREF _Toc422168309 h 2 HYPERLINK l _Toc422168310 第三部分 回歸 PAGEREF _Toc422168310 h 2 HYPERLINK l _Toc422168311 （一）相關(guān) PAGEREF _Toc422168311 h 2 HYPERLINK l _Toc422168312 計量均正態(tài)-Pearson r（r）實例9-5 PAGEREF _Toc422168312 h 2 HYPERLINK l _Toc422168313 計數(shù)秩相關(guān)-Spear

11、man r（rs）實例9-8 PAGEREF _Toc422168313 h 2 HYPERLINK l _Toc422168314 （二）回歸 PAGEREF _Toc422168314 h 2 HYPERLINK l _Toc422168315 1、直線回歸 PAGEREF _Toc422168315 h 2 HYPERLINK l _Toc422168316 一元線性回歸實例9-1 PAGEREF _Toc422168316 h 2 HYPERLINK l _Toc422168317 多元線性回歸實例15-1 PAGEREF _Toc422168317 h 2 HYPERLINK l _

12、Toc422168318 2、曲線回歸 PAGEREF _Toc422168318 h 2 HYPERLINK l _Toc422168319 對數(shù)y=a+blog10(x)實例9-13 PAGEREF _Toc422168319 h 2 HYPERLINK l _Toc422168320 指數(shù)y=exp(a+bx)實例9-14 PAGEREF _Toc422168320 h 2 HYPERLINK l _Toc422168321 拋物線y=c+bx+axx補充 PAGEREF _Toc422168321 h 2 HYPERLINK l _Toc422168322 （三）回歸診斷-略 PAGE

13、REF _Toc422168322 h 2 HYPERLINK l _Toc422168323 （四）邏輯回歸 PAGEREF _Toc422168323 h 2 HYPERLINK l _Toc422168324 1、y二值 PAGEREF _Toc422168324 h 2 HYPERLINK l _Toc422168325 一元無啞變量補充 PAGEREF _Toc422168325 h 2 HYPERLINK l _Toc422168326 一元有啞變量補充 PAGEREF _Toc422168326 h 2 HYPERLINK l _Toc422168327 多元無交互作用補充 PA

14、GEREF _Toc422168327 h 2 HYPERLINK l _Toc422168328 多元有交互作用補充 PAGEREF _Toc422168328 h 2 HYPERLINK l _Toc422168329 多元有交互作用-篩選危險因子實例16-2 PAGEREF _Toc422168329 h 2 HYPERLINK l _Toc422168330 2、y多值-略 PAGEREF _Toc422168330 h 2 HYPERLINK l _Toc422168331 第四部分 其他 PAGEREF _Toc422168331 h 2 HYPERLINK l _Toc42216

15、8332 1、聚類-略 PAGEREF _Toc422168332 h 2 HYPERLINK l _Toc422168333 2、判別-略 PAGEREF _Toc422168333 h 2 HYPERLINK l _Toc422168334 3、主成分分析-略 PAGEREF _Toc422168334 h 2 HYPERLINK l _Toc422168335 4、因子分析補充 PAGEREF _Toc422168335 h 2實例第一部分 計量資料（一）t檢驗1、單樣本實例3-5原始數(shù)據(jù)：112 105 172 157 165 140 137 178 116 151 171 16312

16、9 130 125 135 128 100 126 128 90 113 128 12988 126 96 175 160 116 90 103 162 129 110 127SAS：data eg; input x; cards; 112127; proc univariate normal;/*正態(tài)性檢驗*/ var x; proc ttest h0=140;/*t檢驗，均值和140比較*/ var x; run;答案：H0：服從正態(tài)分布P=0.09880.05不拒絕H0，服從正態(tài)分布H0：等于140t=-2.14，P=0.0397P0.05，拒絕H0，差異有統(tǒng)計學(xué)意義x=130.8314

17、0因此，所有從事鉛工作男工人血紅蛋白含量均數(shù)低于正常人均數(shù)140g/L2、兩組配對實例3-6SAS：data eg; input x1 x2; d=x1-x2;cards;0.840.580.870.506;proc univariate normal; var d;/*對d值做正態(tài)性檢驗*/proc ttest; paired x1*x2;/*配對t檢驗*/proc means n mean; var x1 x2;/*算兩組的均值*/run;答案：H0：d值服從正態(tài)分布P=0.99930.10不拒絕H0，d值服從正態(tài)分布H0：兩種相同t=7.93，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義x1

18、=0.7952x2=0.5228因此，哥特里-羅紫法偏高3、兩組獨立實例3-7SAS：data eg; input c x;cards;1-0.72-2;proc univariate normal; var x; class c;/*按c進(jìn)行分組*/proc ttest; var x; class c;run;答案：兩組正態(tài)性檢驗：P=0.7527P=0.7754均0.10，均正態(tài)H0：方差齊Equality of VariancesP=0.31530.10不拒絕H0，方差齊H0：兩組相同t=-0.64，P=0.52480.05不拒絕H0，差異無統(tǒng)計學(xué)意義因此，尚不能認(rèn)為阿卡波糖與拜糖平對空

19、腹血糖降糖效果不同（二）方差分析1、完全隨機實例4-2SAS：data eg; input c x;cards;13.5343.71;proc univariate normal; var x; class c;proc anova; /*方差分析*/ class c; model x=c; /*x=影響因素*/ means c/hovtest welch;/*方差齊檢驗，不齊用welch*/ means c/snk; /*snk多重比較*/run;/*注意：t檢驗SAS默認(rèn)做方差齊檢驗，F(xiàn)檢驗（方差分析）要自己寫程序means c/hovtest*/答案：各組正態(tài)性檢驗：P=0.1947P=

20、0.4527P=0.1202P=0.2806均0.10，均正態(tài)H0：方差齊P=0.26290.10不拒絕H0，方差齊H0：四組相同F(xiàn)=24.88，P0.0001，拒絕H0差異有統(tǒng)計學(xué)意義，四組不全相同SNK多重比較：1、2、3、4即，安慰劑一組，2.4g與4.8g一組，7.2g一組2、隨機區(qū)組實例4-4SAS：data eg; input yw tz x;cards;110.82350.24;proc anova; class yw tz;/*按藥物、體重分組*/ model x=yw tz;/*x=藥物、體重兩個影響因素*/ means yw/snk;/*對藥物進(jìn)行多重比較*/run;/*注

21、意：本題樣本量太少，略去正態(tài)性檢驗和方差齊檢驗*/答案：H0：三種藥物相同F(xiàn)=11.94，P=0.00400.05拒絕H0，差異有統(tǒng)計學(xué)意義，三種藥物不全相同SNK多重比較：1、2、3因此，A藥一組，B藥和C藥一組H0：五個體重相同F(xiàn)=5.98，P=0.01580.05拒絕H0，差異有統(tǒng)計學(xué)意義，五個體重不全相同3、兩階段交叉實例4-6SAS：data eg; input jd ren yw$ x;/*yw$：藥物數(shù)據(jù)有字符，故加$*/cards;11A760210A803; proc anova; class jd ren yw;/*按階段、人、藥物分組*/ model x=yw jd re

22、n;/*x=階段、人、藥物三個影響因素*/run;答案：H0：兩種閃爍液相同F(xiàn)=4.02，P=0.07990.05拒絕H0，差異有統(tǒng)計學(xué)意義，兩種閃爍液不同H0：兩個階段相同F(xiàn)=9.93，P=0.01360.05拒絕H0，差異有統(tǒng)計學(xué)意義，兩個階段不同H0：10個人相同F(xiàn)=1240.19，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，10個人不全相同4、析因設(shè)計22補充1個例 缺鐵性貧血12例組1組2組3組4一般療法一般甲藥一般乙藥一般甲藥乙藥0.81.30.92.10.91.21.12.20.71.112SAS：data add; input a b x;cards; 000.8112.0;pr

23、oc anova; class a b; model x=a b a*b;/*a*b：考慮ab的交互作用*/ means a b a*b;/*算均數(shù)*/run;答案：一級分解：H0：4組相同F(xiàn)=98.75，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，四組不全相同二級分解：H0：甲的效應(yīng)為0F=168.75，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義用甲=1.2不用甲=0.8，甲的效應(yīng)為1.2-0.8=0.4，甲有效H0：乙的效應(yīng)為0F=90.75，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義用乙=1.0不用乙=0.8，乙的效應(yīng)為1.0-0.8=0.2，乙有效H0：甲乙的效應(yīng)為0F=36.75，P=0.

24、0003拒絕H0，差異有統(tǒng)計學(xué)意義甲乙同用=2.1，甲乙同用的效應(yīng)為2.1-0.8=1.3因此，甲乙的交互作用為1.3-0.4-0.2=0.7，甲乙存在22實例11-1SAS：data eg; input a b x;cards; 11102230;proc anova; class a b; model x=a b a*b; means b;run;答案：一級分解：H0：4組相同F(xiàn)=2.91，P=0.06660.05拒絕H0，差異有統(tǒng)計學(xué)意義，四組不全相同二級分解：H0：兩種縫合方法相同F(xiàn)=0.60，P=0.44990.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為兩種縫合方法不同H0：兩種縫

25、合后時間效果相同F(xiàn)=8.07，P=0.01180.05拒絕H0，差異有統(tǒng)計學(xué)意義，兩種縫合后時間效果不同一個月=26兩個月=48，因此縫合后2個月比1個月神經(jīng)軸突通過率高H0：縫合方法和縫合后時間無交互作用F=0.07，P=0.79950.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為縫合方法和縫合后時間存在交互作用33實例11-2SAS：data eg; input a b x;cards; 1110533160;proc anova; class a b; model x=a b a*b; means a b a*b/snk;run;答案：一級分解：H0：9組相同F(xiàn)=6.92，P=0.000

26、30.05拒絕H0，差異有統(tǒng)計學(xué)意義，9組不全相同二級分解：H0：A藥3種劑量效果相同F(xiàn)=8.47，P=0.00260.05拒絕H0，差異有統(tǒng)計學(xué)意義，A藥3種劑量效果不全相同SNK多重比較：1、2、3即，A藥1.0mg一組，2.5mg與5.0mg一組H0：B藥3種劑量效果相同F(xiàn)=9.05，P=0.00190.05拒絕H0，差異有統(tǒng)計學(xué)意義，B藥3種劑量效果不全相同SNK多重比較：1、2、3即，B藥5ug和15ug一組，30ug一組H0：A藥和B藥無交互作用F=5.07，P=0.00650.05拒絕H0，差異有統(tǒng)計學(xué)意義，A藥和B藥有交互作用最佳組合：A3和B3，即A藥5.0mg、B藥30ug

27、鎮(zhèn)痛效果最好522實例11-3見下頁：實例11-3：SAS：data eg; input b c a x;cards; 1110.252254.8;proc anova; class a b c; model x=a b c a*b a*c b*c a*b*c; means a/snk; run;答案：一級分解：H0：20組相同F(xiàn)=38.71，P0.001拒絕H0，差異有統(tǒng)計學(xué)意義，20組不全相同二級分解：H0：5種軍裝類型相同F(xiàn)=3.03，P=0.02210.05拒絕H0，差異有統(tǒng)計學(xué)意義，5種軍裝類型不全相同SNK多重比較：5、2、1、3、4即，第5種類型保溫，第4種類型散熱，其余三種介于

28、二者之間H0：2種活動環(huán)境相同F(xiàn)=23.08，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，2種活動環(huán)境不同H0：2種活動狀態(tài)相同F(xiàn)=659.10，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，2種活動狀態(tài)不同H0：軍裝類型與活動環(huán)境無交互作用F=1.13，P=0.34720.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為軍裝類型與活動環(huán)境有交互作用H0：軍裝類型與活動狀態(tài)無交互作用F=0.86，P=0.49080.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為軍裝類型與活動狀態(tài)有交互作用H0：活動環(huán)境與活動狀態(tài)無交互作用F=29.52，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，活動環(huán)境與活動狀態(tài)有交互作

29、用H0：軍裝類型與活動環(huán)境、活動狀態(tài)三者無交互作用F=0.94，P=0.44790.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為軍裝類型與活動環(huán)境、活動狀態(tài)三者有交互作用5、正交設(shè)計有空白列，無重復(fù)實例11-4見下頁：實例11-4：SAS：data eg; input a b ab c n1 n2 d x;cards; 111111186221211288;proc anova; class a b c d; model x=a b a*b c d; means a b c d a*b;run;答案：/*P.S.大檢驗P=0.08340.05，不用管，睜只眼閉只眼就直接跳過*/H0：A因素兩種

30、情況相同F(xiàn)=3.20，P=0.21550.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為A因素兩種情況不同H0：B因素兩種情況相同F(xiàn)=7.20，P=0.11530.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為B因素兩種情況不同H0：A因素B因素不存在F=20.00，P=0.04650.05拒絕H0，差異有統(tǒng)計學(xué)意義，A因素B因素存在均值比較，選擇：A2、B1H0：C因素兩種情況相同F(xiàn)=24.20，P=0.03890.05拒絕H0，差異有統(tǒng)計學(xué)意義，C因素兩種情況不同均值比較，選擇：C2H0：D因素兩種情況相同F(xiàn)=1.80，P=0.31180.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為D因素兩種

31、情況不同D1、D2可任選其一，選擇D2因此，最佳試驗條件：A2、B1、C2、D2此外，本次正交試驗未在該條件下進(jìn)行，建議補做無空白列，有重復(fù)補充打印1例 在胃蛋白酶生產(chǎn)過程中，產(chǎn)品質(zhì)量指標(biāo)是殘留蛋白酶，影響因素及其水平如下：因素水平水解溫度A水解時間B加鹽酸量C烘房溫度D1433.5h2.0%552464.0h2.6%603494.5h3.2%65選用L9（34）安排試驗，試驗結(jié)果如下（K3）：試驗號1A2B3C4D殘留蛋白酶111111.51.31.4212220.70.50.2313330.40.50.3421230.80.60.7522311.11.21.0623120.90.90.87

32、31320.70.70.8832130.30.30.4933213.23.03.1SAS：data add; input a b c d x;cards; 11111.533213.1;proc anova; class a b c d; model x=a b c d; means a b c d;run;答案：H0：27組相同F(xiàn)=159.99，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，27組不全相同H0：A因素3種相同F(xiàn)=73.22，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，A因素3種不全相同均值比較，均值越小殘留越少，選擇：A1H0：B因素3種相同F(xiàn)=113.19，P0.0001拒絕H0

33、，差異有統(tǒng)計學(xué)意義，B因素3種不全相同均值比較，選擇：B2H0：C因素3種相同F(xiàn)=85.70，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，C因素3種不全相同均值比較，選擇：C3H0：D因素3種相同F(xiàn)=367.86，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，D因素3種不全相同均值比較，選擇：D3因此，最佳試驗條件：A1、B2、C3、D3此外，本次正交試驗未在該條件下進(jìn)行，建議補做3遍有空白列，有重復(fù)補充打印2例 在正氟醚對大白鼠肝重影響的研究中，觀察指標(biāo)是肝重和體重之比?？紤]的因素及其水平如下：因素水平誘導(dǎo)藥物AB鼠的性別C1生理鹽水不用正氟醚雄性大白鼠2戊巴比妥用正氟醚雌性大白鼠需分析交互作用A

34、B，AC，BC，不考慮ABC，選用L8（27）安排試驗，每號試驗重復(fù)3次，即K3。試驗結(jié)果如下：試驗號1A2B3AB4C5AC6BC7肝重/體重（）111111115.00 5.52 5.38 211122225.26 5.68 5.83 312211226.30 7.02 5.90 412222115.42 5.60 5.70 521212126.13 6.46 5.21 621221215.87 5.50 6.20 722112216.02 5.70 5.48 822121124.64 4.60 5.44 SAS：data add; input a b ab c ac bc n x;car

35、ds;1111111522121125.44;proc anova; class a b c ; model x=a b c a*b a*c b*c; means a b c a*b a*c b*c;run;答案：H0：24組相同F(xiàn)=4.32，P0.0079拒絕H0，差異有統(tǒng)計學(xué)意義，24組不全相同H0：A因素兩種情況相同F(xiàn)=0.48，P0.49810.05拒絕H0，差異有統(tǒng)計學(xué)意義，A因素兩種情況不同H0：B因素兩種情況相同F(xiàn)=0.01，P0.91210.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為B因素兩種情況相同H0：C因素兩種情況相同F(xiàn)=4.97，P0.03950.05拒絕H0，差異有

36、統(tǒng)計學(xué)意義，C因素兩種情況不同H0：A因素B因素不存在F=11.84，P0.00310.05拒絕H0，差異有統(tǒng)計學(xué)意義，A因素B因素存在H0：A因素C因素不存在F=0.33，P0.57610.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為A因素C因素存在H0：B因素C因素不存在F=8.30，P0.01040.05拒絕H0，差異有統(tǒng)計學(xué)意義，B因素C因素存在/*P.S.本題不找最佳試驗條件，把關(guān)系分析出來就行了*/（三）非參1、單樣本類似3-5實例8-2SAS：data eg; input x1; d=x1-45.30; cards; 44.2187.37;proc univariate norm

37、al;var d;run;/*中位數(shù)比較：*/proc means p50;var x1;run;答案：H0：d正態(tài)P=0.0320.1，拒絕H0，非正態(tài)H0：二者相同S=31.5，P=0.00290.05，拒絕H0，差異有統(tǒng)計學(xué)意義，兩組不同中位數(shù)：53.2145.30因此，工廠高于正常人2、兩組配對類似3-6實例8-1SAS：data eg; input x1 x2; d=x1-x2; cards; 60 7695 100;proc univariate normal;var d;run;答案：H0：d正態(tài)P=0.07580.1，拒絕H0，非正態(tài)H0：兩組相同S=-21.5，P=0.056

38、60.05不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為兩組不同因此，尚不能認(rèn)為原法和新法測血清谷丙轉(zhuǎn)氨酶不同3、兩組獨立類似3-7實例8-3SAS：data eg; input x c;cards; 2.78 15.10 2;/*正態(tài)性檢驗：均正態(tài)*/proc univariate normal;var x;class c;run;/*方差齊檢驗：不齊取t的值*/proc ttest;var x;class c;run;/*作者不甘心，直接走非參，取單側(cè)值：*/proc npar1way wilcoxon;var x;class c;run;答案：H0：均正態(tài)P=0.7748P=0.5328均0.

39、1，不拒絕H0，均正態(tài)H0：方差齊F=16.84，P0.0001，拒絕H0，方差不齊H0：兩組相同若取t=0.05640.5，則不拒絕H0，尚不能認(rèn)為兩組不同故取非參單側(cè)值：P=0.04320.05，拒絕H0，差異有統(tǒng)計學(xué)意義，兩組不同因此，肺癌和矽肺RD值不同，由專業(yè)知識可知肺癌RD值高于矽肺4、多組非參實例8-5SAS：data eg; input c x;cards; 1 32.53 24.0;proc npar1way wilcoxon;var x;class c;run;答案：H0：三組全相同29.74，P=0.00770.05，拒絕H0，差異有統(tǒng)計學(xué)意義，三組不全相同秩均值比較：甲

40、=12.6乙=7.6丙=3.8因此三種藥物殺滅釘螺的效果有差別，甲藥死亡率最高，丙藥死亡率最低，乙藥居中第二部分 計數(shù)資料（一）四格表（22）1、獨立四格表基本卡方2實例7-1SAS：data eg; input r c f;/*f：頻數(shù)*/cards;1 1 992 2 21;proc freq;/*調(diào)用freq過程*/tables r*c/chisq nopercent nocol expected;/*chisq卡方檢驗；nopercent無總百分比；nocol無行百分比；expected理論值*/weight f;/* weight：權(quán)重，f為頻數(shù)變量*/run;答案：H0：兩組相同2

41、12.857，P=0.00030.05拒絕H0，差異有統(tǒng)計學(xué)意義，兩組不同因此，兩組降低顱內(nèi)壓的總體有效率有差別連續(xù)校正c2或Fisher精確概率實例7-2SAS：data eg; input r c f;cards; 1 1 462 2 8;proc freq;tables r*c/chisq nopercent nocol expected;weight f;run;答案：H0：兩組相等Fisher精確概率P0.05840.05/*只有一個T(理論值)5，故用Fisher精確概率，也可用連續(xù)校正卡方c2*/不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為兩組不等因此，尚不能認(rèn)為兩種藥物治療腦血管疾

42、病的有效率不等Fisher精確概率實例7-4SAS：data eg; input r c f;cards;1 1 42 2 6;proc freq;tables r*c/chisq nopercent nocol expected;weight f;run;答案：H0：兩組相同F(xiàn)isher精確概率P0.12100.05/*n=3340，故用Fisher精確概率*/不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為兩組不同因此，尚不能認(rèn)為兩組新生兒的HBV總體感染率有差別2、配對四格表Kappa值實例7-3SAS：data eg; input r c f;cards; 1 1 112 2 33;proc

43、freq;weight f;tables r*c/agree;/*agree：配對卡方檢驗*/run;答案：H0：兩組相同P=0.00750.05Kappa值0.45500.5拒絕H0，差異有統(tǒng)計學(xué)意義，兩組不同免疫熒光法陽性率：23/58乳膠凝集法陽性率：13/58因此，兩種方法的檢測結(jié)果有差別，免疫熒光法陽性率較高/*Kappa值：一致性檢驗介于0-1之間0：兩組不一致1：兩組一致區(qū)間：0-0.4較差0.4-0.75一般，最好過0.5，0.5認(rèn)為兩組一致0.75-1較好配對卡方重在檢測兩組的差異性，Kappa檢驗重在檢驗兩組的一致性*/（二）行列表（rc）1、教材三個32率實例7-6SAS

44、：data eg; input r c f;cards; 1 1 1993 2 26;proc freq;weight f;tables r*c/chisq nopercent nocol expected;run;答案：H0：三種相同221.0377，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，三組不全相同校正值：=/(n1)=0.05/4=0.0125H0：物理療法和藥物治療相同數(shù)據(jù)只留1、2組：26.7560，P=0.0093=0.0125拒絕H0，差異有統(tǒng)計學(xué)意義，物理療法和藥物治療不同H0：物理療法和外用膏藥相同數(shù)據(jù)只留1、3組：221.3228，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意

45、義，物理療法和外用膏藥不同H0：藥物治療和外用膏藥相同數(shù)據(jù)只留2、3組：24.5910，P=0.0321=0.0125不拒絕H0，差異無統(tǒng)計學(xué)意義，尚不能認(rèn)為藥物治療和外用膏藥不同/*四格表：SAS會把基本卡方2、連續(xù)校正卡方c2、Fisher精確概率都顯示出來，自個選；但，行列表Fisher精確概率要自己寫程序：/chisq exact*/23構(gòu)成比實例7-7SAS：data eg; input r c f;cards; 1 1 422 3 36;proc freq;weight f;tables r*c/chisq nopercent nocol expected;run;答案：H0：兩組

46、相同27.9127，P0.01910.05拒絕H0，差異有統(tǒng)計學(xué)意義，兩組不同因此，兩組2型糖尿病患者的ACE基因型總體分布有差別43列聯(lián)系數(shù)c實例7-8SAS：data eg; input r c f;cards; 1 1 4314 3 32;proc freq;weight f;tables r*c/chisq nopercent nocol expected;run;答案：H0：無關(guān)2213.1616，P0.0001拒絕H0，差異有統(tǒng)計學(xué)意義，有關(guān)列聯(lián)系數(shù)C=0.1883，關(guān)聯(lián)不大因此，兩種血型系統(tǒng)之間有關(guān)聯(lián)，但關(guān)聯(lián)不大/*列聯(lián)系數(shù)：Contingency CoefficientC介于0

47、-1之間0：獨立無關(guān)聯(lián)1：完全關(guān)聯(lián)*/2、補充三個23列有序-CMH2補充1個例 嬰兒營養(yǎng)狀況與副食品充足與否差中好充足73841不充足4206SAS：data add; input r c f;cards; 1 1 72 3 6;proc freq;weight f;tables r*c/CMH nopercent nocol;run;答案：H0：無關(guān)2CMH=5.9008，P=0.01510.05拒絕H0，差異有統(tǒng)計學(xué)意義，有關(guān)因此，嬰兒營養(yǎng)狀況與副食品充足與否有關(guān)/*CMH2三行結(jié)果：Nonzero Correlation雙向有序Row Mean Scores Differ 列有序 Ge

48、neral Association雙向無序*/33雙向有序-CMH2補充1個例 工齡與鉛中毒-短58144中32102長24128SAS：data add; input r c f;cards; 1 1 583 3 8;proc freq;weight f;tables r*c/CMH nopercent nocol;run;答案：H0：無關(guān)2CMH=7.3072，P=0.00690.05拒絕H0，差異有統(tǒng)計學(xué)意義，有關(guān)因此，工齡與鉛中毒程度有關(guān)43分層卡方-CMH2補充1個例 手術(shù)與術(shù)后并發(fā)癥無輕度中度醫(yī)院1手術(shù)12372手術(shù)220135醫(yī)院2手術(shù)11862手術(shù)213132SAS：data

49、add; input hospital trt effect f;cards; 111232232;proc freq order=data;tables hospital*trt*effect/CMH nopercent nocol;weight f;run;答案：H0：無關(guān)2CMH=4.2472，P=0.03930.05拒絕H0，差異有統(tǒng)計學(xué)意義，有關(guān)因此，剔除醫(yī)院影響之后，手術(shù)類型和并發(fā)癥有關(guān)（三）非參25列有序-非參實例8-4SAS：data eg; input c x f;cards; 1 1 12 5 0;proc npar1way wilcoxon;freq f;var x;cl

50、ass c;run;答案：H0：兩組相同Z=3.6971，P=0.00020.05，拒絕H0，差異有統(tǒng)計學(xué)意義，兩組不同秩均值比較：吸煙=49.15不吸煙=31.07因此，吸煙工人HbCO含量高于不吸煙工人44列有序-非參實例8-7SAS：data eg; input c x f;cards; 1 1 04 4 0;proc npar1way wilcoxon;freq f;var x;class c;run;答案：H0：四組全相同2=15.5058，P=0.00140.05，拒絕H0，差異有統(tǒng)計學(xué)意義，四組不全相同秩均值比較：支氣管擴張=43.50肺水腫=29.10肺癌=24.08病毒性呼吸

51、道感染=22.22因此，四種疾病患者痰液內(nèi)嗜酸性粒細(xì)胞有差別，支氣管擴張肺水腫肺癌病毒性呼吸道感染第三部分 回歸（一）相關(guān)計量均正態(tài)-Pearson r（r）實例9-5SAS：data eg;input x y;cards;43 217.2254 252.08;proc univariate normal;/*normal：正態(tài)性檢驗*/var x y;proc corr;/*corr：相關(guān)(correlation)*/var x y;run;答案：H0：x、y均正態(tài)P=0.8647P=0.5279均0.1，不拒絕H0，x、y均正態(tài)H0：無直線關(guān)系r=0.87543，P0.0001拒絕H0，有

52、直線關(guān)系，而且正相關(guān)/*r0正相關(guān)，r0負(fù)相關(guān)*/*x、y相關(guān)的條件：1、相關(guān)性檢驗P0.052、相關(guān)系數(shù)r絕對值不應(yīng)太小，至少0.45，最好過0.50-0.45：很差0.45-0.75：一般0.75-1：較強*/計數(shù)秩相關(guān)-Spearman r（rs）實例9-8SAS：data eg;input x y;cards;0.03 0.0527.96 9.33;/*計數(shù)資料不做正態(tài)性檢驗*/proc corr spearman;/*spearman相關(guān)系數(shù)(rs)*/var x y;run;答案：H0：無直線關(guān)系rs=0.90506，P0.0001拒絕H0，有直線關(guān)系，而且正相關(guān)（二）回歸1、直線

53、回歸一元線性回歸實例9-1SAS：data eg;input x y;cards;13 3.547 2.65;/*正態(tài)性檢驗：*/proc univariate normal;var x y;/*相關(guān)性檢驗：*/proc corr;var x y;/*回歸：*/proc reg;model y=x;/*畫圖：*/plot y*x;run;答案：H0：x、y均正態(tài)P=0.9332P=0.6122均0.1，不拒絕H0，x、y均正態(tài)H0：x、y無直線關(guān)系r=0.88177，P=0.00380.05，拒絕H0，x、y有直線關(guān)系，且正相關(guān)/*相關(guān)：x與y不分因果*/H0：回歸方程無意義（=0）F=20.

54、97，P=0.00380.05，拒絕H0，有直線回歸關(guān)系，回歸方程有意義/*回歸：y因變量，x自變量*/回歸方程：y=1.661670.13917x，r2=0.7775多元線性回歸實例15-1見下頁：實例15-1：SAS：data eg; input x1-x4 y ;cards;5.681.904.538.2011.203.841.206.459.6010.40;/*一步到位：*/proc reg;model y=x1-x4/stb selection=stepwise sle=0.1 sls=0.1;run;/*分別算x2、x3、x4與y的偏相關(guān)系數(shù)：*/proc corr nosimpl

55、e;var x2 y;partial x1 x3 x4;proc corr nosimple;var x3 y;partial x2 x1 x4;proc corr nosimple;var x4 y;partial x2 x3 x1;run;答案：H0：回歸方程無意義（1=2=3=4=0）F=11.41，P0.0001，拒絕H0，回歸方程有意義H0：2=0t=2.61，P=0.01560.05，拒絕H0，20H0：3=0t=-2.57，P=0.01710.05，拒絕H0，30H0：4=0t=2.88，P=0.00840.05，拒絕H0，40因此，標(biāo)準(zhǔn)化回歸方程：y=0.354x2-0.36x

56、30.413x4，調(diào)整的R2=0.5456偏相關(guān)系數(shù)：x2與y：0.34449x3與y：-0.42922x4與y：0.488162、曲線回歸對數(shù)y=a+blog10(x)實例9-13SAS：data eg;input x y;cards;0.005 34.1125.000 169.98;/*畫散點圖：*/proc gplot;plot y*x;run;/*線性排除：*/proc reg;model y=x;plot y*x;run;/*非線性回歸：*/proc nlin;parms a=0 b=0;model y=a+b*log10(x);run;答案：散點圖：/* 對數(shù)模型：y=a+b*lo

57、g10(x) */H0：x、y無線性關(guān)系F=7.54，P=0.07100.05，不拒絕H0，尚不能認(rèn)為x、y有線性關(guān)系H0：曲線方程無意義F=148.09，P=0.00120.05，拒絕H0，曲線方程有意義回歸方程：y=110.136.1154lgx指數(shù)y=exp(a+bx)實例9-14SAS：data eg;input x y;cards;2 5465 6;/*畫散點圖：*/proc gplot;plot y*x;run;/*線性排除：*/proc reg;model y=x;plot y*x;run;/*非線性回歸：*/proc nlin;parms a=0 b=0;model y=exp

58、(a+b*x);run;答案：散點圖：/* 指數(shù)模型：y=exp(a+b*x) */H0：x、y無線性關(guān)系F=100.62，P0.0001，拒絕H0，x、y有線性關(guān)系，r2=0.8856/*雖然r2很不錯了，但是散點圖看起來更像曲線，指數(shù)函數(shù)：*/H0：曲線方程無意義F=1585.01，P0.0001，拒絕H0，曲線方程有意義回歸方程：y=e4.0708-0.0396x拋物線y=c+bx+axx補充例 時間與紫外線強度時間9101112131415紫外線強度0.470.570.680.730.670.550.38SAS：data add;input x y;cards;9 0.4715 0.3

59、8;/*畫散點圖：*/proc gplot;plot y*x;run;/*線性排除：*/proc reg;model y=x;plot y*x;run;/*非線性回歸：*/proc nlin;parms a=0 b=0 c=0;model y=c+b*x+a*x*x;run;答案：散點圖： /* 拋物線模型：y=c+b*x+a*x*x */H0：x、y無線性關(guān)系F=0.20，P=0.67110.05，不拒絕H0，尚不能認(rèn)為x、y有線性關(guān)系H0：曲線方程無意義F=93.97，P=0.00040.05，拒絕H0，曲線方程有意義回歸方程：y=-0.03240.7657x-3.8176x2（三）回歸診

60、斷-略（四）邏輯回歸1、y二值一元無啞變量補充例 小孩出生體重與支氣管病體重組中值有病無病0-9507504919951-1350115018621351AS：data add;input y x wt;cards;1 750 490 1550 66;proc logistic descending;freq wt;model y=x/rsq cl;/*rsq：r2；cl：95%可信區(qū)間*/output out=a p=yp;proc print data=a;run;/* 模型：ln（p/(1-p)）=x */答案：H0：回歸方程無意義（=0）2=60.9739，P