资源描述:
《用matlab做一元线性回归分析.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、用matlab做一元线性回归分析一元线性回归分析是在排除其他影响因素的假定其他影响因素确定的情况下,分析某一个因素(自变量)是如何影响另外一个事物(因变量)的过程,所进行的分析是比较理想化的。用SPSS可以做一元线性回归分析,但是当回归的自变量比较多的时候,一个一个的输入会比较麻烦,增加了计算量,本文中描述了如何用matlab语言来实现一元线性回归分析。在matlab中,regress命令是用来做回归的。假如有96个SNP,作为自变量,有一个因变量,比如说HDL,LDL等等,将它们以列导入matlab。值得注意的是:自变量前面必须有一列全为1的数据,看下面
2、例子即可理解。 fori=1:96z=[ones(2334,1),x(:,i)];[b,bint,r,rint,stats]=regress(y,z);c(i,:)=stats;end在一元线性回归方程中,回归方程的显著性检验可以替代回归系数的显著性检验,并且F=T2 百度中的一个例子: X=[1146811141721]'Y=[2.493.303.6812.2027.0461.10108.80170.90275.50]'X=[ones(9,1),X][b,bint,r,rint,stats]=regress(Y,X)输出向量b,bint为回归系数估
3、计值和它们的置信区间,r,rint为残差及其置信区间,stats是用于检验回归模型的统计量,有三个数值,第一个是R2,其中R是相关系数,第二个是F统计量值,第三个是与统计量F对应的概率P,当P<α时拒绝H0,回归模型成立 regressMultiplelinearregressionSyntaxb=regress(y,X)[b,bint]=regress(y,X)[b,bint,r]=regress(y,X)[b,bint,r,rint]=regress(y,X)[b,bint,r,rint,stats]=regress(y,X)[...]=regress
4、(y,X,alpha)Descriptionb=regress(y,X)returnstheleastsquaresfitofyonXbysolvingthelinearmodel forβ,where:yisann-by-1vectorofobservationsXisann-by-pmatrixofregressorsβisap-by-1vectorofparametersɛisann-by-1vectorofrandomdisturbances[b,bint]=regress(y,X)returnsamatrixbintof95%confidence
5、intervalsforβ.[b,bint,r]=regress(y,X)returnsavector,rofresiduals.[b,bint,r,rint]=regress(y,X)returnsamatrixrintofintervalsthatcanbeusedtodiagnoseoutliers.Ifrint(i,:)doesnotcontainzero,thentheithresidualislargerthanwouldbeexpected,atthe5%significancelevel.Thisisevidencethattheithob
6、servationisanoutlier. [b,bint,r,rint,stats]=regress(y,X)returnsavectorstatsthatcontainstheR2statistic,theFstatisticandapvalueforthefullmodel,andanestimateoftheerrorvariance.