欢迎来到淘文阁 - 分享文档赚钱的网站! | 帮助中心 好文档才是您的得力助手!
淘文阁 - 分享文档赚钱的网站
全部分类
  • 研究报告>
  • 管理文献>
  • 标准材料>
  • 技术资料>
  • 教育专区>
  • 应用文书>
  • 生活休闲>
  • 考试试题>
  • pptx模板>
  • 工商注册>
  • 期刊短文>
  • 图片设计>
  • ImageVerifierCode 换一换

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

    • 资源ID:29935672       资源大小:61KB        全文页数:8页
    • 资源格式: DOC        下载积分:15金币
    快捷下载 游客一键下载
    会员登录下载
    微信登录下载
    三方登录下载: 微信开放平台登录   QQ登录  
    二维码
    微信扫一扫登录
    下载资源需要15金币
    邮箱/手机:
    温馨提示:
    快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。
    如填写123,账号就是123,密码也是123。
    支付方式: 支付宝    微信支付   
    验证码:   换一换

     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

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

    R软件一元线性回归分析合金钢强度与碳含量的数据序号碳含量/%合金钢强度/107pa10.10 42.0 20.11 43.0 30.12 45.0 40.13 45.0 50.14 45.0 60.15 47.5 70.16 49.0 80.17 53.0 90.18 50.0 100.20 55.0 110.21 55.0 120.23 60.0 这里取碳含量为x是普通变量,取合金钢强度为y是随机变量使用R软件对以上数据绘出散点图程序如下:> x=matrix(c(0.1,42,0.11,43,0.12,45,0.13,45,0.14,45,0.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(EC,data = outputcost)>summary(lm.sol)得到以下结果:Call:lm(formula = E C, data = outputcost)Residuals: Min 1Q Median 3Q Max -2.00449 -0.63600 -0.02401 0.71297 2.32451 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 28.083 1.567 17.92 6.27e-09 *C 132.899 9.606 13.84 7.59e-08 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 Residual standard error: 1.309 on 10 degrees of freedomMultiple R-squared: 0.9503, Adjusted R-squared: 0.9454 F-statistic: 191.4 on 1 and 10 DF, p-value: 7.585e-08由计算结果分析:常数项=28.083,变量(即碳含量)的系数=132.899得到回归方程:=28.083+132.899x由于回归模型建立使用的是最小二乘法 ,而最小二乘法只是一种单纯的数学方法 ,存在着一定的缺陷 ,即不论变量间有无相关关系或有无显著线性相关关系 ,用最小二乘法都可以找到一条直线去拟合变量间关系。所以回归模型建立之后 ,还要对其进行显著性检验 :在上面的结果中sd()=1.567,sd()=9.606。而对应于两个系数的P值6.27e-09和7.59e-08,故是非常显著的。关于方程的检验,残差的标准差=1.309。相关系数的平方R = 0.9503。关于F分布的P值为7.585e-08,也是非常显著的。我们将得到的直线方程画在散点图上,程序如下:> abline(lm.sol)得到散点图及相应的回归直线:下面分析残差:在R软件中,可用函数residuals()计算回归方程的残差。程序如下:> y.res=residuals(lm.sol);plot(y.res)得到残差图从残差图可以看出,第8个点有些反常,这样我们用程序将第8个点的残差标出,程序如下:>text(8,y.res8,labels=8,adj=1.2)这个点可能有问题,下面做简单处理,去掉该样本点,编程如下:>i=1:12;outputcost2=as.data.frame(xi!=8,)lm2=lm(EC,data=outputcost2)summary(lm2)结果输出如下:Call:lm(formula = E C, data = outputcost2)Residuals: Min 1Q Median 3Q Max -1.7567 -0.5067 -0.1308 0.6821 1.6787 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 28.124 1.335 21.06 5.75e-09 *C 131.293 8.217 15.98 6.51e-08 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 Residual standard error: 1.115 on 9 degrees of freedomMultiple R-squared: 0.966, Adjusted R-squared: 0.9622 F-statistic: 255.3 on 1 and 9 DF, p-value: 6.506e-08由结果分析,去掉第8个点之后,回归方程系数变化不大,R2 相关系数有所提高,并且p-值变小了,这说明样本点8可以去掉。所得新模型较为理想。总结程序如下:> x2=matrix(c(0.1,42,0.11,43,0.12,45,0.13,45,0.14,45,0.15,47.5,0.16,49 ,0.18,50,0.2,55,0.21,55,0.23,60),nrow=11,ncol=2,byrow=T,dimnames=list(1:11,c("C","E") >outputcost=as.data.frame(x2)>plot(outputcost$C,outputcost$E)>lm.sol = lm(EC,data = outputcost)>summary(lm.sol)Call:lm(formula = E C, data = outputcost)Residuals: Min 1Q Median 3Q Max -1.7567 -0.5067 -0.1308 0.6821 1.6787 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 28.124 1.335 21.06 5.75e-09 *C 131.293 8.217 15.98 6.51e-08 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 Residual standard error: 1.115 on 9 degrees of freedomMultiple R-squared: 0.966, Adjusted R-squared: 0.9622 F-statistic: 255.3 on 1 and 9 DF, p-value: 6.506e-08>abline(lm.sol)得到最后的散点图和回归直线得到回归方程:=28.124+131.293x8回 分 .程程直归和的 0- , : . 0 : 0. 0 0*0 * -* 0 . . . * . . > . 0 . 0. , = .( = . $ $ . )", ( = 0 . 0 0 . 0 ,0 , . 0 , 0( 下序想为较得掉 样这,-并,数相 不数系,点个去分0- , : 0 , . - : . 0*.0 * .-* - . *0 . . ) ( . . . = 下出 =( = :下下点该,单面,能 . , 下如出差个第程样常有点出以差 () ( .下序程程归算 函,件差残线直应图 下序,点画方到的著常 0 .的布。 0 平系关0=准差验的显是 .和-.值系于而0 ( 中 检显进要 立归所关量合条到都法用 关著有关无间论 的着,法的种只二最,二是立模 .=方 . 数量即量 0 析果0 0 00 .: 00* 0*0. . -*0 . 0 * - . . )> . . .00. 0 果下 = 同,接(归上的录之上线在地精但上这 $ ( ( )""" = , 0 0, . .0 , 0 ., . 0 0, ,. .( 下点散据以变随为强取变普 0 0 00 0 0 0 0 0 0 0 00 度/序数的碳

    注意事项

    本文(R软件一元线性回归分析(非常详细).doc)为本站会员(豆****)主动上传,淘文阁 - 分享文档赚钱的网站仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知淘文阁 - 分享文档赚钱的网站(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    关于淘文阁 - 版权申诉 - 用户使用规则 - 积分规则 - 联系我们

    本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

    工信部备案号:黑ICP备15003705号 © 2020-2023 www.taowenge.com 淘文阁 

    收起
    展开