我正在此模型上的RStudio中运行“ step”功能:
inputData.entry = lmer(height ~ ENTRY_NO + REP + (1|SUB.BLOCK), data=inputData); # our model
这就是我使用“ step”运行的内容:
help.search("step",package="lmerTest");
st <- step(inputData.entry, reduce.fixed=FALSE);
print(st);
以下是输出:
Backward reduced random-effect table:
Eliminated npar logLik AIC LRT Df Pr(>Chisq)
<none> 142 -397.15 1078.3
(1 | SUB.BLOCK) 1 141 -397.47 1076.9 0.63157 1 0.4268
Backward reduced fixed-effect table:
Eliminated Df Sum of Sq RSS AIC F value Pr(>F)
ENTRY_NO 0 138 4238.1 6210.4 844.18 1.9775 5.749e-05 ***
REP 0 1 30.6 2002.9 816.03 1.9720 0.1627
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Model found:
height ~ ENTRY_NO + REP
我的问题是声明----找到模型:
当inputData.entry在“ lmer”中显示结果时,为什么结果列表中不会在模型中显示“ Sub.Block”?
我在做错什么吗?
谢谢您的建议!
答案 0 :(得分:0)
Step函数根据AIC选择模型。找到的模型告诉您根据这些条件选择了哪个模型。在这种情况下,该算法告诉我们Sub.Block不会因其引起的复杂情况而将有意义的信息带入分析,因此不在建议的模型中。