直線回和相關的Stata實現(xiàn)

1、文檔來源為：從網絡收集整理.word版本可編輯.歡迎下載支持第十章直線回歸和相關的 Stata實現(xiàn)本章使用的Stata命令為:直線回歸regress應變量 自變量殘差計算在執(zhí)行regress命令后，執(zhí)行卜列命令：predict 變量名，residualPearson 相關pwcorr 變量 1 變量 n, sigSpearman 相關spearman 變量 1 變量 2例10-1為了研究血清膽固醇含量與舒張壓之間是否存在依存關系，2006年在鄭州某大學隨機抽取10名成年男性，測得他們的血清膽固醇（mg/dL）含量和 舒張壓（mmHg）如表10-1,請作統(tǒng)計分析。表10-1 10名成人的血清膽固

2、醇（mg/dL）含量和舒張壓（mmHg）土匕右編12345678910膽固醇X307259341237254416267320374316舒張壓Y807590707510570858878H。：總體回歸方程不成立（0）H1：總體回歸方程成立（0）=0.05Stata數(shù)據：xy13078022597533419042377052547564161057267708320859374881031678Stata命令為:Model | 928.3431321 928.343132 Prob > F = 0.0000Residual | 114.0568688 14.2571086 R-squa

3、red = 0.8906+ Adj R-squared = 0.8769Total |1042.49 115.822222 Root MSE = 3.7759y | Coef. Std. Err. t P>|t|95% Conf. Interval+x |.177864.02204198.07 0.000.1270352 .2286927_cons |26.62225 6.9169953.85 0.00510.67163 42.57287t=8.07 ,P值0.001(Stata輸出值0.000),構建直線回歸方程將a和b代入式(10-2),可知,本例中，b的統(tǒng)計學意義為：血清膽固醇含量

4、每增加1mg/dL,總體中舒張壓平均增加0.178mmHg?？傮w均數(shù)yx的區(qū)間估計 給定X=X。時，Y的總體均數(shù)的點估計，例 101中，當自變量X取值為307 mg/dL時。個體丫值的容許區(qū)間估計 給定X = X0值時，估計總體中個體 Y值的波動范圍，以例 10-1中第一個樣本點的數(shù)據(307,80)為例。Stata命令：predict yhatpredict stdp, stdppredict stdf, stdfgenerate clm1 = yhat - invttail(7,0.025)*stdpgenerate clm2 = yhat + invttail(7,0.025)*stdp

5、generate clp1 = yhat - invttail(7,0.025)*stdfgenerate clp2 = yhat + invttail(7,0.025)*stdftwoway (lfitci y x, level(95) (scatter y x) (line clm1 clm2 clp1 clp2x, pstyle(p2 p2 p3 p3) sort)結果:xyyhatstdpstdfclm1clm2clp1clp23078081.231.193.9678.484.0571.8690.592597572.691.634.1168.8476.5362.9782.4134190

6、87.271.394.028490.5577.7696.782377068.781.994.2764.0873.4858.6978.872547571.81.74.1467.7775.836281.59416105100.612.644.6194.37106.8689.72111.512677074.111.514.0770.5477.6964.4983.733208583.541.223.9780.6686.4274.1692.923748893.141.864.2188.7497.5583.19103.13167882.831.23.9679.9885.6773.4692.2例10-2某發(fā)

7、熱門診醫(yī)生根據患者就診順序隨機抽取12名2040歲發(fā)熱患者，試探討體溫與脈搏之間的伴隨關系，數(shù)據見表10-4。表10-4 12名發(fā)熱患者的體溫（C）與脈搏（次/分）編號指標 123456789101112體溫X41.037.540.738.438.840.038.039.039.539.840.839.6脈才H Y13585134108120124106117125130138135Stata數(shù)據xy141135237.585340.7134438.4108538.8120640124738106839117939.51251039.81301140.81381239.6135H0：=0H1:

8、W0Stata命令： pwcorr y x,sig 結果： | y x + y | 1.0000I I x |0.91731.0000|0.0000PV0.0001,體溫與脈搏存在線性相關關系。例10 3某醫(yī)生研究新生兒缺氧缺血性腦?。╤ypoxic-ischemic encephalopathy HIE）臨床分度與其血漿 MDA含量（出生后第7天）的關系， 獲得了如下資料，試分析HIE臨床分度（1輕度，2中度，3重度）與血漿MDA 含量的相關方向及相關程度。HIE 臨床分度 1111122223333MDA（ umol/L） 7.4 7.0 6.7 6.3 8.17.2 8.4 8.98.8

9、 11.5 10.6 10.2 11.0本例目的在于分析HIE臨床分度與血漿MDA含量的關系，從數(shù)據特征上看， HIE臨床分度是等級資料，血漿MDA含量是連續(xù)性變量，因此應計算Spearman 等級相關系數(shù)rs。同樣地，人是樣本相關系數(shù)，僅為總體相關系數(shù) s的一個點估計，還應做總體相關系數(shù)s是否為0的假設檢驗，其無效假設與備擇假設分別為：H0： s 0Stata數(shù)據如下:ixy117.4217316.7416.3518.1627.2728.4828.9928.810311.511310.612310.213311Stata命令如下:spearman y x結果：Number of obs =13Spearman's