sas04__sas基本統(tǒng)計(jì)分析功能-文庫(kù)吧
2025-04-20 12:33 本頁(yè)面
【正文】 句作用同前。 【 例 3】 ?某克山病區(qū)測(cè)得 11例克山病人與 13名健康人的血磷值（ mmol/L）如表，據(jù)此判斷該地急性克山病人與健康人的血磷值是否相同？ 輸出為： 分析： ? (1) 先作正態(tài)性檢驗(yàn)如下： ?零假設(shè)為 Ho：患者組的血磷值變量 x服從正態(tài)分布，其中患者組的 shapirowilk檢驗(yàn)的統(tǒng)計(jì)量為 w=，檢驗(yàn)的 p=α=，故應(yīng)接受零假設(shè)，即有 95%把握認(rèn)為 x正態(tài)。 ?零假設(shè)為 Ho：健康組的血磷值變量 x服從正態(tài)分布，其中健康組的 shapirowilk檢驗(yàn)的統(tǒng)計(jì)量為 w=，檢驗(yàn)的 p=α=，故應(yīng)接受零假設(shè)，即有 95%把握認(rèn)為 x正態(tài)。 分析： ? (2) 因此可采用兩獨(dú)立樣本均值的 T檢驗(yàn)。由 ttest過(guò)程輸出先作方差齊性檢驗(yàn)如下： ?Ｈ 0：患者組和健康組來(lái)自方差相等的總體，即，檢驗(yàn)的統(tǒng)計(jì)量 F=， P=α=，故應(yīng)接受零假設(shè)，即有95%把握認(rèn)為患者組和健康組方差滿足齊性。 ?再作 T檢驗(yàn)。 H0：患者組和健康組來(lái)自均值相等的總體，即。選擇方差齊性一行的結(jié)果知 t=,p=α=，故應(yīng)拒絕零假設(shè)，即有 95%把握認(rèn)為患者組與健康組血磷值的均值有顯著差異，且是患者組比健康組的均值要高。 注： ?使用 SAS軟件中的“分析家”，打開(kāi)數(shù)據(jù)集后，利用菜單“統(tǒng)計(jì)” → “假設(shè)檢驗(yàn)” → “均值的雙樣本 T檢驗(yàn)”，可以進(jìn)行單邊和雙邊檢驗(yàn)。 ?如果數(shù)據(jù)不服從正態(tài)分布可以采用非參數(shù)檢驗(yàn)，檢驗(yàn)兩個(gè)獨(dú)立樣本的中心位置是否相同的非參數(shù)方法有 wilxocon秩和方法，使用 npar1way過(guò)程加 wilcoxon選項(xiàng)，如教材 P123。 ?【 例 4】 檢驗(yàn)數(shù)據(jù)集 、女生的 gpa 分?jǐn)?shù)有無(wú)顯著差異？ 分析： ?在例 1中我們討論過(guò)變量 gpa是非正態(tài)分布，故要采用wilxocon秩和非參數(shù)檢驗(yàn)方法。 ?零假設(shè)為 H0：男生和女生來(lái)自中心位置相同的總體。由輸出結(jié)果的 wilcoxon秩和檢驗(yàn)中用正態(tài)近似得到的雙邊檢驗(yàn)的 z=， p=α=，故應(yīng)接受零假設(shè)，即有95%把握認(rèn)為男生和女生的 gpa無(wú)顯著性差異。 兩相關(guān)樣本均值檢驗(yàn)（ Univariate過(guò)程） ? ?相關(guān)樣本的有兩種情況 : ?一種是將研究對(duì)象按一定的條件先配對(duì)，每對(duì)中的兩個(gè)對(duì)象隨機(jī)分配到實(shí)驗(yàn)組和對(duì)照組，一個(gè)試驗(yàn)由若干對(duì)組成，稱為配對(duì)試驗(yàn)設(shè)計(jì)； ?另一種情況是同一批研究對(duì)象經(jīng)過(guò)某種處理前后的指標(biāo)值比較，或者是同一批樣品經(jīng)過(guò)兩種不同方法的測(cè)定結(jié)果的比較。此時(shí)這兩個(gè)變量不再獨(dú)立，而是相關(guān)的。 兩相關(guān)樣本均值檢驗(yàn)（ Univariate過(guò)程） ?檢驗(yàn)兩個(gè)相關(guān)變量的均值是否相等，等價(jià)于檢驗(yàn)這兩個(gè)變量間的差值變量的均值是否為零。當(dāng)差值變量服從正態(tài)分布時(shí)，可用配對(duì)樣本 T檢驗(yàn)的統(tǒng)計(jì)量為 ?其中 X為兩個(gè)樣本的差值變量， 、 S分別 為差值變量的均值和標(biāo)準(zhǔn)差。 0 /XtSn??X ?為了檢驗(yàn)兩個(gè)相關(guān)樣本的均值是否有顯著差異，先用一個(gè)數(shù)據(jù)步計(jì)算差值，然后對(duì)差值變量用 univariate過(guò)程可以實(shí)現(xiàn)檢驗(yàn)差值變量的均值是否顯著為零。 ?【 例 5】 用克矽平霧化吸入治療矽肺患者 7人，沒(méi)得治療前后的血清粘蛋白（ mg/L） 7對(duì)觀測(cè)值如表，據(jù)此能否認(rèn)為治療會(huì)引起血清蛋白的變化？（ α =） 輸出部分結(jié)果： 分析： ?先作正態(tài)性檢驗(yàn)如下： ?零假設(shè)為 Ho：差值變量 x服從正態(tài)分布，由輸出結(jié)果知shapirowilk檢驗(yàn)的統(tǒng)計(jì)量為 w=，檢驗(yàn)的p=α=，故應(yīng)接受零假設(shè)，即有 95%把握認(rèn)為 x正態(tài)。 ?故可采用兩相關(guān)樣本均值 T 檢驗(yàn) 。 ? H0：治療前后的差值變量 x的均值為 0。由輸出結(jié)果知 T檢驗(yàn)的統(tǒng)計(jì)量 t=，雙邊檢驗(yàn)的 p值為 α=，故拒絕原假設(shè)，即有 95%的把握認(rèn)為治療后血清蛋白有下降。 注： ?使用 SAS軟件中的“分析家”，打開(kāi)數(shù)據(jù)集后，利用菜單“統(tǒng)計(jì)” → “假設(shè)檢驗(yàn)” → “均值的雙樣本成對(duì) T檢驗(yàn)”，除了可以進(jìn)行單邊和雙邊檢驗(yàn)。 ?如果數(shù)據(jù)不服從正態(tài)分布可以采用非參數(shù)檢驗(yàn)，如符號(hào)檢驗(yàn)、符號(hào)秩檢驗(yàn)，在 SAS中可用 univariate過(guò)程實(shí)現(xiàn)。 ?【 例 6】 為了檢驗(yàn)一種新的復(fù)合肥料和原來(lái)使用的肥料相比是否顯著地提高了小麥的產(chǎn)量，在一個(gè)農(nóng)場(chǎng)中選擇了 10塊田地，每塊等分為兩部分，其中任指定一部分使用新的復(fù)合肥料，另一部分使用原肥料，小麥成熟后稱得各部分小麥的產(chǎn)量（單位： kg）如表，用符號(hào)檢驗(yàn)法檢驗(yàn)新復(fù)合肥是否會(huì)顯著提高小麥產(chǎn)量？（ α =） 分析： ?先作正態(tài)性檢驗(yàn)如下： ?零假設(shè)為 Ho：差值變量 diff服從正態(tài)分布，由輸出結(jié)果知shapirowilk檢驗(yàn)的統(tǒng)計(jì)量為 w=，檢驗(yàn)的p=α=，故應(yīng)拒絕零假設(shè)，即有 95%把握認(rèn)為差值變量 diff不服從正態(tài)分布。 ?故采用符號(hào)秩非參數(shù)檢驗(yàn)。 ? H0：差值變量 diff的均值為 0。由輸出結(jié)果知符號(hào)秩檢驗(yàn)的統(tǒng)計(jì)量 S=，雙邊檢驗(yàn)的 p值為 α=，故拒絕原假設(shè)，即有 95%的把握認(rèn)為新復(fù)合肥會(huì)顯著提高小麥產(chǎn)量。 作業(yè)： P157習(xí)題 1 , 3 , 5, 7 方差分析 ? ? ? ?思路 :檢驗(yàn)多組獨(dú)立樣本均值有無(wú)顯著性差異 ， 等價(jià)于檢驗(yàn)這個(gè)因素的各個(gè)取值水平會(huì)不會(huì)影響到指標(biāo)的取值 。 ? ?單因素方差分析模型 : ?其中 為分類變量（因子 A）的第 i個(gè)水平的效應(yīng)。 ?零假設(shè) ?平方和分解式： ?即總偏差平方和 =誤差的偏差平方和 +因子 A的偏差平方和 ?統(tǒng)計(jì)量 1, 1 , 2 , .. . , 1 , 2 , .. .0ij i ijriiy i r j n? ? ???? ? ? ? ???i?0. ..: 210 ???? rH ??? 0, .. .,: 211 不全為rH ???Aet SSS ??),1(~)()1( rnrFrnSrSFeA????? ?前提條件是獨(dú)立性、正態(tài)性和方差齊性。 ?當(dāng)數(shù)據(jù)滿足這些條件時(shí)，可采用 anova過(guò)程來(lái)進(jìn)行單因素方差分析，一般格式為： Proc anova data=數(shù)據(jù)集名； Class 因素變量； Model 指標(biāo)變量 =因素變量效應(yīng)表； Means 因素變量 / hovtest=levene ； Run； ?【 例 1】 （摘自魏宗舒 《 概率論與數(shù)理統(tǒng)計(jì) 》 例 ）為尋求適應(yīng)本地區(qū)的高產(chǎn)油菜品種，今選了五種不同品種進(jìn)行試驗(yàn)，每一品種在四塊試驗(yàn)田上試種，得到在每一塊田上的畝產(chǎn)量如表。 ?試問(wèn)不同品種的油菜的平均畝產(chǎn)是否相同？（ α =） A1 A2 A3 A4 A5
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1