【正文】 X Source DF Sum of Squares Mean Square F Value Pr F Model 3 Error 12 Total 15 RSquare . Root MSE X Mean Source DF Anova SS Mean Square F Value Pr F A 1 B 1 A*B 1 不同方法、不同化合物間均有差別，且交互作用也有統(tǒng)計學(xué)意義。 新法高于舊法，乙物高于甲物，以新法 +乙物為最佳組合。 八、交叉試驗設(shè)計的方差分析 研究新藥催醒寧對麻醉劑的作用，將 6只大鼠配成 3對，隨機決定每對中 1只接受處理的順序（催醒寧或生理鹽水），等動物完全清醒 3天后將實驗組和對照組動物互換，同法重復(fù)試驗一次。結(jié)果如下： （ Y 藥物組 D 對照組） 動物編號 1 2 3 4 5 6 第 1次 Y 15 D 31 Y 17 D 30 D 25 Y 5 第 2次 D 27 Y 25 D 28 Y 14 Y 18 D 26 DATA AA13。 DO CX=1 TO 2。 DO DW=1 TO 6。 INPUT YW$ X 。 OUTPUT。 END。 END。 CARDS。 Y 15 D 31 Y 17 D 30 D 25 Y 5 D 27 Y 25 D 28 Y 14 Y 18 D 26 。 PROC ANOVA。 CLASS CX DW YW。 MODEL X=CX DW YW。 RUN。 Dependent Variable: X Source DF Sum of Squares Mean Square F Value Pr F Model 7 Error 4 Total 11 RSquare . Root MSE X Mean Source DF Anova SS Mean Square F Value Pr F CX 1 DW 5 YW 1 藥物間 P，拒絕 H0，接受 H1，可以認為實驗組與對照組不同，催醒寧有效。 方差分析 —— 非編程方式（分析員模塊） 菜單欄 Globals Analyze Statistical Analysis 菜單欄 File Open by SAS name (1)單因素方差分析 菜單欄 Statistics ANOVA Oneway ANOVA 變量 X選入 Dependent窗口，分組 G選入 Independent 窗口。 (2)兩因素、三因素及交叉設(shè)計的方差分析 菜單欄 Statistics ANOVA Factorial ANOVA 變量 X選入 Dependent窗口，分組選入 Independent 窗口。 需要作均數(shù)間兩兩比較時選 Means， 在 Comparison Mathod窗口中選兩兩比較的方法， 在 Main Effacts窗口選中分組變量， Add OK OK 析因設(shè)計的方差分析：交互作用分析選 Model, 在 Standard Models窗口中選 Twoway(兩因素 ) 或 Threeway(三因素 ) 第五節(jié) 定性資料假設(shè)檢驗 一、樣本率與總體率比較 沒有 SAS過程可以直接調(diào)用。 可以按照二項分布或 Poisson分布的概率函數(shù)在數(shù)據(jù)步中編程計算。 己知某藥治愈率為 60%。某醫(yī)生隨機抽取 10名病人試用此藥加 VC同服，結(jié)果有 8人治愈，問加服 VC是否有增效作用？ 該題己知 π=, n=10, x=8, 求 n≥8的概率。 可以用二項分布函數(shù) PROBBNML(π， n, x) DATA BB。 Q=PROBBNML(,10,7)。 P=1Q。 PROC PRINT。 RUN。 結(jié)果： OBS Q P 1 P，不拒絕 H0，尚不能認為 VC有增效作用 (1)二項分布直接概率法 (2)Poisson分布直接概率法 國家規(guī)定飲用水中大腸菌數(shù)平均不得超過 3個 /L?，F(xiàn)從某水源地抽取水樣 1L測得大腸菌 4個，問該水源是否符合飲用水條件？ 該資料服從 Poisson分布，己知 μ=3，求 X≥3的概率?？梢杂?Poisson分布函數(shù) POISSON (μ,x)求得。 DATA BB。 Q=POISSON(3,3)。 P=1Q。 PROC PRINT。 RUN。 結(jié)果： OBS Q P 1 P，不拒絕 H0，尚不能認為細菌數(shù)超標(biāo)。 (3) 大樣本正態(tài)近似法 — u檢驗 一般胃潰瘍患者中有 20%會發(fā)生胃出血癥狀。 某醫(yī)院觀察 65歲以上潰瘍病人 304例，有 96人發(fā)生胃出血。問老年胃潰瘍患者是否容易出血？ 己知 π=， n=304， x=96，可以用 u檢驗。 DATA BB。 INPUT P N X 。 S=SQRT(P*(1P)/N)。 U=ABS(X/NP)/S。 P1=1PROBNORM(U)。 P=2*P1。 U=ROUND(U,)。 P=ROUND(P,)。 IF P THEN P=。 FILE PRINT。 PUT 2 10 u= 25 P= 。CARDS。 304 96 。 RUN。 結(jié)果： u= P= 二、樣本率之間比較 — χ2檢驗 用兩種藥物治療十二腸潰瘍療效比較 藥物 愈合 未愈合 合計 愈合率 呋喃硝胺 54 8 62 甲氰咪胍 44 20 64 合計 98 28 126 DATA BB。 DO I=1 TO 2。 DO J=1 TO 2。 INPUT F。OUTPUT。 END。 END。 CARDS。 54 8 44 20 。 PROC FREQ。 TABLES I*J/CHISQ。 WEIGHT F。 RUN。 WEIGHT 語句指明變量 F 是頻數(shù) 加選項 EXACT可以對 2 2表作 Fisher確切概率 EXPECTED可以顯示每個格子的理論值 CMH可以計算隊列研究和病例對照研究的 RR值 NOROW、 NOCUL、 NOPERCENT 不顯示相應(yīng)百分數(shù) NOPRINT 不顯示表格，只顯示計算的統(tǒng)計量 χ2檢驗用 FREQ過程加 CHISQ選項 STATISTICS FOR TABLE OF I BY J Statistic 統(tǒng)計量 DF Value Prob ChiSquare χ2 1 Likelihood Ratio ChiSquare 似然比 χ2 1 Continuity Adj. ChiSquare 校正 χ2 1 MantelHaenszel ChiSquare MH χ2 1 Fisher‘ s Exact Test 確切概率 (Left) (Right) (2Tail) Phi Coefficient (Phi系數(shù) ) Contingency Coefficient (列聯(lián)系數(shù) ) Cramer‘ s V (克萊姆系數(shù) ) Sample Size = 126 由原始數(shù)據(jù)計算： 如 258人資料，三個定性變量 A文化程度（ 1=文盲 2=初中以下 3=高中 4=大學(xué)） B工作強度（ 1=輕 2=中 3=重） C月收入 （ 1=300 2=300 3=500 4=1000) 試分析文化程度與月收入之間的關(guān)系 DATA AA4。 INPUT A B C。 CARDS。 3 2 3 4 1 4 4 1 3…… 。 PROC FREQ。 TABLES A*C/CHISQ； RUN； 三、配對定性資料的 χ2檢驗 用甲乙兩種血清學(xué)方法檢查鼻咽癌患者血清 410份。問兩法的檢出率有無差別？ 甲法 乙法 + 合計 + 261 87 348 31 31 62 合計 292 118 410 DATA BB。 DO I=1 TO 2。 DO J=1 TO 2。 INPUT F。OUTPUT。 END。 END。 CARDS。 261 87 31 31 。 PROC FREQ。 TABLES I*J/AGREE。 WEIGHT F。 RUN。 AGREE選項作定性資料的一致性檢驗，計算 χ2和Kappa統(tǒng)計量， Kappa， 。 配對四格表資料計算 McNemar‘ s χ2統(tǒng)計量 配對 R C表 (方表 )計算 Bowker對稱性檢驗 QB統(tǒng)計量， QB近似服從 DF=R(R1)/2的 χ2分布。 當(dāng) P，拒絕 H0時，結(jié)論為不服從對稱分布。 結(jié)果： McNemar39。s Test Statistic = DF = 1 Prob = Simple Kappa Coefficient 95% Confidence Bounds Kappa = ASE = Sample Size = 4