ch5-方差分析I-完全隨機化析因設計

《ch5-方差分析I-完全隨機化析因設計》由會員上傳分享，免費在線閱讀，更多相關內容在教育資源-天天文庫。

1、上海財經大學統(tǒng)計與管理學院方差分析I-平衡設計問題和數據描述研究高血壓治療方案效果比較方法介紹方差分析、多重比較方差齊性檢驗參考資料汪仁官，陳榮昭（譯），蒙哥馬利（著），1998，實驗設計與分析，中國統(tǒng)計出版社。SAS分析tabulate過程anova過程glm過程gplot過程版權所有：吳純杰制表方差分析（平衡設計）廣義線性模型(不平衡數據的方差分析)作效應圖數據描述數據描述幾種高血壓治療方案效果研究。數據集：hyper變量：drug-X,Y,Z三種治療藥物biofeed-心理護理Present有，Absent無diet-特殊膳食Y有，N無

2、問題提出幾種治療方案是否有差異？治療、護理和膳食三因子是否存在交互效應？上海財經大學統(tǒng)管院數據描述上海財經大學統(tǒng)管院圖：hyper數據集數據錄入和整理datahyper;/*Chapter5*/infile'd:hypertension.dat';inputn1-n12;if_n_<4thenbiofeed='P';elsebiofeed='A';if_n_in(1,4)thendrug='X';if_n_in(2,5)thendrug='Y';if_n_in(3,6)thendrug='Z';arraynall{12}n1-n12;doi

3、=1to12;ifi>6thendiet='Y';elsediet='N';bp=nall{i};cell=drug

4、

5、biofeed

6、

7、diet;output;end;dropin1-n12;run;上海財經大學統(tǒng)管院注：原始數據為6行×12列條數據，即共有6條行觀測；其中1-3行對應著進行心理護理，而4-6行則無；而三種藥物分別為X在1和4行，Y在2和5行，Z在3和6行。每行的12條治療后的血壓數據使用數組nall存儲。通過循環(huán)將每行12條數據轉換成12行血壓數據，并且前6條對應輔助特殊膳食方案，后6條則無。單元cell存儲不同水平的組合

8、，共12種。注：_n_代表觀測值的序號；

9、

10、為字符串連接命令；output將讀入的數據和產生的新變量數據錄入到當前數據集；drop舍棄錄入數據所產生的輔助變量i和數組n1-n12。帶交互效應的方差分析模型上海財經大學統(tǒng)管院模型檢驗(F檢驗)獨立性正態(tài)性方差齊次性制表-tabulate過程procprintdata=hyper;run;proctabulatedata=hyper;classdrugdietbiofeed;varbp;tabledrug*diet*biofeed,bp*(meanstdn);run;注：對血壓bp按照drug,d

11、iet,Biofeed水平分類，并在表格中顯示各分類的均值、標準差和數據個數。思考：過程語句添加選項order=data，結果有何不同？上海財經大學統(tǒng)管院方差齊性檢驗procanovadata=hyper;classcell;modelbp=cell;meanscell/hovtest;run;注：對12個水平組合的方差通過/hovtest選項進行齊性檢驗；此處方差分析僅對不同水平處理的差異性進行檢驗。上海財經大學統(tǒng)管院Cell方差分析結果上海財經大學統(tǒng)管院模型顯著單元顯著差異cell方差齊性檢驗結果上海財經大學統(tǒng)管院p值>0.05,未通過方

12、差齊性檢驗，說明各水平組合觀察值的方差相等hyper方差分析procanovadata=hyper;classdietdrugbiofeed;modelbp=diet

13、drug

14、biofeed;meansdiet*drug*biofeed;odsoutputmeans=outmeans;run;procprintdata=outmeans;run;procsortdata=outmeans;bydrug;run;注：model語句等價于modelbp=dietdrugdiet*drugbiofeeddiet*biofeeddrug*biofe

15、eddiet*drug*biofeed上海財經大學統(tǒng)管院方差分析模型：做三因子主效應以及所有可能的交互效應模型（包括一階或二階三水平交互）。運算符

16、為邏輯或操作。means語句生成三個水平不同組合即cell的均值和方差等描述性結果。odsoutput為通過輸出傳遞系統(tǒng)將means語句所得結果輸出到outmeans數據集。方差分析結果上海財經大學統(tǒng)管院方差分析結果上海財經大學統(tǒng)管院結果顯示：三因子diet,drug和biofeed的主效應顯著（水平=0.05)，其中一階兩水平交互效應并不顯著，但三因子交互效應顯著，可能是由于一因子固定前提下，另

17、外兩因子交互效應存在很大差異造成。outmean輸出結果上海財經大學統(tǒng)管院按drug排序輸出繪制交互效應圖procsortdata=outmeans;bydrug;