R软件一元线性回归分析(非常详细)

R软件一元线性回归分析(非常详细)

ID:41402464

大小:62.51 KB

页数:8页

时间:2019-08-24

R软件一元线性回归分析(非常详细)_第1页
R软件一元线性回归分析(非常详细)_第2页
R软件一元线性回归分析(非常详细)_第3页
R软件一元线性回归分析(非常详细)_第4页
R软件一元线性回归分析(非常详细)_第5页
资源描述:

《R软件一元线性回归分析(非常详细)》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、R软件一元线性回归分析合金钢强度与碳含量的数据序号碳含量/%合金钢强度/107pa10.1042.020.1143.030.1245.040.1345.050.1445.060.1547.570.1649.080.1753.090.1850.0100.2055.0110.2155.0120.2360.0这里取碳含量为x是普通变量,取合金钢强度为y是随机变量使用R软件对以上数据绘出散点图程序如下:>x=matrix(c(0.1,42,0.11,43,0.12,45,0.13,45,0.14,45,0.

2、15,47.5,0.16,49,0.17,53,0.18,50,0.2,55,0.21,55,0.23,60),nrow=12,ncol=2,byrow=T,dimnames=list(1:12,c("C","E")))>outputcost=as.data.frame(x)>plot(outputcost$C,outputcost$E)很显然这些点基本上(但并不精确地)落在一条直线上。下面在之前数据录入的基础上做回归分析(程序接前文,下同)>lm.sol=lm(E~C,data=outputcos

3、t)>summary(lm.sol)得到以下结果:Call:lm(formula=E~C,data=outputcost)Residuals:Min1QMedian3QMax-2.00449-0.63600-0.024010.712972.32451Coefficients:EstimateStd.ErrortvaluePr(>

4、t

5、)(Intercept)28.0831.56717.926.27e-09***C132.8999.60613.847.59e-08***---Signif.codes:

6、0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’1Residualstandarderror:1.309on10degreesoffreedomMultipleR-squared:0.9503,AdjustedR-squared:0.9454F-statistic:191.4on1and10DF,p-value:7.585e-08由计算结果分析:常数项=28.083,变量(即碳含量)的系数=132.899得到回归方程:=28.083+132.899x由于回归模型建立使用的是最小二

7、乘法,而最小二乘法只是一种单纯的数学方法,存在着一定的缺陷,即不论变量间有无相关关系或有无显著线性相关关系,用最小二乘法都可以找到一条直线去拟合变量间关系。所以回归模型建立之后,还要对其进行显著性检验:在上面的结果中sd()=1.567,sd()=9.606。而对应于两个系数的P值6.27e-09和7.59e-08,故是非常显著的。关于方程的检验,残差的标准差=1.309。相关系数的平方R=0.9503。关于F分布的P值为7.585e-08,也是非常显著的。我们将得到的直线方程画在散点图上,程序如下

8、:>abline(lm.sol)得到散点图及相应的回归直线:下面分析残差:在R软件中,可用函数residuals()计算回归方程的残差。程序如下:>y.res=residuals(lm.sol);plot(y.res)得到残差图从残差图可以看出,第8个点有些反常,这样我们用程序将第8个点的残差标出,程序如下:>text(8,y.res[8],labels=8,adj=1.2)这个点可能有问题,下面做简单处理,去掉该样本点,编程如下:>i=1:12;outputcost2=as.data.frame(

9、x[i!=8,])lm2=lm(E~C,data=outputcost2)summary(lm2)结果输出如下:Call:lm(formula=E~C,data=outputcost2)Residuals:Min1QMedian3QMax-1.7567-0.5067-0.13080.68211.6787Coefficients:EstimateStd.ErrortvaluePr(>

10、t

11、)(Intercept)28.1241.33521.065.75e-09***C131.2938.21715.98

12、6.51e-08***---Signif.codes:0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’1Residualstandarderror:1.115on9degreesoffreedomMultipleR-squared:0.966,AdjustedR-squared:0.9622F-statistic:255.3on1and9DF,p-value:6.506e-08由结果分析,去掉第8个点之后,回归方程系数变化不大,R2相关系数有所提高

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。