(6) *--计算置信区间
use 计量经济学本科李子奈书数据\d3p72.dta,clear
reg y x1 x2
*--第一步计算X0*inv(X'*X)*X0',记为b1
gen cons = 1
mkmat cons x1 x2, mat(X) //把解释变量x1,x2,和常数项构建成矩阵X.
mat list X //显示X
matrix X0 = (1,20000,14000) //定义解释变量预测值所构成的矩阵X0,注意不是Xo
mat list X0
mat b1 = X0*inv(X'*X)*X0'
mat list b1~~~~~
*--2.直接在模型输出中找到MSE=385.92
*---最后3.计算置信区间
gen y0=143.3266+0.5556438*20000+0.2500854*14000
list y0
*--右端点:
genc2=14757.398+2.048*385.92*sqrt(1+0.30867531) //y的点预测值+/-为5%显著性水平t统计量*se*(1+矩阵计算值)的开方
*--其中14757.398为2006年因变量y的点预测值,
*--0.30867531为上面计算出来的X0*inv(X'*X)*X0',2.048为5%显著性水平t统计量,查表得到.(这个t值是怎么得到的???x1、x2不是应该有两个t值吗???)
*--左端点:
genc1=14757.398-2.048*385.92*sqrt(1+0.30867531)
list c1 c2