各位路过的大侠，帮帮忙，怎样在stata运行完命令加权最小二乘法vwls后，得到该方程的可决系数呢？
我自己有试过编程算出可决系数，但那结果好像是错的，下面是我编的算那可决系数的程序
vwls y x ,sd(x2)  /* 权重为1/x2 */
predict yhat
sum y
scalar py=r(mean)
display py
gen ts=(y-py)^2      /*总离差平方和*/
gen es=(yhat-py)^2   /*回归平方和 */
sum ts
scalar tss=r(sum)
display tss   /*总离差平方和*/
sum es
scalar ess=r(sum)
display ess   /*回归平方和 */
display ess/tss  /*可决系数*/