统计学基于r语言第六章和第七章答案

#6.1(1)
load('C:/exercise/ch6/exercise6_1.RData')
par(mfrow=c(1,2),cex=0.8,mai=c(0.7,0.7,0.2,0.1))
qqnorm(exercise6_1$零件误差,xlab="期望正态值",ylab="观测值",datax=TRUE,main="正态Q-Q图")
qqline(exercise6_1$零件误差,datax=T,col='red')
##提出假设：H0:误差服从正态分布 H1：误差不服从正态分布
#shapiro.test(exercise6_1$零件误差)
##结论：由于p<0.01 拒绝原假设 有证据证明绝对误差不服从正态分布

#(2)
#提出假设: H0：有显著降低  H1:无显著降低
install.packages('BSDA')
library(BSDA)
z.test(exercise6_1$零件误差,mu=1.35,sigma.x=sd(exercise6_1$零件误差),alternative="less",conf.level=0.95)
#结论：在该检验中，z=-2.6061，P=0.004579，由于P<0.01，拒绝H0，新机床加工的零件尺寸的平均误差与就机床相比无显著降低

#6.2（1）
load('C:/exercise/ch6/exercise6_2.RData')
#S-W:
#提出假设：H0:重量服从正态分布 H1：重量不服从正态分布
shapiro.test(exercise6_2$重量)
#结论：p=0.7684>0.05 有证据表明原假设成立重量服从正态分布

#K-S:
#提出假设：H0:重量服从正态分布 H1：重量不服从正态分布
ks.test(exercise6_2$重量,"pnorm",mean(exercise6_2$重量),sd(exercise6_2$重量))
#结论：p=0.9539>0.05 有证据表明原假设成立重量服从正态分布

#(2)
#提出假设：H0:金属板符合要求 H1：金属板不符合要求
z.test