所以我有一个数据集：
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之间的可能性。哪个包::功能可以做到这一点？
答案 0 :(得分：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统计”标度（减去平均值并除以sd）和（2）计算x<lower_bound和x<upper_bound的累积概率
probs <- pt((c(850,900)-m)/s,df=length(dd)-1)
值落在该范围内的概率是这两个值的差异。
diff(probs) ## 0.08805229