(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