- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
应用回归分析·上机作业三
学号: 200930980106 姓名: 何斌 年级专业: 10 级统计 1 班 指导老师:丁仕虹
1. 用主成分分析建模分析习题 5.9
执行程序如下:
proc import out=xt59 /* 使用 import 过程导入数据并输出到数据集 xt59*/
datafile=d:\xt5.9.xls
dbms=excel replace;
getnames=yes; /* 首行为变量名 */
run;
proc princomp data=xt59 out=ppp;/* 运行主成分分析过程 */
var x1-x6;/* 分析 x1-x6 的主成分 */
run;
得到结果如下:
图 1-1 相关系数矩阵的特征根及对应的特征向量 ( 主成分 )
图 1-1 中有 6 个主成分的特征值( Eigenvalues) ,最大的是 ????=5.1325 ,最小的是
????=0.0002 。由于????=0.0002 近似接近于 0 或者由于条件数 k= ????/ ????10,所以可以认
为自变量之间存在共线性。
两个主成分的方差贡献率已达到了 97.11%,故选择两个主成分进行回归分析。
执行程序如下:
proc reg data=ppp;
model y=prin1 prin2;
run;
得到结果如下:
图 1-2 主成分回归模型
从图 1-2 可以看出,回归模型通过了显著性检验, Adj R-Sq=0.9784 拟合优度高,
回归系数均通过显著性检验, 用 y 对前两个主成分做普通最小二乘回归的回归方程为:
??= ????????. ??????????+ ????????.????????????????????- ??????. ????????????????????(1-1 )
由
??????????= ??. ??????????????+ ??. ????????????????+ ??. ????????????????+ ??. ????????????????+ ??. ??????????????+ ??. ????????????????
??????????= -?? . ??????????????- ??. ????????????????- ??.????????????????+ ??. ????????????????- ??.??????????+ ??.????????????????
将 ??????????和??????????代入( 1-1 ),可得回归方程为:
??= ????????.??????????- ??.??????????????- ??.??????????????- ??.??????????????- ??. ??????????????+ ??. ??????????????- ??.??????????????
2. 用岭回归分析建模分析习题 5.9
执行程序如下:
proc reg data=xt59 outest=rid;
model y=x1-x6/ridge=0.1 to 1 by 0.1;
/* ridge=0.1 to 1 by 0.1 要求分别取岭回归的 k 为 0.1 ,0.2 ,…, 0.9 ,1 建立岭
回归方程 */
plot/ridgeplot; /* 绘制岭迹图以选择 k 值*/
run;
得到结果如下:
y = 1348.2 -0.641 x1 -0.317 x2 -0.4127 x3 -0.0021 x4 +0.6711 x5 -0.0075 x 6
0.40
文档评论(0)