利用R语言进行线性拟合.docxVIP

下载本文档

60
0
约1.57千字
约 3页
2017-11-19 发布于浙江
举报
版权申诉

利用R语言进行线性拟合.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

利用R语言进行线性拟合

利用R语言进行线性/非线性回归拟合实例生成一组数据vectorfloat xxvec;vectorfloat yyvec;ofstream fout(data2.txt);for (inti =1;i200;i++){float x = i*0.8;float randdnum = rand()%10 * 10;float randomflag = (rand()%10)%2==0?(1):(-1);float y = 3 * x*x + 2*x + 5 + randomflag*randdnum;foutx yendl;xxvec.push_back(x);yyvec.push_back(y);}fout.close();将生成的数据存为txt文件，命名为“data1”2、线性拟合#-------------------------------------------------------------#载入数据 fire - read.table(D:/data.txt, header = TRUE)#-------------------------------------------------------------#回归分析 plot(fire$y ~ fire$x) fire.reg - lm(fire$y ~ fire$x, data = fire) #回归拟合 data1.regCall:lm(formula = data1$y ~ data1$x, data = data1)Coefficients:(Intercept) data1$x 6.202 3.003 summary(data1.reg) #回归分析表Call:lm(formula = data1$y ~ data1$x, data = data1)Residuals: Min 1Q Median 3Q Max -93.345 -42.929 -1.948 46.717 88.793 Coefficients: Estimate Std. Error t value Pr(|t|) (Intercept) 6.202084 3.352055 1.85 0.0646 . data1$x 3.002826 0.007259 413.66 2e-16 ***标红数字即为线性回归系数，由于生成数据时加了一个随机数，所以拟合出来的直线为：y=3.002826 x+6.202084---Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1Residual standard error: 52.93 on 997 degrees of freedomMultiple R-squared: 0.9942,Adjusted R-squared: 0.9942 F-statistic: 1.711e+05 on 1 and 997 DF, p-value: 2.2e-16 anova(data1.reg) #方差分析表Analysis of Variance TableResponse: data1$yDf Sum Sq Mean Sq F value Pr(F) data1$x 1 479462873 479462873 171112 2.2e-16 ***Residuals 997 2793641 2802 Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 abline(data1.reg, col = 2, lty = 2) #拟合直线