建模与R软件课程作业

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