r - 从 varComp 包中的 varCompFixEf 对象中提取 p 值（在 var comp 包中提取 p 值）

Question

我正在尝试从“varComp”包中指定的线性混合效应模型中提取固定效应的 p 值（使用 varComp()）函数。

使用 summary(model)$fixef

使我能够提取固定效果的输出表：

                 Estimate  Std. Error        Lower         Upper     t value Scale        Df     Pr(>|t|)
(Intercept) -10.713652206 0.117801777 -10.94843869 -10.478865719 -90.9464396     1  72.85480 9.330543e-77
snp1GG       -0.077577464 0.082086180  -0.24805847   0.092903541  -0.9450734     1  21.46851 3.551374e-01
SexM          0.198189186 0.071772595   0.05363726   0.342741109   2.7613490     1  45.06217 8.298709e-03
age.years    -0.007443209 0.004365628  -0.01603946   0.001153045  -1.7049574     1 261.52372 8.938987e-02
Average.BCS  -0.065434158 0.030612086  -0.12568251  -0.005185811  -2.1375269     1 291.94797 3.338572e-02
Seasonwet    -0.140436443 0.049207634  -0.23732376  -0.043549130  -2.8539564     1 265.21902 4.658585e-03

但是，我想提取 snp1GG 的 p 值，因为输出表是一个 varCompFixEf 对象，我不能使用诸如summary(model)$fixef[2,8]提取它的代码。

有谁知道如何从 varCompFixEf 对象中的一个固定效果中获取 p 值？

标签： rp-valuevarcomp