【原创】R语言模型中的加总偏误与内生性 一种数值模拟方法数据分析报告论文(含代码数据)
- 格式:docx
- 大小:155.00 KB
- 文档页数:4
咨询QQ:3025393450
有问题百度搜索“”就可以了
欢迎登陆官网:/datablog
R语言模型中的加总偏误与内生性:一种数值模拟方法数据分析报告
来源:大数据部落| 有问题百度一下“”就可以了
引言本文中主题是内生性,它可能严重偏向回归估计。
我将专门模拟由遗漏变量引起的内生性。
在本系列的后续文章中,我将模拟其他规范问题,如异方差性,多重共线性和对撞机偏差。
数据生成过程
考虑一些结果变量的数据生成过程(DGP)
1 2 3 4五6 7 8 9# simulation parameters
set.seed(144);
a=50; b=.5; c=.01;
x=rnorm(n=ss,mean=1000,sd=50); z=d+h*x+rnorm(ss,0,10)
咨询QQ:3025393450
有问题百度搜索“”就可以了
欢迎登陆官网:/datablog
模拟
模拟将估计下面的两个模型。
第一个模型是正确的,它包含实际DGP中的所有术语。
但是,第二个模型省略了DGP中存在的变量。
相反,变量被误入了误差
1
2 3 4 6 7 8 9
1 0
1 1
1 2sim=function(endog){
e=rnorm(n=ss,mean=0,sd=10)
# Select data generation process
if(endog==TRUE){ fit lm(y~x) }else{ fit=lm(y~x+z)}
return(fit$coefficients)
}
sim_results_endog=t(replicate(trials,sim(endog=TRUE)))
略,仿真得到的红色采样分布,围绕0.5895居中。
它偏离.5895的真实值。
此外,偏差采样分布的方差远小于周围的真实方差能力。
咨询QQ:3025393450
有问题百度搜索“”就可以了
欢迎登陆官网:/datablog
造正相关并且它们的斜率系数是正的,所以偏差将是正的。
根据模拟的参数,应以.0895为中心,非常接近真实的偏差值。
咨询QQ:3025393450
有问题百度搜索“”就可以了
欢迎登陆官网:/datablog
上述推导还可以让我们确定从知道的相关偏差的方向面的案例很一般,但有特殊的应用。
例如,如果我们认为个人的收入是教育年限和工作年经验的函数,那么省略一个变量将偏向另一个变量的斜率估计。