R语言自主练习题.docx

1,某工厂生产一批滚珠,其直径服从正态分布N(μ,σ2),现从某天的产品中随机抽出六件,测得直径为:,,,,,,。若σ2=,求μ的置信区间。()解:x<-c(,,,,,)sigema<-sqrt()alpha<-<-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=,修正样本标准差s*=,。解:alpha<-<-<-<-sqrt(11)*Sdx/sqrt(qchisq(1-alpha/2,11))T2<-sqrt(11)*Sdx/sqrt(qchisq(alpha/2,11))list(T1,T2)():调用格式:(数统P138,例6-3)x<-c(,,,,,,,,,)α<-<-(x,mu=,alternative="",=1-α)#x是一个服从正态分布的总体,mu是均值μ#alternative是指备择假设,“”(缺省)指双侧(H1:μ≠μ0),less表示单边检验(H1:μ<μ1),greater表示单边检验(H1:μ>μ1)#-αsolutionif(solution$>α){print("接受H0")}else{print("拒绝H0,接受H1")}#如果p-value>α,则可以认为接受H0,否则拒绝H0,接受H1例题:水泥厂用自动包装机包装水泥,每袋额定重量是50kg,某日开工后随机抽查了9袋,称得重量如下:,,,,,,,,,问包装机工作是否正常(α=)?x<-c(,,,,,,,,)α<-<-(x,mu=50,alternative="",=1-α)solutionif(solution$>α){print("接受H0")}else{print("拒绝H0,接受H1")}p-value>α,接受H0认为包装机工作正常。例题:,有一个实验检测了该公司所制造的6套电池,得如下的寿命小时数:19,18,22,20,16,25这些结果是否表明这种类型的电池与该公司宣称的寿命不同?(α=)x<-c(19,18,22,20,16,25)α<-<-(x,mu=,alternative="",=1-α)solutionif(solution$>α){print("接受H0")}else{print("拒绝H0,接受H1")}p-value>α,接受H0认为这种类型的电池与该公司宣称的寿命相同例题:据长期的经验和资料分析,某砖瓦厂所生产的砖的“抗断强度”X服从正态分布,方差σ^2=,今从该厂所生产的一批砖中,随便抽取6块,测得抗断度如下:,,31,64,,,:?("BSDA")library("BSDA")x<-c(,,,,,)α<-<-(x,alternative="",==1-α)solutionif(solution$>α){print("接受H0")}else{print("拒绝H0,接受H1")}例2:片剂车间生产的一种药片,片重服从N(μ,),今抽取8片,测得每片