建模与R软件课程作业_第1页
建模与R软件课程作业_第2页
建模与R软件课程作业_第3页
建模与R软件课程作业_第4页
建模与R软件课程作业_第5页
已阅读5页,还剩5页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、餐馆售货机数量咖啡销售量(/杯)餐馆售货机数量咖啡销售量(/杯)10508.183697.520498.494755.331568.2104758.941577.3115787.652651.7125792.162657.0136841.473713.4146831.8 建模与R软件课程作业 学号:姓名:江俊逸成绩:1.一位饮食公司的分析人员想调查自助餐馆中的自动咖啡售货机数量与咖啡销售量的关系,他选择了14家餐馆来进行实验。这14家餐馆在营业额、顾客类型和地理位置方面都是相近的。放在试验餐馆的自动售货机数量从0(这里咖啡由服务员端来)到6不等,并且是随机分配到每个餐馆的。表1所示的是关于试验

2、结果的数据。(1)作线性回归模型;(2)作多项式回归模型表1自动咖啡售货机数量与咖啡销售量数据(3)画出数据的散点图和拟合曲线。提示:参考统计建模与R软件第6章的相关内容。解:#输入数据,作线性回归模型X-data.frame(x=c(0,0,1,1,2,2,3,3,4,4,5,5,6,6),y=c(508.1,498.4,568.2,577.3,651.7,657.0,713.4,697.5,755.3,758.9,787.6,792.1,841.4,831.8)lm.sol|t|) (Intercept) 523.800 8.474 61.81 2e-16 *x 54.893 2.350

3、23.36 2.26e-11 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1Residual standard error: 17.59 on 12 degrees of freedomMultiple R-squared: 0.9785, Adjusted R-squared: 0.9767 F-statistic: 545.5 on 1 and 12 DF, p-value: 2.265e-11回归系数与回归方的检验都是显著的,故回归方程为:#画散点图及拟合曲线plot(X$x,X$y)abline(lm.sol)结果如下: #作多项式拟合

4、lm.sol_1|t|) (Intercept) 502.5560 4.8500 103.619 2e-16 *x 80.3857 3.7861 21.232 2.81e-10 *I(x2) -4.2488 0.6063 -7.008 2.25e-05 *-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1Residual standard error: 7.858 on 11 degrees of freedomMultiple R-squared: 0.9961, Adjusted R-squared: 0.9953 F-statistic: 13

5、91 on 2 and 11 DF, p-value: 5.948e-14回归系数与回归方的检验都是显著的,故回归方程为:#画散点图及拟合曲线xfit-seq(0,6,len=200);yfit-predict(lm.sol_1,data.frame(x=xfit)plot(X$x,X$y)lines(xfit,yfit)结果如下2.为考查学生的学习情况。学校随机地抽取12名学生的5门课程期末考试成绩,其中数据见表2,试用因子分析的方法对这种数据进行分析。(1)找出5门课程的公共因子,并进行合理的解释;(2)用回归方法或Bartlett方法计算样本的因子得分,画出因子得分的第1第序号政治(X1

6、)语文(X2)外语(X3)数学(X4)物理(X5)199949310010029988969997310098819610049388889996510091729678690788275977757388978989384836888987736076841095829062391176724367781285755034372公共因子的散点图,通过散点图来分析这12名学生的学习情况。 提示:参考统计建模与R软件第9章的相关内容。解:#输入数据X-data.frame(x1=c(99,99,100,93,100,90,75,93,87,95,76,85),x2=c(94,88,98,88,9

7、1,78,73,84,73,82,72,75),x3=c(93,96,81,88,72,82,88,83,60,90,43,50),x4=c(100,99,96,99,96,75,97,68,76,62,67,34),x5=c(100,97,100,96,78,97,89,88,84,39,78,37)fa-factanal(X,factors=2); fa显示结果Call:factanal(x = X, factors = 2)Uniquenesses: x1 x2 x3 x4 x5 0.005 0.141 0.494 0.005 0.346 Loadings: Factor1 Factor

8、2x1 0.992 0.104 x2 0.854 0.360 x3 0.497 0.509 x4 0.284 0.956 x5 0.132 0.798 Factor1 Factor2SS loadings 2.059 1.950Proportion Var 0.412 0.390Cumulative Var 0.412 0.802Test of the hypothesis that 2 factors are sufficient.The chi square statistic is 0.28 on 1 degree of freedom.The p-value is 0.598 #相关矩阵r=cor(X)fa由以上两种方法可得相关系数矩阵和用原矩阵所得到的结果是一致的,在计算结果中,因子f1前两个变量(x1,x2)的载荷因子接近于1,这类归为“文科类”,因此称之为文科因子;因子f2后两个变量(x4,x5)的载荷因子接近于1,这类归为“理科类”,因此称之为理科因子。fa-factanal(X,factors=2,sc

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论