单个、双个正态分布区间估计--R code

最新推荐文章于 2024-07-11 09:07:14 发布

Vac_cat

最新推荐文章于 2024-07-11 09:07:14 发布

阅读量4.5k

点赞数 2

分类专栏： R 文章标签：数据分析 r语言

本文链接：https://blog.csdn.net/Vac_cat/article/details/78866494

版权

这篇博客分享了如何使用R语言进行单个和双个正态分布的区间估计，内容源于《数据分析与统计》课程的附加题，供读者参考并欢迎指出错误。

摘要由CSDN通过智能技术生成

代码是选修课《数据分析与统计》上附加题，可直接应用。欢迎纠错。

#*****************
#本程序利用R分布别做出单个、双个正态分布总体
#区间估计的均值和双侧、单侧置信区间
#包含4个函数
#code by Vac 2017.12.8
#*****************

#单个正态分布
#u的区间估计，sigma已知
onenorm_u_sigma<-function(X,sigma,alpha=0.99){
    n<-length(X);
    tmpx<-sigma*qnorm(1-alpha/2)/sqrt(n);
    a<-mean(X)-tmpx;b<-mean(X)+tmpx;#双侧
    cat("其",alpha*100,"%双侧置信区间是：[",a,",",b,"]","\n")

    bu<-mean(X)+sigma*qnorm(1-alpha)/sqrt(n)#单侧上限
    ad<-mean(X)-sigma*qnorm(1-alpha)/sqrt(n)#单侧下限
    cat("其",alpha*100,"%单侧置信区间上限是：",bu,"\n下限是",ad,"\n")
}

#u的区间估计，sigma未知
onenorm_u_sd<-function(X,alpha=0.99){
    n<-length(X);
    tmpx<-sd(X)*qt(1-alpha/2,n-1)/sqrt(n);
    a<-mean(X)-tmpx;b<-mean(X)+tmpx;#双侧
    cat("其",alpha*100,"%双侧置信区间是：[",a,",",b,"]","\n")

    bu<-mean(X)+sd(X)*qt(1-alpha,n-1)/