数学建模多元线性回归内生性问题的蒙特卡洛模拟matlab代码
回归分析是数据分析中最基础也是最重要的分析工具,绝大多数的数据分析问题,都可以使用回归的思想来解决。回归分析的任务就是,通过研究自变量X和因变量Y的相关关系,尝试去解释Y的形成机制,进而达到通过X去预测Y的目的。
常见的回归分析有五类:线性回归、0-1回归、定序回归、计数回归和生存回归,其划分的依据是因变量Y的类型。本讲我们主要学习线性回归。
%% 蒙特卡洛模拟:内生性会造成回归系数的巨大误差
times = 300; % 蒙特卡洛的次数
R = zeros(times,1); % 用来储存扰动项u和x1的相关系数
K = zeros(times,1); % 用来储存遗漏了x2之后,只用y对x1回归得到的回归系数
for i = 1: timesn = 30; % 样本数据量为nx1 = -10+rand(n,1)*20; % x1在-10和10上均匀分布,大小为30*1u1 = normrnd(0,5,n,1) - rand(n,1); % 随机生成一组随机数x2 = 0.3*x1 + u1; % x2与x1的相关性不确定, 因为我们设定了x2要加上u1这个随机数u = normrnd(0,1,n,1); % 扰动项u服从标准正态分布y = 0.5 + 2 * x1 + 5 * x2 + u ; % 构造yk = (n*sum(x1.*y)-sum(x1)*sum(y))/(n*sum(x1.*x1)-sum(x1)*sum(x1)); % y = k*x1+b 回归估计出来的kK(i) = k;u = x2 + u; % 因为我们忽略了x2,所以扰动项要加上x2r = corrcoef(x1,u); % 2*2的相关系数矩阵R(i) = r(2,1);
end
plot(R,K,'*')
xlabel("x_1和u'的相关系数")
ylabel("k的估计值")% % 注意:代码文件仅供参考,一定不要直接用于自己的数模论文中
% % 国赛对于论文的查重要求非常严格,代码雷同也算作抄袭
% % 更多优质数模资料可在我的微店获取:https://weidian.com/?userid=1372657210
% % 数学建模讨论群获取地址:http://note.youdao.com/noteshare?id=4997251d8219a45d56631e412b1e9392
内生性将会导致回归系数的不一致哦~
数学建模多元线性回归内生性问题的蒙特卡洛模拟matlab代码相关推荐
- 数学建模——多元线性回归 /多元线性拟合 (工具:matlab)
文章目录 一.多元线性回归是什么 二.样例及代码 第一步:画散点图进行分析 第二步:拟合 三.总结: 四.参考附录: 技巧篇
- 【数学建模】常用微分方程模型 + 详细手写公式推导 + Matlab代码实现
文章目录 一.学习内容 二.学习时间 三.学习产出 3.1 微分方程基本概念 3.2 微分方程在数学建模中的应用 3.3 微分方程常用模型 3.3.1 人口增长模型 3.3.1.1 指数增长模型(马尔 ...
- 【数学建模】单、多因素试验的方差分析(Matlab代码实现)
目录 1 概述 2 单因素方差分析 2.1 语法 2.2 算例 1 2.2.1 算例 2.2.2 Matlab代码 2.2.3 结果 2.3 算例2 2.3.1 算例 2.3.2 Matlab代码 2 ...
- 数学建模——一维、二维插值模型详解Python代码
数学建模--一维.二维插值模型详解Python代码 一.一维插值 # -*-coding:utf-8 -*- import numpy as np from scipy import interpol ...
- 2021年数学建模国赛C题问题三详细思路和代码
2021年数学建模国赛C题问题三详细思路和代码 1 问题分析 问题三要求压缩成本.减少转运及仓储的成本.降低转运损耗率,并制定对应的订购方案和转运方案,最大的特点在于没有供应商数量限制,故每周的供应商 ...
- 2021年数学建模国赛C题问题二详细思路和代码
2021年数学建模国赛C题问题一详细思路和代码 说明: 本思路仅为小编个人思路,其中关于本文中给出的AdaBoost回归预测模型于今年的 国赛参考评分标准有所出入,所以本题思路为个人创新思路,不一定符 ...
- 数学建模——智能优化之粒子群模型详解Python代码
数学建模--智能优化之粒子群模型详解Python代码 import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplo ...
- 【2022年华为杯数学建模E题赛后总结加思路详细介绍配代码----10月11号写的总结】
提示:下文将介绍2022年华为杯数学建模E题赛后总结加思路详细介绍配代码 傻逼队友,傻逼队友,傻逼队友一定要看好人在进行组队,这是劝告. 这里有几点总结进行描述: 第一,图一定要尽量多,对图的解释要多 ...
- 数学建模、运筹学之整数规划(原理、例题、代码)
数学建模.运筹学之整数规划(原理.例题.代码) 一.何为整数规划? 二.整数规划与线性规划 三.分枝定界法 四.0-1整数规划 一.何为整数规划? 整数规划具有深远的实际应用价值,因为现实生活中的规划 ...
最新文章
- 高考631能上什么好的计算机学校,2021年高考630分能上什么大学 可以报哪些学校...
- SunDay天气——开放源代码
- VTK:IO之GenericDataObjectReader
- 无向图:查找最小环集合(最短路径回溯算法)
- 时隙aloha协议仿真程序_工控ModbusTCP/IP协议仿真环境搭建
- 最大化 AIX 上的 Java 性能,第 3 部分: 更多就是更好
- Windows Server 2008 配置使用动态IP和备用地址
- 使用拷贝的方式(adb push) 绕过Android系统和adb install直接安装APK
- python Web开发你要理解的WSGI uwsgi详解
- 结网读书笔记-从产品经理的角度看产品
- 智力问答选择题_2018智力问答题大全及答案:智力问答题大全及答案
- 【官方】下载最新adb及安装驱动的方法
- 2019最新python下载安装详解
- 大家崇拜凯文.米特尼克吗?
- HDU 6638 Snowy Smile(线段树)
- 嵌入式音频架构 - AudioWeaver模块开发
- RFSoC应用笔记 - RF数据转换器 -09- RFSoC关键配置之RF-DAC内部解析(三)
- Win10 + GTX1650显卡 安装 Cuda-10.1 + CuDNN-7.6 + Tensorflow-gpu-1.15.0 一步到位
- 不要混淆事件的互不相容和互相独立
- 计算机windows无效映像,应用程序或dll 为无效的windows映像怎么办
热门文章
- 拓扑关系——九交模型
- delphi 控件的安装
- 如何制作生成gif?快速生成制作生成的方法来了
- Python软件目录结构规范
- 一名前端开发者的2020年终总结
- android 平台的RNDIS功能
- ShareLatex+Overflow:PDF Rendering Error Something went wrong while rendering this PDF问题解决
- matlab R2011a汉化包,MATLAB R2011a (7.12)发布了,MATLAB R2011a最新功能
- C# windows media player 双窗口同步播放视频
- 仙人掌之歌——投石问路(2)