回歸分析實(shí)驗課-實(shí)驗

PAGEPAGE4實(shí)驗報告八實(shí)驗課程：回歸分析實(shí)驗課專業(yè)：統(tǒng)計學(xué)年級：姓名：學(xué)號：指導(dǎo)教師：完成時間：得分：教師評語：學(xué)生收獲與思考：26 92 021 175 030 31 022 104 00 277 012 210 019 120 04 290 016 238 028 164 115 272 111 295 138 68 131 85 121 224 120 166 113 305 130 124 114 246 1;run;procregdata=xt103;modely=x1x2;run;研制一種新型玻璃，對其做耐沖實(shí)驗。用一個小球從不同的高度h對玻璃做自由落體撞擊，玻璃破碎記為y=1，玻璃未破碎記y=0.數(shù)據(jù)見表22.是對表中數(shù)據(jù)建立玻璃耐沖性對高度h的logistic回歸，并解釋回歸方程的含義。模型信息：模型解出的是y=0的概率。由三個檢驗中，統(tǒng)計量的P值都小于0.05，可以認(rèn)為模型是顯著的。由Wald檢驗的顯著性概率及其P值，可以看出，h變量對方程的影響是顯著的。由極大似然估計，各個參數(shù)系數(shù)也通過檢驗。因此模型有效。二元logit模型為模型意義為，小球掉落高度為h，則玻璃未破碎的概率為p,而y=0表示玻璃未破碎。也就是說，該種新型的玻璃，用小球?qū)ζ渥矒?，?dāng)小球的掉落高度為h時，玻璃未破碎的概率就是，那么，玻璃會破碎的概率就為1-p(y=0)，這也可以看成是一種比例，就是大量實(shí)驗中，同個高度h，玻璃會被擊破的比例。SAS程序：datawjz;inputhy;/*引入虛擬變量，將公司類型的互助設(shè)為0，股份設(shè)為1*/cards;1.50 01.52 01.54 01.56 01.58 11.60 01.62 01.64 01.66 01.68 11.70 01.72 01.74 01.76 11.78 01.80 11.82 01.84 01.86 11.88 11.90 01.92 11.94 01.96 11.98 12.00 1;run;proclogisticdata=wjz;modely=h;run;proclogisticdata=wjz;classh;modely=h/link=glogitaggregatescale=none;run;某學(xué)校對本科畢業(yè)生的去向做了一個調(diào)查，分析影響畢業(yè)去向的相關(guān)因素，結(jié)果見表23.其中畢業(yè)去向“1”=工作，“2”=讀研，“3”=出國留學(xué)。性別“1”=男生，“0”=女生。用多類別的Logisitic回歸分析影響畢業(yè)去向的因素。專業(yè)課x1英語x2性別x3月生活費(fèi)x4畢業(yè)去向y兩個統(tǒng)計量的P值均大于0.05，說明模型擬合的較好。檢驗全局零假設(shè):BETA=0無效假設(shè)檢驗結(jié)果(似然比,評分）的結(jié)果P值均小于0.01，具有顯著統(tǒng)計學(xué)意義。三個變量中，有兩個是不顯著的變量，x3,x2,剔除x3:兩個統(tǒng)計量的P值均大于0.05，說明模型擬合的較好。檢驗全局零假設(shè):BETA=0無效假設(shè)檢驗結(jié)果(似然比,評分,wald）的結(jié)果P值均小于0.01，具有顯著統(tǒng)計學(xué)意義。三個變量都是顯著的。以x4=“1”，即參加工作，為參照。由模型可以看出：從參數(shù)估計表中，與參加工作的同學(xué)相比，讀研的（y=2）的同學(xué)相比，讀研的同學(xué)其專業(yè)課成績更好（x1的P值=0.003），而外語成績（x2的p值=0.356）和經(jīng)濟(jì)狀況（x4的P值=0.184）沒有顯著差異；出國留學(xué)的（y=3）學(xué)生其專業(yè)課成績和參加工作的沒有顯著差異，外語成績和經(jīng)濟(jì)狀況則更好。Sas程序：dataa;inputx1x2x3x4y;cards;95 65.0 1 600 263 62.0 0 850 182 53.0 0 700 260 88.0 0 850 372 65.0 1 750 185 85.0 0 1000 395 95.0 0 1200 292 92.0 1 950 263 63.0 0 850 178 75.0 1 900 190 78.0 0 500 182 83.0 1 750 280 65.0 1 850 383 75.0 0 600 260 90.0 0 650 375 90.0 1 800 263 83.0 1 700 185 75.0 0 750 273 86.0 0 950 286 66.0 1 1500 393 63.0 0 1300 273 72.0 0 850 186 60.0 1 950 276 63.0 0 1100 196 86.0 0 750 271 75.0 1 1000 163 72.0 1 850 260 88.0 0 650 167 95.0 1 500 186 93.0 0 550 163 76.0 0 650 186 86.0 0 750 276 85.0 1 650 182 92.0 1 950 373 60.0 0 800 182 85.0 1 750 275 75.0 0 750 172 63.0 1 650 181 88.0 0 850 392 96.0 1 950 2;run;procprint;run;proclogistic;classx3;modely(ref='3')=x1x2x3x4/link=glogitaggregatescale=none;run;proclogistic;classx3;modely(ref='3')=x1x2x4/link=glogitaggregatescale=none;run;proclogistic;classx3;modely(ref='1')=x1x2x4/link=glogitaggregatescale=none;run;六．收獲與思考七.思考題當(dāng)自變量是定性變量的時候，我們需要引進(jìn)虛擬變量進(jìn)行數(shù)量化，當(dāng)定性變量有n個水平的時候，我們該引進(jìn)多少的虛擬變量，否則會怎樣？不妨試試在sas中試試會出現(xiàn)什么問題。答：當(dāng)定性變量有n個水平時應(yīng)該引進(jìn)n-1個虛擬變量。否則最后一個虛擬變量無法用最小二乘估計計算出來。例：X1-X3為虛擬變量。Dataa;inputx1x2x3xy@@;cards;1001.26751001.35771001.40781001.58820101.71650101.76660101.80680101.85700011.22680011.35690011.46700011.4472;procregdata=a;modely=x1-x3x;run;X3沒有參數(shù)估計結(jié)果。因為x1x2x3出現(xiàn)完全共線性，x1x2均為0時即代表了x3為1.表21iyx1公司類型117151互助22692互助321175互助43031互助522104互助60277互助712210互助819120互助94290互助1016238互助1128164股份1215272股份1311295股份143868股份153185股份1621224股份1720166股份1813305股份1930124股份2014246股份表22序號h（m）y序號h（m）y11.500141.76121.520151.78031.540161.80141.560171.82051.581181.84061.600191.86171.620201.88181.640211.90091.660221.921101.681231.940111.700241.961121.720251.981131.740262.001表23序號專業(yè)課x1英語x2性別x3月生活費(fèi)x4畢業(yè)去向y19565.01600226362.00850138253.00700246088.00850357265.01750168585.001000379595.001200289292.01950296363.008501107875.019001119078.005001128283.017502138065.018503148375.006002156090.006503167590.018002176383.017001188575.007502197386.009502208666.0115003219363.0013002227372.008501238660.019502247663.0011001259686.00