如何计算非标准频带的置信度？

2018-06-20 13:02:35

所以我有一个数据集： 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可以做到这一点？

我不知道内置函数，但使用pt()函数计算并不太难，它是Student's 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统计量表”（减去平均值和除以sd）和（2）计算x<lower_bound和x<upper_bound的累积概率

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

值落入该范围的概率是这两个值的差异。

diff(probs)  ## 0.08805229

链接地址: http://www.djcxy.com/p/57763.html