研究生数理统计多元线性回归及显著性检验Matlab程序完美版.doc

《研究生数理统计多元线性回归及显著性检验Matlab程序完美版.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、多元线性回归及显著性检验Matlab程序（完美版）一、说明：1、本程序是研究生教材《数理统计》（杨虎、刘琼、钟波编著）例4.4.1（P133）的Matlab编程解答程序。教材上的例题只做了回归方程显著性分析和一次回归系数显著性分析（剔除x1后没有再检验x2和x3）。2、本程序在以上的基础之上，还分别检验了x2和x3，并且计算精度更高。3、本程序可根据用户的需要，在输入不同的显著性水平α之下得到相应的解答。4、本程序移植性强，对于其他数据，只需要改变excel中的数据即可。5、本程序输出的可读性强，整洁美观。二、数据入下（将数据存入excel表

2、格，文件名为jc_p133_example.xls。注意数据是按x1，x2，…，xk，y这样来列来存储。若不是3个变量，则相应增减数据列就行。）：218504.330279403.6485514464.483123435.5468120645.497312403.1125317645.118265393.875978374.67023554.9536316605.006018495.270184505.3772614515.4849021514.596314515.6645712566.0795160483.2194616455.807601

3、5524.730690404.680546323.2172017472.610490443.7174216393.894696392.7066125515.6314613415.8152127475.1302024615.391512374.4533415494.6569020454.5212616424.865417485.3566104484.6098414362.3815513363.874698514.5919613545.1588581005.4373511443.99686634.397213554.062278502.29054

4、10454.7115105404.531317645.3637415726.0771三、完整程序如下：%----------------------------byggihhimm----------------------------%《数理统计》杨虎、刘琼、钟波编著例4.4.1多元线性回归及显著性检验完整解答%输入需要的显著水平α（默认α=0.02），计算出不同结果（见运行结果）%该程序也适合其他维数的数据分析（只需改变excel表格中的数据即可）%----------------------------byggihhimm-------

5、---------------------clear;clc;data=xlsread('jc_p133_example.xls','sheet1');xi=data(:,1:end-1);[n,k]=size(data);k=k-1;index_of_xi_array=ones(1,k);X=[ones(n,1)xi];Y=data(:,end);fprintf('第1次计算结果：r')beta_mao=((X'*X)X'*Y)';fmt_str0='';fori0=1:k+1fmt_str0=[fmt_str0'β'num2str(i

6、0-1)'=%0.4fr'];endfprintf(fmt_str0,beta_mao)fprintf('r')%%检验回归方程的显著性x_ba=mean(xi);y_ba=mean(Y);St_square=sum(Y.^2)-n*y_ba^2;lxy=sum((xi-ones(n,1)*x_ba).*((Y-y_ba)*ones(1,k)));Sr_square=sum(beta_mao(2:end).*lxy);Se_square=St_square-Sr_square;c_flag=Sr_square/Se_square;F_al

7、pha=input('>>>>>>请输入您要求的显著性水平(0<α<1)α=');while~(isscalar(F_alpha)&&F_alpha<1&&F_alpha>0)F_alpha=input('您的输入有误，请重新输入一个大于0，小于1的数，α=');endF_fenweidian=finv(1-F_alpha,k,n-k-1);c=k/(n-k-1)*F_fenweidian;ifc_flag>cfprintf(['r--------------------回归方程显著性检验（H0：β1=β2=...=βk=0）'...'--

8、------------------r经过计算：拒绝H0，原假设不成立。'])elsefprintf(['r--------------------回归方程显著