Question

我有一个数据集，我使用来自41个人的10个估算器来估算家庭范围大小。我想测试这些估算器是否彼此显着不同，所以我使用nlme拟合线性混合效果模型如下：

mod3fML1b<-lme(size~band*nobsS+Sex,random=~1|snake,data=HR.compare,weights=vf6,method='REML')

大小是家庭范围大小（对数转换为满足常态假设），带是10个级别的因子代表我的10个估计量，nobsS是采样强度的度量（按抽样天数标准化的个别位置的数量），性是男/女，而蛇是我个体的随机效应。摘要输出如下：

> summary(mod3fML1b)$tTable
                              Value  Std.Error  DF     t-value      p-value
(Intercept)             4.099809927 0.24824427 351 16.51522508 9.632758e-46
bandBPIdiagonal        -0.344448847 0.05724340 351 -6.01726718 4.462854e-09
bandBPIfull            -0.303881612 0.05593369 351 -5.43289002 1.038615e-07
bandHREFdiagonal       -0.053639559 0.06749969 351 -0.79466377 4.273461e-01
bandHREFfull           -0.131436107 0.06471844 351 -2.03089130 4.301931e-02
bandLCVdiagonal        -0.186017321 0.11520159 351 -1.61471137 1.072717e-01
bandLCVfull            -0.186017321 0.11520176 351 -1.61470908 1.072722e-01
bandLCVsingle          -0.181687618 0.11940300 351 -1.52163366 1.290012e-01
bandSCVdiagonal        -0.163761675 0.05816320 351 -2.81555466 5.144255e-03
bandSCVfull            -0.120439828 0.05720672 351 -2.10534398 3.597148e-02
nobsS                   0.406335759 1.16315897  38  0.34933811 7.287640e-01
SexMale                 1.457327832 0.23373517  38  6.23495313 2.711438e-07
bandBPIdiagonal:nobsS  -0.415077222 0.24845863 351 -1.67060896 9.569032e-02
bandBPIfull:nobsS      -0.442855108 0.24277399 351 -1.82414559 6.898016e-02
bandHREFdiagonal:nobsS -0.143832274 0.29297492 351 -0.49093716 6.237777e-01
bandHREFfull:nobsS     -0.007949556 0.28090319 351 -0.02829999 9.774390e-01
bandLCVdiagonal:nobsS  -2.072417895 0.50001973 351 -4.14467228 4.270077e-05
bandLCVfull:nobsS      -2.072417895 0.50002044 351 -4.14466640 4.270182e-05
bandLCVsingle:nobsS    -2.171179830 0.51825545 351 -4.18940086 3.542319e-05
bandSCVdiagonal:nobsS  -0.745834058 0.25245092 351 -2.95437248 3.344486e-03
bandSCVfull:nobsS      -0.766847848 0.24829943 351 -3.08839960 2.172757e-03

p值表明几个级别与参考级别显着不同，并且模式与我基于原始数据的预期一致。但是，我希望以95％置信区间（非预测区间）显示平均对数（原点范围大小），以比较不同估计值之间的效果大小。为此，我使用predict.lme函数计算预测值，使用先前R线程中的代码计算SE： https://stat.ethz.ch/pipermail/r-sig-mixed-models/2010q1/003320.html 然而，我的95％CI广泛重叠，远远超过我根据我的lme模型的p值所预期的。我复制了我获得的值，显示了95％CI中的广泛重叠：

data<-data.frame(band=c('BPId','BPIf','REFs','REFd','REFf','LCVs','LCVd','LCVf','SCVd','SCVF'),
             log.mean.size=c(4.99,5.03,5.41,5.33,5.28,4.88,4.89,4.89,5.13,5.17),
             lci=c(4.62,4.65,5.03,4.95,4.90,4.90,4.50,4.50,4.75,4.79),
             uci=c(5.38,5.41,5.79,5.71,5.66,5.28,5.29,5.29,5.51,5.54))

有没有人建议为什么p值显示估算器级别之间的显着差异但CI重叠？我为没有提供任何数据而道歉，但如果有帮助，我很乐意通过电子邮件发送我的原始数据和完整代码的副本。

谢谢， JBauder

nlm p值和置信区间之间的差异

0 个答案: