最新最新R语⾔⾃主练习题语⾔⾃主练习题
1,某⼯⼚⽣产⼀批滚珠,其直径服从正态分布N(µ,σ2),现从某天的产品中随机抽出六件,测得直径为:
15.1,14.8,15.2,14.9,14.6,15.1,。若σ2 =0.06,求µ的置信区间。(置信度为0.95)
解:
x<-c(15.1,14.8,15.2,14.9,14.6,15.1)
sigema<-sqrt(0.06)
alpha<-0.05
xbar<-mean(x)
n<-length(x)
t1<-xbar-qnorm(1-alpha/2)*sigema/sqrt(n)
t2<-xbar+qnorm(1-alpha/2)*sigema/sqrt(n)
list(t1,t2),
2,某某⾃动包装机包装洗⾐粉,其重量ζ~N(µ,σ2),其中µ,σ未知。今随机抽取⼗⼆袋测得其重量,经计算得样本均值
为xbar=1000.25,修正样本标准差s*=2.6329,试求总体标准差σ的置信⽔平为0.95的置信区间。
解:
alpha<-0.05
Xbar<-1000.25
Sdx<-2.6329
T1<-sqrt(11)*Sdx/sqrt(qchisq(1-alpha/2,11))
T2<-sqrt(11)*Sdx/sqrt(qchisq(alpha/2,11))
list(T1,T2)
使⽤t.text函数进⾏⽅差未知的均值假设检验
t检验t.test():调⽤格式:(数统P138,例6-3)
x <- c(11.6,11.5,11.3,11.2,11.4,11.7,11.5,11.6,11.4,11.3)
α<- 0.05
solution <- t.test(x,mu=11.4,alternative="two.sided",conf.level = 1-α)
#x是⼀个服从正态分布的总体,mu是均值µ
#alternative是指备择假设,“two.sided”(缺省)指双侧(H1:µ≠µ0),less表⽰单边检验(H1:µ<µ1),greater表⽰单边检验(H1:µ>µ1)
#conf.level指置信度即1-α
solution
if(solution$p.value>α){
print("接受H0")
}else{
print("拒绝H0,接受H1")
}
#如果p-value>α,则可以认为接受H0,否则拒绝H0,接受H1