位置: IT常识 - 正文

数学建模-回归分析(Stata)(数学建模回归模型例题)

编辑:rootadmin
数学建模-回归分析(Stata)

推荐整理分享数学建模-回归分析(Stata)(数学建模回归模型例题),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:数学建模回归模型论文,数学建模回归模型,数学建模回归分析案例,数学建模回归模型,数学建模回归分析案例,数学建模回归模型论文,数学建模回归分析,数学建模回归分析模型案例,内容如对您有帮助,希望把文章链接给更多的朋友!

 注意:代码文件仅供参考,一定不要直接用于自己的数模论文中国赛对于论文的查重要求非常严格,代码雷同也算作抄袭 如何修改代码避免查重的方法:https://www.bilibili.com/video/av59423231   //清风数学建模

 一、基础知识1.简介

X是自变量,Y是因变量。目的是通过X去预测Y。

一般处理模型像:期末成绩分析,Y是成绩,X是性别、是否是班干部、平时作业完成度等自变量。银行借贷成功率分析等问题。

2.不同数据类型的处理方法

 3.一元线性回归

①扰动项u与x均不相干,模型有外生性;否则存在内生性。包含了所有与y相关,但未添加到回归模型中的变量,如果这些变量和我们已经添加的自变量相关,则存在内生性。

②内生性的蒙特卡罗模拟:

相关系数绝对值越大,代表内生性越大(matlab代码)%% 蒙特卡洛模拟:内生性会造成回归系数的巨大误差times = 300; % 蒙特卡洛的次数R = zeros(times,1); % 用来储存扰动项u和x1的相关系数K = zeros(times,1); % 用来储存遗漏了x2之后,只用y对x1回归得到的回归系数for i = 1: times n = 30; % 样本数据量为n x1 = -10+rand(n,1)*20; % x1在-10和10上均匀分布,大小为30*1 u1 = normrnd(0,5,n,1) - rand(n,1); % 随机生成一组随机数 x2 = 0.3*x1 + u1; % x2与x1的相关性不确定, 因为我们设定了x2要加上u1这个随机数 % 这里的系数0.3我随便给的,没特殊的意义,你也可以改成其他的测试。 u = normrnd(0,1,n,1); % 扰动项u服从标准正态分布 y = 0.5 + 2 * x1 + 5 * x2 + u ; % 构造y k = (n*sum(x1.*y)-sum(x1)*sum(y))/(n*sum(x1.*x1)-sum(x1)*sum(x1)); % y = k*x1+b 回归估计出来的k K(i) = k; u = 5 * x2 + u; % 因为我们回归中忽略了5*x2,所以扰动项要加上5*x2 r = corrcoef(x1,u); % 2*2的相关系数矩阵 R(i) = r(2,1);endplot(R,K,'*')xlabel("x_1和u'的相关系数")ylabel("k的估计值")

③主要证明:核心变量与u不相干

核心解释变量:我们最感兴趣的变量,因此我们特别希望得到对其系数的一致估计(当样本容量无限增大时,收敛于待估计参数的真值 )。 控制变量:我们可能对于这些变量本身并无太大兴趣;而之所以把它们也放入回归方程,主要是为了 “控制住” 那些对被解释变量有影响的遗漏因素。

 ④什么时候取对数

对于什么时候取对数还没有固定的规则,但是有一些经验法则: (1)与市场价值相关的,例如,价格、销售额、工资等都可以取对数; (2)以年度量的变量,如受教育年限、工作经历等通常不取对数; (3)比例变量,如失业率、参与率等,两者均可; (4)变量取值必须是非负数,如果包含0,则可以对y取对数ln(1+y); 取对数的好处:(1)减弱数据的异方差性(2)如果变量本身不符合正态分布,取了对数后可能渐近服从正态分布(3)模型形式的需要,让模型具有经济学意义。

⑤三(四)种模型

