統(tǒng)計(jì)學(xué)5相關(guān)與回歸

醫(yī)學(xué)統(tǒng)計(jì)學(xué)實(shí)習(xí)課STATA軟件的應(yīng)用(5)--相關(guān)與回歸相關(guān)與回歸分析線性相關(guān)分析線性回歸分析等級(jí)相關(guān)分析列聯(lián)相關(guān)分析

21線性相關(guān)分析例1以下資料選自Galton的一項(xiàng)研究。探討成年時(shí)的身高是否與2歲時(shí)的身高存在相關(guān)。(ex_1.dta)單位：英寸

32歲身高3930323435363630成年身高71636367686870641線性相關(guān)分析考察變量間線性相關(guān)趨勢(shì)繪制兩個(gè)變量間的散點(diǎn)圖scatter變量1變量2scatteryx

42歲身高與成年身高的散點(diǎn)圖

51線性相關(guān)分析計(jì)算變量間的線性相關(guān)系數(shù)correlate變量列表correlateyx

61線性相關(guān)分析線性相關(guān)系數(shù)的假設(shè)檢驗(yàn)pwcorr變量列表,[選項(xiàng)]pwcorryx,sigstar(0.05)sig /*顯示假設(shè)檢驗(yàn)P值*/star(0.05) /*星號(hào)標(biāo)記P≤0.05者*/

71線性相關(guān)分析correlate變量名[，選擇項(xiàng)]注意：這里的變量可以是兩個(gè)或多個(gè)，corrlate命令給出的是變量間的兩兩的簡(jiǎn)單相關(guān)系數(shù)矩陣。選擇項(xiàng)：

means /*同時(shí)輸出均數(shù)、標(biāo)準(zhǔn)差等統(tǒng)計(jì)量

covariance /*不輸出相關(guān)系數(shù)矩陣，而輸出協(xié)方差

矩陣pwcorr變量名,[選擇項(xiàng)]sig打印假設(shè)檢驗(yàn)之P值star(#)如果相關(guān)系數(shù)的假設(shè)檢驗(yàn)之P值小于#，

則在相關(guān)系數(shù)旁打星號(hào)

2線性回歸分析例2請(qǐng)按例1提供的數(shù)據(jù)。能否利用2歲時(shí)的身高預(yù)測(cè)成年后的身高。已經(jīng)判斷兩個(gè)變量間存在線性相關(guān)建立模型建立2歲身高與成年身高間的回歸模型模型預(yù)測(cè)將2歲身高納入回歸模型預(yù)測(cè)成年身高

92線性回歸分析建立模型regress因變量自變量regressyx

102線性回歸分析

112線性回歸分析模型預(yù)測(cè)predict新變量,[選項(xiàng)]predictyhat,xb /*線性預(yù)測(cè)值*/predictresidual,residuals /*殘差*/predictse,stdp /*預(yù)測(cè)值標(biāo)準(zhǔn)誤*/predictsd,stdf /*觀察值標(biāo)準(zhǔn)差*/

122線性回歸分析估計(jì)預(yù)測(cè)值的95%可信區(qū)間genCI_L=yhat-invttail(6,0.025)*segenCI_U=yhat+invttail(6,0.025)*se估計(jì)觀察值的95%容許區(qū)間genR_L=yhat-invttail(6,0.025)*sdgenR_U=yhat+invttail(6,0.025)*sd

132線性回歸分析繪制回歸分析五線譜sortxscatteryx||lineyhatx||lineCI_Lx||lineCI_Ux||lineR_Lx||lineR_Ux

14回歸分析五線譜

153等級(jí)相關(guān)分析例3某研究所用野百合治療白血病，并作抗白指數(shù)及療效的分析，結(jié)果見(jiàn)表，問(wèn)抗白指數(shù)與臨床療效間有無(wú)關(guān)系？(ex_3.dta)

163等級(jí)相關(guān)分析spearmanxy

174列聯(lián)相關(guān)（定性資料的相關(guān)）tabvar1var2[fw=頻數(shù)變量][,V]4列聯(lián)相關(guān)例4某縣防疫站觀察三種藥物驅(qū)鉤蟲(chóng)的療效，在服藥后7天得糞檢鉤蟲(chóng)卵陰轉(zhuǎn)率（%）如下，問(wèn)三藥療效是否不同？(ex_4.dta)3種不同藥物與是否轉(zhuǎn)陰之間有無(wú)相關(guān)？tabdrugeffect[fw=f],V

19藥物陰轉(zhuǎn)例數(shù)未陰轉(zhuǎn)例數(shù)合計(jì)陰轉(zhuǎn)率(%)復(fù)方敵百蟲(chóng)片2893775.7純敵百蟲(chóng)片18203847.4滅蟲(chóng)靈10243429.4合計(jì)565310951.4直線回歸與t檢驗(yàn)、F檢驗(yàn)的關(guān)系(ex_5.dta)方案血紅蛋白增加量(g)A(n=20)1.80.52.33.72.42.01.52.71.10.91.41.22.30.70.51.41.73.03.22.5B(n=20)5.00.20.50.31.91.02.4-0.42.01.62.00.01.63.01.60.03.00.71.20.7用A和B組研究t檢驗(yàn)和直線回歸的關(guān)系ttestx,by(group)regxgroup

22練習(xí)(lx.dta)44名女醫(yī)學(xué)生的肺活量(L)和身高(cm)資料，試作如下分析:（1）以肺活量作為應(yīng)變量Y，身高作為自變量X，做回歸分析；（2）求相當(dāng)于各身高值的肺活量總體均數(shù)的95%可信區(qū)間及個(gè)體值的95%容許區(qū)間，并繪制回歸“五線譜”。

scatteryxregyxpredictyhatpredictseyhat,stdpgenl1=yhat-invttail(42,0.025)*seyhatgenl2=yhat+invttail(42,0.025)*seyhatpredictsey,stdfgenl3=yhat-invttail(42,0.025)*seygenl4=yhat+invtt