使用 lmer 函数获取混合模型运行的 p 值

时间:2020-12-26 01:35:06

标签: r lme4

我使用 lmer 运行了一些混合模型,但它们没有给出 p 值。我想知道是否有办法获得这些模型的 p 值。有人建议使用 afex 包。我已经调查过了,感到困惑和不知所措。例如,在 https://rdrr.io/rforge/afex/man/mixed.html 处,它给出了看起来非常复杂和复杂的代码;它是如此势不可挡,让我怀疑这是否真的是我需要做的!下面是我运行的混合模型的示例;我想获得固定效应的 p 值和固定效应的相关性。任何帮助将不胜感激!

Linear mixed model fit by REML ['lmerMod']
Formula: score ~ group + condition + (1 | subject) + (1 | token_set) + (1 | list)
Data: EN_JT_1

REML criterion at convergence: 744.9

Scaled residuals:
Min 1Q Median 3Q Max
-3.5860 -0.0364 0.2183 0.5424 1.6575

Random effects:
Groups Name Variance Std.Dev.
subject (Intercept) 0.006401 0.08000
token_set (Intercept) 0.001667 0.04083
list (Intercept) 0.000000 0.00000
Residual 0.084352 0.29043
Number of obs: 1704, groups: subject, 71; token_set, 24; list, 2

Fixed effects:
Estimate Std. Error t value
(Intercept) 0.99796 0.02425 41.156
groupHS -0.08453 0.02741 -3.084
groupSB -0.03103 0.03034 -1.023
conditionEN-GJT-D-ENG -0.10329 0.01990 -5.190
conditionEN-GJT-D-NNS -0.01288 0.02617 -0.492
conditionEN-GJT-D-NTR -0.19250 0.02596 -7.415

Correlation of Fixed Effects:
(Intr) gropHS gropSB cEN-GJT-D-E cEN-GJT-D-NN
groupHS -0.452
groupSB -0.409 0.361
cEN-GJT-D-E -0.410 0.000 0.000
cEN-GJT-D-NN -0.531 0.000 0.000 0.380
cEN-GJT-D-NT -0.535 0.000 0.000 0.383 0.700
optimizer (nloptwrap) convergence code: 0 (OK)
boundary (singular) fit: see ?isSingular

0 个答案:

没有答案
相关问题