- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
回归分析的模型
线性回归模型
假设1 解释变量是非随机的或固定的,且各之间互不相关。
假设2 随机误差项具有零均值、同方差及不序列相关性。
假设3 解释变量与随机项不相关。
假设4 随机变量满足正态分布。
模型的建立与求解
回归方程的建立
设是一个可观测的随机变量,它受到个非随机因素和随机误差的影响。且有如下线性关系:
且
要估计未知参数(),需要进行次观测,得到下式
, ,
的最小二乘估计
为的一个无偏估计。
则最终的回归方程为
线性回归关系的检验
(1)建立方差分析表
离差平方和 ,
残差平方和
回归平方和
假设
其中为的观测值
检验的影响是否显著
的置信度为的置信区间为
同理的置信度为的置信区间为
残差分析原模型的合理性
对残差作分析,用Matlab做出残差图,观察观测值是否有异常,
程序:
x1=[53 47 57 61 55 56 51 58 64 61 74 62 59 50 54 63 57 56 54 55 54 72 57 58 53 85 72 63 47 44];
x2=[89 83 80 92 104 97 97 123 111 111 115 109 111 110 112 109 89 95 94 87 82 92 97 99 101 120 121 110 102 85];
x3=[76 88 51 81 96 99 121 157 127 159 145 143 131 136 124 159 145 137 111 91 82 116 119 112 106 156 236 149 120 96];
x4=[19 29 31 28 34 30 31 54 47 51 65 59 45 44 44 60 61 48 39 41 53 75 66 50 32 52 73 57 40 40 ];
x5=[30 8 13 8 8 10 28 9 8 24 4 27 38 27 9 5 15 36 26 12 11 8 24 32 43 19 18 30 62 23];
y=[90 143 58 142 175 215 250 309 273 329 299 299 246 261 260 295 282 262 204 179 227 277 242 226 173 266 426 307 230 201];
X=[ones(length(y),1),x1,x2,x3,x4,x5];
Y=y;
[b,bint,r,rint,stats]=regress(Y,X);
b,bint,stats
运行结果为:
b =
-12.2923
-1.5271
0.8327
1.3775
2.0068
-0.3086
bint =
-106.3464 81.7618
-2.9355 -0.1186
-0.4217 2.0871
0.8554 1.8995
1.0254 2.9882
-1.0032 0.3860
stats =
0.9284 62.2206 0.0000 474.4773
因此可得出,
r2=0.9284,F=62.2206,p=0.0000 则p0.05,回归模型为:
在Matlab命令中输入
rcoplot(r,rint)
可得到残差图如下图所示
从残差图中和以上分析可以看出,此回归方程效果良好
作图
Z=
按同样步骤对全部数据进行了回归分析,运行后的结果为
b =
-32.4586
0.1718
0.4463
0.8737
2.2106
-0.3352
bint =
-46.4503 -18.4669
-0.2300 0.5736
0.1428 0.7499
0.6943 1.0531
1.7915 2.6297
-0.4650 -0.2053
stats =
0.8373 238.7281 0 927.5409
同样可得出,
r2=0.8373,F=238.7281,p=0.0000 则p0.05,回归模型为:
在Matlab命令中输入
coplot(r,rint)
可得到残差图如下图所示
从残差图中和以上理论分析可以看出,此回归方程效果也亦良好
文档评论(0)