首页 文章

如何计算非标准波段的置信度?

提问于
浏览
0

所以我有一个数据集: dt<-c(1106,1402, 827, 781,876, 1134,1014, 964, 848, 814, 772, 912, 923, 996, 569, 774, 1389, 900) 让我们假设一条正常曲线并运行 t.test(dt) 我得到了

一个样本t检验数据:dt t = 19.057,df = 17,p值= 6.579e-13备选假设:真实均值不等于0 95%置信区间:839.9344 1049.0656样本估计值:平均值x 944.5

这是95%置信区间的所有标准票价 .

但是我想要找到的是850到900之间特定范围内的信心 . 因为我想知道下一个数据点落在850到900之间的可能性 . 哪个package :: function可以做到这一点?

1 回答

  • 4

    我不太难以使用 pt() 来计算,这是学生t的累积分布函数:

    dd <- c(1106,1402, 827, 781,876, 1134,1014, 
            964, 848, 814, 772, 912, 923, 996, 569, 774, 1389, 900)
    m <- mean(dd)
    s <- sd(dd)
    

    现在我们(1)将所需范围转换为"t statistic" scale(减去均值并除以sd)和(2)计算 x<lower_boundx<upper_bound 的累积概率

    probs <- pt((c(850,900)-m)/s,df=length(dd)-1)
    

    值落在该范围内的概率是这两个值的差值 .

    diff(probs)  ## 0.08805229
    

相关问题