資源描述:
《計量經(jīng)濟學(xué)2014秋 ans》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、計量2014秋作業(yè)1上交日期為2014/10/131.證明:因此,A不是無偏估計量,因此A是一致估計量。附加證明:證明:注意到,因此,因為,所以,對于(*)式,有2.在STATA里輸入命令setobs100,然后setseed2,然后genx=A*rnormal(),其中A是你學(xué)號的后兩位數(shù)字。a)生成y使得下面條件符合:把y回歸到x上,x的系數(shù)≈1,并且se(b1)≈0.2。題目取A=36為例。7代碼:clearsetobs100setseed2genx=36*rnormal()geny=x+72*rnormal()a)用reg命令驗證y符合條件。從回歸結(jié)果可以看出,x的系
2、數(shù)約為1,并且se(b1)約為0.2,符合條件。代碼:regyxb)運用reg命令得出的SST、SSE和SSR分別是多少,手動計算R方,其是否與輸出R方相同,解釋R方的含義。由回歸結(jié)果可知,SST=537078.194,SSE=118762.434,SSR=418315可以看出,手動算出的R方與回歸輸出的R方相同。R方的含義:表示變量y的22.11%變化可以被變量x解釋。代碼:di118762.434/537078.1947a)驗證OLS的三個代數(shù)性質(zhì)。性質(zhì)1:殘差和為零。性質(zhì)2:樣本與殘差的乘積和為零。注意到,和為零也及均值為零。故從下表可以看出,符合條件。性質(zhì)3:均值點在
3、擬合的直線上。如下圖所示,該性質(zhì)得到驗證。該題的另外做法是算出b0+b1*xbar,顯示其等于ybar代碼:regyxpredictfv,xbpredictu,residualegensum_u=sum(u)genxu=x*uegensum_xu=sum(xu)susum_usum_xusuyxtwoway(scatteryx)(lfityx),///xline(-1.240947)yline(6.724633)legend(off)///title("Therelationshipofyandx")7a)利用STATA的運算功能手動計算b1,并比較其與reg給出的系數(shù)的相似
4、性(提示:按照公式輸出的結(jié)果為:0.96332865,與回歸結(jié)果一致。代碼:egenxmean=mean(x)egenymean=mean(y)genX=x-xmeangenY=y-ymeangenXY=X*YgenXX=X*Xegensigma_XY=sum(XY)egensigma_XX=sum(XX)disigma_XY/sigma_XX2.homework.dta是一份來自于2010微觀實際調(diào)查數(shù)據(jù)的節(jié)選,請使用該數(shù)據(jù),探討收入影響因素:ytincome:個人去年全年的總收入,單位RMB元y_birth:出生年份年a)請用summarize命令列出相關(guān)變量的均值、標準
5、差、最大最小值。代碼:suy_birthsuytincomeb)微觀調(diào)查數(shù)據(jù)的處理中最需要注意的就是常常原始數(shù)據(jù)中被訪者未填答或者拒絕回答將用特殊的方式表達,這份數(shù)據(jù)中的ytincome中就是以9999997代表了不適用;9999998代表了不知道;9999999代表了拒絕回答。請你用柱狀圖命令hist大致看一下yincome的分布,用tab看一下這些不適用、不知道或者拒絕回答的樣本的個數(shù),并將他們的ytincome值設(shè)置為缺失。代碼:histytincometabytincomeifytincome<07replaceytincome=.ifytincome<0a)生成年齡
6、變量age,即樣本在被調(diào)查時候的年齡。注意缺失值的處理。代碼:replacey_birth=.ify_birth<0genage=2010-y_birthb)用stata的reg命令做ytincome和age的一元回歸分析,先寫出回歸模型再解釋方程估計系數(shù)b0和b1。b0表示:在剛出生,個人的與年齡無關(guān)的賺錢稟賦為年收入33487.4元。b1表示:在其他因素都不變的情況下,年齡每增長1歲,個人年收入會減少約302元。代碼:regytincomeagec)對上述回歸,輸出擬合值和殘差。請問有缺失值的樣本還會參與STATA的回歸么?有擬合值嗎?有殘差嗎?任何變量有缺失值的樣本不會
7、參與stata的回歸。解釋變量(age)有缺失的樣本沒有擬合值也沒有殘差;而被解釋變量(ytincome)有缺失的樣本有擬合值,但是沒有殘差。代碼:regytincomeagepredictfv,xbpredictu,residual7a)如果收入的單位改成了萬元,請預(yù)測b0和b1的系數(shù)將會有什么樣的改變?請用STATA驗證你的預(yù)測b0和b1都會縮小到原來的1/10000。預(yù)測得到了驗證。代碼:gennewincome=ytincome/10000regnewincomey_birthb)在實際處理數(shù)據(jù)