多元统计分析及R语言建模（第五版）——第5章 广义与一般线性模型课后习题

第5章 广义与一般线性模型

文章会用到的数据请在这个网址下下载多元统计分析及R语言建模（第五版）数据
1）现有甲、乙、丙三个工厂生产同一种零件，为了了解不同工厂的零件强度有无明显的差异，现分别从每一个工厂随机抽取部分零件测定其强度，…试问三个工厂的零件的平均强度是否相同？

x1 = c(103,101,98,110)
x2 = c(113,107,108,116,115,109)
x3 = c(82,92,84,86,88) 

mean(x1)

mean(x2)

mean(x3)

library(openxlsx)
d5.1 <- read.xlsx('mvexer5.xlsx',sheet = 'E5.1')
d5.1

anova(lm(零件强度 ~ factor(工厂),data = d5.1))

P<0.05,说明不同工厂的零件强度有明显差异，根据上述结果可以看出三个工厂的零件的平均强度不相同
2）生产某种化工产品时，要比较四种不同的配方对生产率的影响。考虑到生产率随生产日期的不同而变动较大，所以把实验日期也选为因子。实验分四天进行。配方因子和