在R语言中如何修改假设检验的显著性水平alpha值?
R中获取alpha=0.01对应t临界值的实现方法
R内置的t分布分位数函数qt()可以直接实现自定义alpha水平的t值查询,不需要依赖默认的0.05阈值,具体用法如下:
基础参数说明
qt()的核心入参有两个:
p:对应t分布的左侧累积概率,需要根据你使用的检验类型(单侧/双侧)调整df:t检验对应的自由度
不同检验场景的用法
双侧检验(最常用)
alpha=0.01对应的双侧检验,左右两侧各分摊0.005的拒绝域,因此累积概率取1 - 0.01/2 = 0.995。
示例:查询自由度为10时,alpha=0.01双侧检验的t临界值,代码如下:
qt(p = 0.995, df = 10)
运行输出结果约为3.169,和标准t界值表的对应数值完全一致。
单侧检验
alpha=0.01对应的单侧检验,拒绝域全部在分布一侧,因此累积概率取1 - 0.01 = 0.99。
示例:查询自由度为10时,alpha=0.01单侧检验的t临界值,代码如下:
qt(p = 0.99, df = 10)
t检验直接指定alpha=0.01水平
如果你是在执行t检验时需要输出对应alpha=0.01的置信区间、检验判定结果,直接在t.test()函数中指定conf.level = 0.99参数即可,示例如下:
# 独立样本t检验指定alpha=0.01水平 t.test(x = 实验组数据, y = 对照组数据, conf.level = 0.99)
内容的提问来源于stack exchange,提问作者Yellow Singh
相关产品推荐
相关产品推荐