1、一元线性回归:𝑦 = 𝑎 + 𝑏𝑥 + 𝜇,x每增加1个单位,y平均变化b个单位;2、双对数模型:𝑙𝑛𝑦 = 𝑎 + 𝑏𝑙𝑛𝑥 + 𝜇,x每增加1%,y平均变化b%;3、半对数模型:𝑦 = 𝑎+ 𝑏𝑙𝑛𝑥 + 𝜇,x每增加1%,y平均变化b/100个单位;4、半对数模型:𝑙𝑛𝑦= 𝑎+ 𝑏𝑥 + 𝜇,x每增加1个单位,y平均变化(100b)%。

 ⑥自变量有定性变量(如男女)引入虚拟变量x(女x=1 男x=0)

二、实战演练(Stata)1.题目

 打卡stata导入,勾上第一个

 2.定量数据和定性数据

定性数据是与描述有关的数据,可以观察到但无法计算。 相反,定量数据是一个专注于数字和数学计算,可以计算和计算的。

数学建模-回归分析(Stata)(数学建模回归模型例题)

比如数据中价格、评价量和重量都是可以拿数字直接表示,而产品名产地这种无法计算的只能定性处理。

summarize 团购价元 评价量 商品毛重kg定量:summarize 变量1 变量2 ... 变量n

评价变量,想评价啥就输入评价量

定性:tabulate 变量名,gen(A) 返回对应的这个变量的频率分布表,并生成对应的虚拟变量(以 A开头)。// 定性变量的频数分布,并得到相应字母开头的虚拟变量tabulate 配方,gen(A)tabulate 奶源产地 ,gen(B)tabulate 国产或进口 ,gen(C)tabulate 适用年龄岁 ,gen(D)tabulate 包装单位 ,gen(E)tabulate 分类 ,gen(F)tabulate 段位 ,gen(G)

 之后右侧会生成A1A2A3B1B2等变量

3.回归(默认OLS:普通最小二乘估计法)regress y x1 x2 … xkregress 评价量 团购价元 商品毛重kg

 加入虚拟变量回归

小技巧全选:鼠标在右边变量框下滑选中最后一个点击(文章是G4),再shift点击第一个即可全选,复制即可得到所有变量名。

regress 评价量 团购价元 商品毛重kg A1 A2 A3 B1 B2 B3 B4 B5 B6 B7 B8 B9 C1 C2 D1 D2 D3 D4 D5 E1 E2 E3 E4 F1  F2 G1 G2 G3 G4

注意只要定量变量+字母的定性变量

Stata会自动检测数据的完全多重共线性问题。

 上图的R-squared即R方,越高越好。

 4.标准化回归系数regress y x1 x2 … xk,beta

regress 评价量 团购价元 商品毛重kg,b

 5.拓展1.异方差;

横截面数据容易出现异方差的问题,时间序列数据容易出现自相关的问题。

异方差的检验:

在回归结束后运行命令: rvfplot      (画残差与拟合值的散点图) rvpplot x     (画残差与自变量x的散点图)2.拟合值出现负数的原因

 3.异方差的假设检验

BP法

Stata命令(在回归结束后使用): estat hettest ,rhs iid

原假设:扰动项不存在异方差 P值小于0.05,说明在95%的置信水平下拒绝原假设,即我们认为扰动项存在异方差怀特检验怀特检验原假设: 不存在异方差 Stata命令(在回归结束后使用): estat imtest,white

 4.处理异方差

使用OLS + 稳健的标准误

regress y x1 x2 … xk,robust

 5.检验多重共线性

(vif>10认为存在严重的多重共线性)

Stata计算各自变量VIF的命令(在回归结束后使用):estat vif

