【正文】 。　　　　　　　　　　　　　　DO s=1 TO 4。 DO t=1 TO 5。DO s=1 TO 4。　　　　　　　　　　　　　　　 INPUT y @@。 OUTPUT?！　NPUT y1y5?！　　　　　　　　　　　?END。 END。 END。 END?！　UTPUT?！　　　　　　　　　　　　　?PROC GLM。END。 END。 END?！　　　　　　　　　　　　LASS a b s t。PROC GLM?！　　　　　　　　　　　　　　　ODEL y=a b a*b s(a b) t a*tCLASS a b?！　　　　　　　　　　　　　　　　　　　　*t a*b*t / ss3。MODEL y1y5=a b a*b / NOUNI。　　　　　　TEST H=a b a*b　E=s(a b)。REPEATED t 5 / SUMMARY PRINTE。　　　　　MEANS a*b a*t b*t a*b*t。RUN?！　　　　　　　　　　　　　　　　　UN?！　程序修改指導(dǎo)]　先將上表中16行5列試驗數(shù)據(jù)存在A盤上, 用　 INFILE語句調(diào)用試驗數(shù)據(jù)。在左邊的程序中,用了REPEATED語句, 仿上例可對因素t的5個水平之間實現(xiàn)5種兩兩比較。寫上PRINTE是為了對資料進行球性檢驗。當(dāng)此檢驗的結(jié)果為P,可用右邊的程序, 它可得出各交互作用項所對應(yīng)的均值等有關(guān)信息?！　∮嬎憬Y(jié)果的形式及結(jié)果的解釋,與上例類似,故從略。3. 具有2個重復(fù)測量的2因素設(shè)計　　[]　在“人混合淋巴細胞(MLC)培養(yǎng)試驗”中,用甲、乙２人的淋巴細胞混合后分別接受絲裂霉素和60Co處理,測出細胞反應(yīng)CPM值(CPM為每分鐘脈沖數(shù)),同時在相同條件下還測出甲細胞單獨的CPM值,用甲、乙細胞混合后的 CPM值除以單獨甲細胞的CPM值, 得一個相對的比值,用此比值作為觀測指標(biāo),將同一個樣品均分為２份,分別接受上述２種處理,然后,都在0、3周后4富同時間點上觀測相對比值,共用了6個樣品作試驗,設(shè)計和試驗數(shù)據(jù)如下表所示,試作分析。　　　　　　　　　　　人混合淋巴細胞經(jīng)絲裂霉素和60Co處理　　　　　　　后不同時間點上觀測的CPM的相對比值────────────────────────────────────────　樣　　　　　　　　　　　　CPM 值的　相　對　比　值　品　　───────────────────────────────────　號　　　　A1(T1　　　T2　　　T3　　 T4)　　A2(T1　　 T2　　 T3　　T4)────────────────────────────────────────　　1　　　　 　 　 　　　　 　　　 　　2　　　　 　 　 　 　　 　　　 　　3　　　　 　 　 　 　　 　　　 　　4　　　　 　 　 　 　　 　　　 　　5　　　　 　 　 　 　　　　　 　　6　　　　　 　 　 　　　　　 ────────────────────────────────────────　注: 因素A(處理),A1(絲裂霉素),A2(60Co)。因素T(時間),T1至T4分別為處理后0、3周觀測結(jié)果。[SAS程序] []　[]DATA a1?！　　　　　　　　ATA a2。　　　　　　　PROC GLM。INFILE 39。a:39。　　　 INFILE 39。a:39。　 CLASS a t s。DO s=1 TO 6。　　　　　　　DO s=1 TO 6?！　　　　ODEL y=a t a*s t*s a*t / ss3?！　NPUT y1y8?！　　　　?DO a=1 to 2?！　　　　EST H=a E=a*s。 TEST H=t E=t*s?！　UTPUT。　　　　　　　　DO t=1 to 4?！　　　　ONTRAST 39。(t1t4)/a139。END?！　　　　　　　　　　　?INPUT y @@?！　　　 1 1　0　0 a*t 1 1　0　0,PROC GLM?！　　　　　　　　　?OUTPUT?！　　　　　 1　0 1　0 a*t 1　0 1　0,MODEL y1y8= / NOUNI?！　?END。　　　　　　　　　t 1　0　0 1 a*t 1　0　0 1。REPEATED a 2, t 4?！　　　ND?！　　　　　　　　ONTRAST 39。(t1 vs t4)/a139。　　　　　/ PRINTE。RUN。　　END?！　　　　　　　　 1 0 0 1 a*t 1 0 0 1。 RUN。(程序一)　　　　　　　　　(程序二的第1部分)　　 (程序二的第2部分)　　[程序一修改指導(dǎo)]　A、T２因素都與重復(fù)測量有關(guān),因此,在過程步中不需寫CLASS語句,MODEL語句中等號右邊也就沒有分組變量了。當(dāng) REPEATED語句中有２個及２個以上因素時,把水平改變最慢的寫在前面,次慢的寫在第2位,依次類推,并在每個因素后面寫上水平數(shù), 用逗號與后面的因素隔開。因本設(shè)計中2個因素都與重復(fù)測量有關(guān),程序?qū)⒎謩e對每個因素及其交互作用項進行球性檢驗?！　程序二修改指導(dǎo)]　程序一不便對某因素的各水平進行兩兩比較,在程序二中就可方便的實現(xiàn)這一目的,但前提是資料應(yīng)通過球性檢驗。第1個CONTRAST語句是在A1條件下,對 T1～T4４個水平間作總的方差分析。第2個CONTRAST語句是在A1的條件下, 對T1與T4２水平間進行比較,進一步的修改法參見“析因設(shè)計及其資料的統(tǒng)計分析”一節(jié)。　　計算結(jié)果的形式及結(jié)果的解釋,與本節(jié)第2例類似,故從略。4. 具有2個重復(fù)測量的3因素設(shè)計　　[]　設(shè)第1個試驗因素A(藥物):A1(生理鹽水)、A2(H303藥,劑量 )、A3(復(fù)方丹參片,)。區(qū)組因素S(大鼠)。第2個試驗因素B(處理):B1(給藥)、B2(接扎冠脈前降支造成缺血)。第3個試驗因素C(時間,秒): C1(0～1)、C2(3)、C3(5～10)。觀測指標(biāo)為心率(分/次),表中給出的是以0時刻為基準(zhǔn)的相對心率。每1行中的6個數(shù)據(jù)重復(fù)測自同1只大鼠,試作分析(分析此資料的SAS程序見[])?！　　　　　　　　?A,B,C３個試驗因素對大鼠心率影響的測定結(jié)果────────────────────────────────────────　藥物　　鼠號　　　　　　　　　　以0時刻為基準(zhǔn)的相對心率(分/次)　　　　　　　 ──────────────────────────────　 A　　　 S　　B1(C1　　　C2　　　　 C3)　　　　B2(C1　　 C2　　　　C3)────────────────────────────────────────A1(生理　　1　　　　　 　　 　　　　　鹽水)　 2　　　　　 　　　　　　　　　　　 3　　　　　 　　　 　　　　　　　 4　　　　　 　　 　　　　A2(H303　　5　　　 　　 　　 　　　　　 　 6　　　　　 　　 　　　　　 /ml)　　7　　　 　　 　　 　　　　　　　　　 8　　　 　　 　　 　　　　A3(復(fù)方　　9　　　　　　　 　　　　　丹參　 10　　　 　　 　　 　　　　　 　11　　　 　　 　　 　　　　　 /ml)　 12　　　　　　　 　　　 ────────────────────────────────────────　注:。5. 具有2個重復(fù)測量的4因素設(shè)計　　[]　設(shè)第1個試驗因素A(藥物):A1(H303)、A2(復(fù)方丹參片)。第2個試驗因素B(藥物劑量,mg/ml):B1()、B2(2)。區(qū)組因素S(大鼠)。第3個試驗因素C(處理): C1(給藥)、C2(接扎冠脈前降支造成缺血)。第4個試驗因素D(時間,秒): D1(0～1)、D2(3)、D3(5～10)。觀測指標(biāo)為心率(分/次), 表中給出的是以0時刻為基準(zhǔn)的相對心率。每1行中的6個數(shù)據(jù)重復(fù)測自同1只大鼠,試作分析(分析此資料的SAS程序見[])。　　　　　　　　　 A,B,C,D４個試驗因素對大鼠心率影響的測定結(jié)果────────────────────────────────────────　藥物　　劑量　鼠號　　　　　　以0時刻為基準(zhǔn)的相對心率(分/次)　　　　　　　　　　 ─────────────────────────────　 A　　　 B　　　S　C1(D1　　　D2　　　　 D3)　　C2(D1　　 D2　　　　D3)────────────────────────────────────────A1(H303) B1()　1　　 　 　 　 　 　　　　　　　　　2　　　 　 　 　 　　　　　　　　　3　　 　 　 　 　 　　　　　　　　　4　　 　 　 　 　 　　　　 B2(2)　　5　　 　 　　 　 　　　　　　　　　6　　 　 　 　　　　　　　　　　　7　　 　 　 　 　 　　　　　　　　　8　　 　 　 　 　 A2(復(fù)方　B1()　9　　　　 　 　 　丹參)　　　　 10　　 　 　 　 　 　　　　　　　　 11　　 　 　 　 　 　　　　　　　　 12　　　　 　 　　　　　 B2(2)　 13　　 　 　 　 　 　　　　　　　　 14　　 　 　 　 　 　　　　　　　　 15　　 　 　 　 　 　　　　　　　　 16　　 　 　 　 　 ────────────────────────────────────────注:。　　[SAS程序]　[]　　　　　　　[]DATA a1?！　　　　　　　　　　　　　　　ATA a2。INFILE 39。a:39?！　　　　　　　　?INFILE 39。a:39。DO a=1 TO 3?！　　　　　　　　　　　　　O a=1 TO 2。DO s=1 TO 4。　　　　　　　　　　　　　　DO b=1 TO 2?！　NPUT y1y6。　　　　　　　　　　　　 DO s=1 TO 4?！　UTPUT?！　　　　　　　　　　　　　　　?INPUT y1y6。END。 END。　　　　　　　　　　　　　　　　　OUTPUT。PROC GLM?！　　　　　　　　　　　　　　　ND。 END。 END。CLASS a s?！　　　　　　　　　　　　　?PROC GLM。MODEL y1y6=a s / NOUNI SS3。　　　　　　CLASS a b s。REPEATED b 2, c 3　　　　　　　　　　　 MODEL y1y6=a b s / NOUNI SS3。　　　　　/ SUMMARY PRINTE?！　　　　　　EPEATED c 2, d 3RUN?！　　　　　　　　　　　　　　　　　　　　　? SUMMARY PRINTE?！UN?！　說明]　上述2個例子的程序修改指導(dǎo)及輸出結(jié)果的解釋從略,必要時可參看本節(jié)開始的幾個例子進行。