关于Rmisc包summarySE()返回置信区间(ci)数值异常的疑问

你混淆了正态分布（z分布）和t分布的临界值使用场景，这是导致结果差异的核心原因：

• 你手动计算时用的1.96是正态分布下95%置信区间的临界值，但summarySE()默认采用t分布计算置信区间——当样本量较小（每组n=10，自由度df=10-1=9）时，t分布更适合小样本的置信区间估计。

验证计算过程 #

对于constructionA组：

1. 先计算自由度为9的95%置信区间t临界值：

qt(0.975, df=9) # 输出约2.262157

2. 用该t值乘以标准误（se=0.1261062）：

2.262157 * 0.1261062 ≈ 0.2852721

这个结果和summarySE()输出的ci值完全一致。

补充说明 #

summarySE()的底层逻辑是基于t检验的置信区间计算，而非正态分布的z值。如果你的场景确实需要使用z分布计算，可以手动计算ci值，或者自定义函数实现该逻辑。

内容的提问来源于stack exchange，提问作者Lawrence Lam