第二次 R 语言作业——P144-P145练习题1，2，3

1） 现有甲乙丙三个工厂生产同一种零件，为了了解不同工厂的零件的度有无明显的差异。现分别从每一个工厂随机抽取部分零件，测定其强度。数据如图表 5~9 所示：

• 试问三个工厂的零件的平均强度是否相同？
  
  由上图可得知

P<0.05,说明不同工厂的零件强度有明显差异，可以看出三个工厂的零件的平均强度不相同。

生产某种化工产品时，要比较四种不同配方对生产率的影响，考虑到生产率随生产日期不同而变动较大，所以把实验日期也选为因子。实验分四天进行，配方因子和日期因子分别用 ab 表示数据如图。5-10 所示。

综上分析得：

P(A)<0.05，说明不同配方因子 A 对生产率有影响
P(B)<0.05，说明不同生产日期 B 对生产率有影响

**

3） 某银行从历史贷款客户中随机抽取 16 个样本，根据设计的指标体系分别计算他们的商业信用支持度 X1 和市场竞争地位等级 X2。类别变量记中 1 代表贷款成功，2代表贷款失败数据，如表所示。

（1） 为了给正确贷款提供决策支持，请建立 Logistic 模型进行分析。已知表格内容：

第一步读入数据：
输入命令：
> E5.3 <- read_excel("G:/学习 app/作业合集/R 语言 /test/test01/练习题3.xlsx","E5.3")#查找并读取数据对应的表格

第二步 Logistic 回归模型：

> 输入命令分别为： logit.glm = glm(G ~ x1 + x2,data = E5.3)#建立 Logistic 回归模型结果
> summary(logit.glm)#Logistic 回归模型结果

分析：

由图可知 logistic 拟合的模型 x1（商业信用支持度）的 P=0.5712>0.05，x2（市场竞争地位等级）的 P=0.0013<0.05，说明商业信用支持度对贷款没有显著影响，市场竞争地位等级对贷款有影响。

（2） 根据建立的模型判定是否给某客户 X1=131， X2=-2 提供贷款。

第一步逐步筛选法进行变量选择：
输入命令：

> logit.step <- step(logit.glm,direction = 'both')#逐步筛选法变量选择

逐步筛选法进行变量选择结果：
输入命令：

> summary(logit.step)#结果

第二步进行标椎化：
输入命令：

>pre1 <- predict(logit.glm,data.frame(x1 = 131,x2 = -2))#Logistic 回归模型预测
>pre1#输出查看结果
>pre2 <- predict(logit.step,data.frame(x1 = 131,x2 = -2))#逐步筛选法预测
>pre2#输出查看结果

由上述结果可知，两种方法的预测值都为 2，则表示贷款失败，所以不会提供贷款。

5） 考虑一个化学反应过程,有两个因素,因素 A 为反应物的浓度，它有两个水平,A0(15%),A1(25%).因素 B 为催化剂是否使用.它有两个水平 B0(不用),B1(用),每种组合做 3 次试验结果如表 5-13：

(1) 试写出其一般线性模型的矩阵表示；

处理因素 A 有 G 个水平，单位组 B 有 n 个水平，分别产生 A 的 G 个哑变量和单位白能量后，将实验结果 yij 表示成：

其中 u 为平均数，ai 为处理因素 A 的第 i 个水平的效应；bj 为第 J 个单位组的效应，cij为误差项也被称之为残差。最后可得矩阵表达式为：

其一般线性模型的矩阵表示为：
其中 X 时设计阵，元素为 0 或 1，是误差向量，Y 为观察结果向量（处理因素为反应物浓度 A，单位组是催化剂 B。
(2) 分析因此 A,B 和交互作用 A*B 对化学反应的影响.
导入数据：

 >library(readxl) #导入包
> X5_13 <- read_excel("C:/Users/13623/Desktop/5-13.xlsx")#导入数据

有图中数据可知信息如下所示：

P（A）<0.05，说明反应物的浓度对该化学反应有显著影响
P(B)<0.05，说明催化剂对该化学反应有显著影响
P(A：B)>0.05，说明反应物的浓度和催化剂交互作用对该化学反应没有影响

7）磁鼓电机是彩色录像机磁鼓组件的关键部件之一，按质量要求其输出力矩应大于 210g.cm。某生产厂过去这项指标的合格率较低，从而希望通过试验找出好的条件，以 提高磁鼓电机的输出力矩，表 5-15 是磁鼓电机试验数据。

对上述正交设计进行方差分析。
导入数据：

> library(readxl)#导入包
> X5_15 <- read_excel("C:/Users/13623/Desktop/X5_15.xlsx")#导入数据
> > X5_15#查看数据

> anova(lm(实验结果输出~充磁量+定位角度+钉子线圈+钉子线圈*充磁量*定位角度+充磁量*定位角度+充磁量*钉子线圈+定位角度*钉子线圈,data = X5_15))

由上图可知数据：

P（充磁量）>0.05，说明充磁量对输出力矩没有显著影响。 2、P（定位角度）>0.05，说明定位角度对输出力矩没有显著影响。
P（充磁量：钉子线圈）>0.05，说明充磁量和钉子线圈对输出力矩没有显著影 响。
P（定位角度：钉子线圈）<0.05，说明定位角度和钉子线圈对输出力矩有显著影响。
P（充磁量：钉子线圈:定位角度）>0.05，说明充磁量、定位角度和钉子线圈对输 出力矩没有显著影响。