1、Ming-chi Chen社會統計Stata教學Regression and CorrelationMing-chi Chen社會統計資料散佈圖scatter diagram 打開Chinalife.dta。 討論中國各省市的平均餘命（2000年）和其他地區變數（2004）的關係。 試畫出各省市每萬人病床數和平均餘命關係的scatter diagram。Ming-chi Chen社會統計Scatter DiagramIVDVMing-chi Chen社會統計PRC各省平均餘命與萬人病床數65707580平均餘命1020304050人均床位數（每萬人）Ming-chi Chen社會統計Stata

2、畫迴歸線迴歸線從迴歸線從0 0畫畫到極大值（到極大值（. .）Ming-chi Chen社會統計在散佈圖中畫出迴歸線6570758001020304050人均床位數（每萬人）Fitted values平均餘命教育投資有助於延長生命嗎？Ming-chi Chen社會統計Y-bar65707580平均餘命05001000150020002500人均教育支出前圖怎麼畫出來的？ twoway (scatter lifeexp pcedu) (lfit lifeexp pcedu, range(0 .), ytitle(平均餘命平均餘命, margin(medsmall) xtitle(人均教育支出人均

3、教育支出, margin(medsmall) legend(off) yline(71.24387) text(71.24387 -60 Y-bar, place(w)Ming-chi Chen社會統計散佈圖迴歸線Y軸標題中小大小的圖邊不要圖例加一條平均餘命的平均值（值是用summary計算出來的在（60, 71.24387）的位置左方（W西方）加上Y-bar的字樣可否用下拉式選單完成上圖？ 部分可以，但非全部Ming-chi Chen社會統計Ming-chi Chen社會統計Overlaid twoway graph：散佈圖Ming-chi Chen社會統計先選散佈圖IVDVOverlaid

4、 twoway graph：迴歸線Ming-chi Chen社會統計迴歸線從X0起到最大值Overlaid twoway graph：Y軸標題 Ming-chi Chen社會統計選擇中小型圖邊Overlaid twoway graph:X軸標題Ming-chi Chen社會統計Overlaid twoway graph:不要圖例Ming-chi Chen社會統計Overlaid twoway graph:產生圖表Ming-chi Chen社會統計65707580平均餘命05001000150020002500人均教育支出圖表缺少了代表平均餘命平均值的線，也沒有相應的註記，這可能是下拉式選單的缺

5、陷，也可能是我還沒找到正確的方法。所以要在command視窗裡加註。下拉式選單到指令視窗 把下拉式選單的指令在指令視窗裡叫回來。Ming-chi Chen社會統計在Review視窗裡有所有執行過的指令，點一下指令就會出現在指令視窗裡。在指令最後的地方加上yline和text的指令。自訂一條線並加上文字註記 twoway (scatter lifeexp pcedu) (lfit lifeexp pcedu, range(0 .), ytitle(平均餘命平均餘命, margin(medsmall) xtitle(人均教育支出人均教育支出, margin(medsmall) legend(off

6、) yline(71.24387) text(71.24387 -60 Y-bar, place(w)Ming-chi Chen社會統計-60是自己試驗出來的位置，因為如果設定為0的話，字會和軸線（不是在0上）重疊，看不清楚。迴歸分析Ming-chi Chen社會統計迴歸分析Ming-chi Chen社會統計依變數自變數迴歸分析：統計結果Ming-chi Chen社會統計解讀Stata迴歸分析結果Ming-chi Chen社會統計Stata以ANOVA表和迴歸係數表以及其他重要統計量來呈現迴歸統計結果：依變數Coef.Stf. Err.tp|t|係數的95%信賴區間自變數迴歸係數迴歸係數的標準

7、誤迴歸係數的檢定統計量tP值常數項截距截距的標準誤截距的檢定統計量 P值SourceSSdfMSNumber of obs=樣本數ModelSSRdf1=k-1MSRF(df1,df2)=模型F值ResidualSSEdf2=n-kMSER-squared=R2TotalSSTdf3Adj R-squared=調整R2=1 - (1 R2)(N - 1) /( N - k - 1). Root MSE=MSE的平方根ANOVA表迴歸係數表迴歸分析：詮釋統計結果Ming-chi Chen社會統計迴歸方程式：lifeexp=69.178+0.003587(pcedu)人均教育支出每增加人民幣一元，

8、該省市的平均餘命增加0.003587歲。R2=0.2965，即迴歸模型可以解釋29.65的總變異，FMSR/MSE12.22（t2），P值0.0015，我們可以拒絕整個迴歸模型沒有解釋力的虛無假設，也就是R2不等於0。迴歸係數與截距的顯著水準均在0.01以下迴歸係數的95信賴區間為（0.001489, 0.005685）。Ming-chi Chen社會統計迴歸模型DV的預測值在data editor中檢視新增加的預測值Ming-chi Chen社會統計預測值加入原資料中母體期望值E(Y|X)或y的95信賴區間 無法在下拉式選單中點選，必須在指令欄裡下指令。 對某特定Xp值而言，所預測的母體期望

9、值的估計變異數Ming-chi Chen社會統計)()(1)(222xxxxnsVaripey se2=MSE=7.4109， 值而不同。會隨著不同的pyXVar)(Stata中求對應特定Xp預測母體值的標準差Ming-chi Chen社會統計222)(1xnxxxnspe 在Stata裡用predict 新變數名稱, stdp這個指令來求對應數值。 predict stderr, stdpStata中求對應特定Xp預測母體值的標準差Ming-chi Chen社會統計母體期望值E(Y|X)或y的95信賴區間 先求出t值，要知道自由度在31-229下，=0.05的t值。 在Stata中，用inv

10、ttail(29, .05/2) generate yhatll=yhat-stderr*invttail(29, .05/2)這是信賴下界 generate yhatul=yhat+stderr*invttail(29, .05/2)這是信賴上界Ming-chi Chen社會統計母體期望值E(Y|X)或y的95信賴區間Ming-chi Chen社會統計圖形Ming-chi Chen社會統計圖形Ming-chi Chen社會統計母體預測值的信賴區間 從實際觀察到的樣本來推測母體，除了可以求出對應特定X值的母體期望值的信賴區間之外，我們還可以求母體預測值的信賴區間。 對某特定Xp值而言，母體預測

11、值的估計變異數Ming-chi Chen社會統計)()(11 )(222xxxxnsYVaripepStata中求對應特定Xp預測母體值的標準差Ming-chi Chen社會統計222)(11xnxxxnspe 在Stata裡用predict 新變數名稱, stdf這個指令來求對應數值。這裡和前面略有不同。 predict stderrf, stdfStata中求對應特定Xp預測母體值的標準差Ming-chi Chen社會統計 predict stderrf, stdf list母體期望值Y-hat的95信賴區間 先求出t值，要知道自由度在31-229下，=0.05的t值。 在Stata中，用

12、invttail(29, .05/2) generate yhatllf=yhat-stderrf*invttail(29, .05/2)這是信賴下界 generate yhatulf=yhat+stderrf*invttail(29, .05/2)這是信賴上界 list v1 lifeexp pcedu yhat stderrf yhatllf yhatulf Ming-chi Chen社會統計母體期望值Y-hat的95信賴區間Ming-chi Chen社會統計圖示Ming-chi Chen社會統計母體預測值95信賴區間圖示Ming-chi Chen社會統計Ming-chi Chen社會統計

13、Covariance and Correlation 如果只是要知道兩變數之間的相關程度，可以計算共變數或相關係數。Ming-chi Chen社會統計Covariance and CorrelationMing-chi Chen社會統計CorrelationMing-chi Chen社會統計Correlation . correlate lifeexp pcedu (obs=31) | lifeexp pcedu -+- lifeexp | 1.0000 pcedu | 0.5445 1.0000Ming-chi Chen社會統計CovarianceMing-chi Chen社會統計Covariance . correlate lifeexp pcedu, covariance (obs=31) | lifeexp pcedu -+- lifeexp | 10.1837 pcedu | 841.87 234698CovVarianceMing-chi Ch