6. 逐步回归分析向前逐步回归Forward selection:stepwise regress y x1 x2 … xk, pe(#1)(显著才加入模型中).向后逐步回归Backward elimination:stepwise regress y x1 x2 … xk, pr(#2)(不显著就剔除出模型).stepwise reg 评价量 团购价元 商品毛重kg A1 A3 B1 B2 B3 B4 B5 B6 B7 B9 C1 D1 D2 D3 D4 E1 E2 E3 F1 G1 G2 G3, r pe(0.05)如果你觉得筛选后的变量仍很多,你可以减小#1或者#2   (#1 #2一般取0.05)如果你觉得筛选后的变量太少了,你可以增加#1或者#2 注: (1)x1 x2 … xk之间不能有完全多重共线性(和regress不同哦) (2)可以在后面再加参数b 和 r,即标准化回归系数或稳健标准误

一般向前向后逐步回归都跑一遍放论文里即可。

本文链接地址:https://www.jiuchutong.com/zhishi/298745.html 转载请保留说明!

上一篇:vue系列(三)——手把手教你搭建一个vue3管理后台基础模板(vue系列教程)

下一篇:【网络应用开发】实验1--Servlet技术及应用(网络应用开发技术)

  • 3dmax怎么把参考图放进去(3dmax怎么把参考图锁定)

    3dmax怎么把参考图放进去(3dmax怎么把参考图锁定)

  • 抖音直播怎么开始(抖音直播怎么开声音)

    抖音直播怎么开始(抖音直播怎么开声音)

  • 苹果微信如何关闭夜间模式(苹果微信如何关掉感应黑屏模式)

    苹果微信如何关闭夜间模式(苹果微信如何关掉感应黑屏模式)

  • 苹果手机搜不到蓝牙耳机(苹果手机搜不到家里的wifi)

    苹果手机搜不到蓝牙耳机(苹果手机搜不到家里的wifi)

  • 携程取消抢票会扣钱吗(携程取消抢票会怎么样)

    携程取消抢票会扣钱吗(携程取消抢票会怎么样)

  • 怎么通过快手号查电话号(怎么通过快手号找到微信号)

    怎么通过快手号查电话号(怎么通过快手号找到微信号)

  • 京东物流是京东自营吗(京东物流是京东快递吗)

    京东物流是京东自营吗(京东物流是京东快递吗)

  • qq号是什么(qq号是什么时候申请的怎么查)

    qq号是什么(qq号是什么时候申请的怎么查)

  • vivoy85可以单独换外屏吗(vivoy85能装几个手机卡)

    vivoy85可以单独换外屏吗(vivoy85能装几个手机卡)

  • 华为mate9pro上市时间(2021年华为mate9pro)

    华为mate9pro上市时间(2021年华为mate9pro)

  • 驱动精灵检测不到声卡(驱动精灵检测不到显卡)

    驱动精灵检测不到声卡(驱动精灵检测不到显卡)

  • 苹果手机老是自动删除联系人(苹果手机老是自动暗屏怎么弄)

    苹果手机老是自动删除联系人(苹果手机老是自动暗屏怎么弄)

  • iphone7plus有红外吗(iPhone7plus有红外功能吗)

    iphone7plus有红外吗(iPhone7plus有红外功能吗)

  • app未完成付款什么意思(app付款显示未完成付款)

    app未完成付款什么意思(app付款显示未完成付款)

  • 苹果7怎么屏蔽垃圾短信(苹果7怎么屏蔽手机通讯录联系人的信息)

    苹果7怎么屏蔽垃圾短信(苹果7怎么屏蔽手机通讯录联系人的信息)

  • 朋友圈翻译对方能看到吗(微信朋友圈翻译对方能看到吗)

    朋友圈翻译对方能看到吗(微信朋友圈翻译对方能看到吗)

  • 注销icloud账户的后果(注销icloud出错是怎么回事)

    注销icloud账户的后果(注销icloud出错是怎么回事)

  • 美团拼团成功后一个人用了另一个能退吗(美团拼团成功后一方取消订单)

    美团拼团成功后一个人用了另一个能退吗(美团拼团成功后一方取消订单)

  • 网页打开速度慢(网页打开速度慢的原因有哪些 手机)

    网页打开速度慢(网页打开速度慢的原因有哪些 手机)

  • airpods2买无线还是有线(airpods2无线充电需要单独买无线充电器吗)

    airpods2买无线还是有线(airpods2无线充电需要单独买无线充电器吗)

  • vivoiqoo在哪里设置返回键(vivo iqoo怎样)

    vivoiqoo在哪里设置返回键(vivo iqoo怎样)

  • tt语音怎么发布房间(tt语音怎么发布游戏房间)

    tt语音怎么发布房间(tt语音怎么发布游戏房间)

  • 小爱同学怎么调节音量大小(小爱同学怎么调音量)

    小爱同学怎么调节音量大小(小爱同学怎么调音量)

  • s1电池容量多大(五号电池多大)

    s1电池容量多大(五号电池多大)

  • 苹果手机话筒没声音(苹果手机话筒没声音重启又好了)

    苹果手机话筒没声音(苹果手机话筒没声音重启又好了)

  • 百度网盘如何加入好友(百度网盘如何加密)

    百度网盘如何加入好友(百度网盘如何加密)

  • 手机有信号但是打不出去电话(手机有信号但是没有4G怎么回事)

    手机有信号但是打不出去电话(手机有信号但是没有4G怎么回事)

  • Win11测试版Build 22000.184更新补丁KB5005642(附更新修复内容汇总+安装)(win11测试版和正式版区别)

    Win11测试版Build 22000.184更新补丁KB5005642(附更新修复内容汇总+安装)(win11测试版和正式版区别)

  • 在快速启动中不再显示最近使用的文件(要求快速启动)

    在快速启动中不再显示最近使用的文件(要求快速启动)

  • 怎样给笔记本的cpu降温(科普笔记本电脑怎样)

    怎样给笔记本的cpu降温(科普笔记本电脑怎样)

  • 奖金个人所得税税率
  • 公司之间借款如何做账
  • 出口退税进项税
  • 汇算清缴时房企毛利额有差异如何调整
  • 小规模纳税人开普票免税
  • 对公账户余额和实际不符
  • 其他收益是什么类
  • 可以先注销银行账户吗
  • 新会计准则 预付房租怎么做账
  • 结转完工产品生产成本会计科目
  • 银行回单卡是什么卡
  • 进口贴息对企业的好处
  • 购买商标权发生损失能税前扣除吗?
  • 捐赠利得会计分录
  • 公司属于亏损状态,但是法人个人收入算公司账吗
  • 材料采购定期结算表
  • 发票是不是都一样
  • 无法执行合同的说明函
  • 税务登记时要财务确认吗
  • 普通发票和专用发票的税点一样吗
  • 购车时服务费用怎么算
  • 房地产开发企业资质管理规定
  • 存货跌价准备怎么计提
  • 一个法人可以办几个公司
  • 小规模企业出口可以不报关吗
  • 小规模核定销售额是什么意思
  • 发票金额与实际付款金额不一致
  • 政府购买缴纳社保吗
  • 如何进行网速测试设置
  • 现金流量表与资产负债表不符
  • 公司车辆违章怎么办
  • 集团公司股份
  • 如何安全的处置电子邮件
  • Linux怎么在目录中创建文件
  • 进程间死锁
  • 委托加工物资加工费计入
  • 公司为实习生买保险列支什么科目?
  • 加利福尼亚riverside
  • phpif函数的使用方法
  • 卡特迈国家公园地图
  • phppoem
  • js表单验证代码
  • php.ini详解
  • 批发企业商品盘点发生短缺时
  • 纳税人数字签名怎么填
  • 淘宝镜像网站
  • 董事长报销应该是怎样的流程
  • 用税盘开发票税费扣除是怎么回事?
  • 财政收入的收缴具体包括
  • 财政拨付注册资本金说明
  • 销售部门发生的商品维修费
  • 工业企业的会计核算程序
  • 转账进公户
  • 28栏分次预缴税额是手动填写吗?
  • 母子公司往来款属于借款吗
  • 职工福利费计提标准
  • 固定资产如何入账处理
  • 建筑业异地预交税款最新规定
  • 员工医药费属于什么费用
  • 生产车间的保险有哪些
  • 差旅费报销单是谁填写的
  • 建筑企业外地预缴
  • 开了20万销项进项为0交多少税
  • 车辆保险证明怎么开
  • sql有什么
  • win10edge浏览器默认主页网址
  • xp系统英文版怎么改成中文
  • solaris命令大全
  • win8.1 升级
  • win10打开软件时卡住
  • windows7显示桌面的操作方法
  • linux无法运行yum
  • 逐一介绍班子成员表现
  • eclipse安转
  • 安卓自定义app
  • 高新技术企业认定
  • 广西空置房物业费减免政策吗
  • 纳税申报的期限是多少年
  • 深圳国税税发票在哪里开
  • 税率分几种 分别是什么
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设