测试发帖

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

发帖测试功能

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：测试

回归分析常用公式

补充参考

课上讲义安排

su R001-Dep08                               
                               
su R001-Dep08, detail
findit covar
corr R001 E0005
corr R001 E0005,c
findit coefficient
tabstat R001, stats(mean sd cv)       
help tabstat
tabstat R001, s(cv)
histogram R001
help histogram
hist R001,bin(6) start(0)
hist R001,bin(6) start(0) percent
hist R001,bin(15) start(0) percent
graph pie prefecture, sort
graph pie, over(prefecture)
gr pie, over(prefecture)
gr pie R001, over(prefecture)
gr pie R001, over(prefecture) plabel(_all percent)
gr pie R001, over(prefecture) plabel(_all percent) sort des ang(0)gr pie gr pie R001, over(prefecture) sort des ang(0) plabel(_all percent)                       
                               

// 清除先前变量的内存
clear all

// 开启日志记录
cap log use

// 当数据过长时关闭more按钮，防止程序中断
set more off

// 导入Excel数据
import excel "Z:\Personal\MBA\课程\数据模型与决策 DMD\作业\1127作业\上课用数据4.xlsx", sheet("Sheet1") firstrow clear

// 删除第一行（Delete the 1st row）
drop in 1

// 改变为数据格式
destring, replace

// 有条件地删除行
drop if code==.

// 给变量加标签

la var         code         "        当年行政区划编码        "
la var         county         "        当年行政区域名称        "
la var         code_2009         "        2009年行政区划编码        "
la var         county_2009         "        对应的2009年行政区域名称        "
la var         province         "        所在省代码        "
la var         prefecture         "        所在地市代码        "
la var         R001         "        一般预算收入合计        "
la var         R0002         "        增值税        "
la var         R0003         "        营业税        "
la var         R0005         "        个人所得税        "
la var         R0007         "        城市建设税        "
la var         R0021         "        农业五税        "
la var         R0024         "        企业收入        "
la var         E001         "        一般预算支出合计        "
la var         E0001         "        基本建设支出        "
la var         E1008         "        支援农村生产支出        "
la var         E1009         "        农林水气事业费        "
la var         E0005         "        教育事业费        "
la var         E1005         "        社会保障补助支出        "
la var         E1006         "        行政管理费        "
la var         E1007         "        公检法司支出        "
la var         R01         "        一般预算收入总计        "
la var         R111         "        返还性收入合计        "
la var         R1101         "        增值税、消费税税收返还补助        "
la var         R201         "        一般性转移支付收入小计        "
la var         R2001         "        原体制补助        "
la var         R301         "        专项补助        "
la var         R2002         "        转移支付补助        "
la var         R2011         "        增发国债补助        "
la var         R2004         "        增发工资补助        "
la var         R2008         "        各项结算补助        "
la var         R2024         "        调整收入任务增加或减少补助        "
la var         R401         "        上年结余收入        "
la var         R501         "        调入其他资金        "
la var         E01         "        一般预算支出总计        "
la var         E201         "        一般性转移支付上解合计        "
la var         E2001         "        原体制上解        "
la var         E301         "        专项上解        "
la var         E401         "        增设预算周转金        "
la var         E501         "        调出资金        "
la var         SP01         "        滚存结余        "
la var         SP02         "        其中：净结余        "
la var         FR001         "        收入合计        "
la var         FR0001         "        工业交通部门基金        "
la var         FR0002         "        文教部门基金         "
la var         FR0004         "        农业部门基金         "
la var         FR0005         "        土地有偿使用收入        "
la var         FR0007         "        政府住房基金收入        "
la var         FR0006         "        地方财政税费附加        "
la var         FE001         "        支出合计        "
la var         FE0001         "        工业交通部门基金         "
la var         FE0002         "        商贸部门基金        "
la var         FE0003         "        文教部门基金        "
la var         FE0005         "        农业部门基金        "
la var         FE0006         "        土地有偿使用支出        "
la var         FE0008         "        政府住房基金支出        "
la var         FE0007         "        地方财政税费附加        "
la var         FR01         "        收入总计        "
la var         FR201         "        补助收入        "
la var         FR301         "        上年结余        "
la var         FE01         "        支出总计        "
la var         FE301         "        上解支出        "
la var         FSP01         "        滚存结余        "
la var         VAT75         "        增值税75%        "
la var         Excise         "        消费税        "
la var         GDP         "        国内生产总值　        "
la var         Eco01         "        农业总产值　        "
la var         Eco02         "        工业总产值　        "
la var         POP         "        人口（万人)        "
la var         APOP         "        其中：农业人口        "
la var         Dep01         "        财政供养人口小计        "
la var         Dep02         "        在职人员        "
la var         Dep03         "        国家职工        "
la var         Dep04         "        集体职工        "
la var         Dep05         "        离休        "
la var         Dep06         "        退休        "
la var         Dep07         "        财政预算拨款及财政补助开支的年末人数        "
la var         Dep08         "        自收自支单位年末人数        "
la var         Inde         "                "

// 数据预处理完毕

// 以增值税为自变量，一般收入合计为因变量，进行回归分析；并令其经过原点（无常数项）
reg R001 R0002, nocons

//  R-squared 为可决系数，Root MSE 为残差对应的标准误 std. err.
// 当最小二乘法解决不了问题时，会采取 MLE GMM 或者 Kernel 等方法
// vce 为方差修正的方式

// // 绘散点图并拟合回归曲线
// gr twoway scatter R001 R0002, ///
//         || lfit R001 R0002

reg R001 R0002, r
/*修正误差项和X的标准差*/

reg R001 R0002

predict e, residual
/*生成了一个新的变量，预测前面估计值就是xb，预测残差就是residual*/

// scatter e R0002
// /*画残差对照X的散点图*/

// 对残差e升序排序
sort e

pnorm e

//
predict e, residual   //生成变量e作为残差
predict e, r // 提取残差
predict yhat, xb // 把x的系数算出来
scatter e R0002, xtitle(R0002) mc(red)  // 残差和R0002的关系，点改为红色
scatter e R0002 if R0002<=2000  //去掉大于2000的点
scatter R001 R0002 || lfit R001 R0002  散点图和直线拟合
reg R001 R0002, r 修正方差
corr R001 R0002 R0003 //计算几个变量的相关系数  结果是一个对称矩阵