R语言的BP检验是否有效？
发布：szc412247754 | 分类：R语言培训
关于本站
人大经济论坛-经管之家：分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!
获取电子版《CDA一级教材》
完整电子版已上线CDA网校，累计已有10万+在读~ 教材严格按考试大纲编写，适合CDA考生备考，也适合业务及数据分析岗位的从业者提升自我。
TOP热门关键词
《计量经济学（第二版）》庞浩主编第五章案例分析在用R语言做这个例题时，发现bptest()对加权回归模型的检验无效，哪位大神能够帮一下忙，告诉我如何做加权回归模型的White检验？？？而且加权回归所得到的结果与Evie ...
免费学术公开课,扫码加入 |
《计量经济学（第二版）》 庞浩主编 第五章 案例分析
在用R语言做这个例题时，发现bptest()对加权回归模型的检验无效，哪位大神能够帮一下忙，告诉我如何做加权回归模型的White检验？？？而且加权回归所得到的结果与Eviews做出来的结果不同。。。
| 地区 | 人口数（万人） | 医疗机构数（个）Y | 地区 | 人口数（万人） | 医疗机构数（个）Y |
| X | X | ||||
| 成都 | 1013.3 | 6304 | 眉山 | 339.9 | 827 |
| 自贡 | 315 | 911 | 宜宾 | 508.5 | 1530 |
| 攀枝花 | 103 | 934 | 广安 | 438.6 | 1589 |
| 泸州 | 463.7 | 1297 | 达州 | 620.1 | 2403 |
| 德阳 | 379.3 | 1085 | 雅安 | 149.8 | 866 |
| 绵阳 | 518.4 | 1616 | 巴中 | 346.7 | 1223 |
| 广元 | 302.6 | 1021 | 资阳 | 488.4 | 1361 |
| 遂宁 | 371 | 1375 | 阿坝 | 82.9 | 536 |
| 内江 | 419.9 | 1212 | 甘孜 | 88.9 | 594 |
| 乐山 | 345.9 | 1132 | 凉山 | 402.4 | 1471 |
| 南充 | 709.2 | 4064 |
####普通回归
>rd<-read.delim("P142.txt")
>lmsol<-lm(Y~X,data=rd)
> summary(lmsol)
Call:
lm(formula = Y ~ X, data = rd)
Residuals:
Min 1Q Median 3Q Max
-700.2-436.3 -163.6 624.11422.6
Coefficients:
EstimateStd. Errort valuePr(>|t|)
(Intercept) -562.9074291.5642-1.931 0.0686 .
X 5.3728 0.6442 8.3408.99e-08***
---
Signif. codes:0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’1
Residual standard error: 623 on 19 degreesof freedom
Multiple R-squared: 0.7854, Adjusted R-squared: 0.7741
F-statistic: 69.55 on 1 and 19 DF,p-value: 8.986e-08
####White检验
>bptest(lmsol,varformula=~X+I(X^2),data=rd)
studentized Breusch-Pagan test
data: lmsol
BP = 18.0748, df = 2,p-value = 0.0001189 # 存在异方差
####加权回归，权数为：weights=1/ (rd$ X)^2
>lmsol2<-lm(Y~X,data=rd,weights=1/ (rd$ X)^2)
> summary(lmsol2)
Call:
lm(formula = Y ~ X, data = rd, weights =1/(rd$X)^2)
Residuals:
Min 1Q Median 3Q Max
-1.42205-0.75313 -0.368250.02241 3.11812
Coefficients:Estimate Std. Error tvalue Pr(>|t|)
(Intercept) 384.6123 87.90444.375 0.000325 ***
X 2.7236 0.43346.284 0.00000494 ***
---
Signif. codes:0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’1
Residual standard error: 1.273 on 19degrees of freedom
Multiple R-squared: 0.6752, AdjustedR-squared: 0.6581
F-statistic: 39.49 on 1 and 19 DF,p-value: 0.000004935
# 回归模型系数和可决系数与书上的结果不同。我自己用Eviews7做了一次，得到的结果与书上的结果不同大致相同。
####加权回归后White检验
>bptest(lmsol2,varformula=~X+I(X^2),data=rd)
studentized Breusch-Pagan test
data: lmsol2
BP = 18.0748, df = 2,p-value = 0.0001189
# 加权回归后White检验与未加权回归White检验结果相同，为何？
「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
您可能感兴趣的文章
人气文章
1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责；
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性，不作出任何保证或承若；
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
