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