在这个关于期刊报告的问题之后,我想知道是否有适合使用的贝叶斯模型的总体效应大小指数stan_lmer
?
在频率论框架中,例如有一个伪 R2(由这个包计算),它返回边际(仅由固定因子解释的方差比例)和条件(由两者解释的方差比例)固定和随机因素)R2。是否有任何等价物可以帮助我们量化和限定效果/预测能力大小?
谢谢你。
对于stan_lm
或stan_glm
,family = gaussian(link = "identity")
您可以计算 R 平方的后验分布
R2 <- rowSums(posterior_linpred(post)^2) / rowSums(posterior_predict(post)^2)
orpost
返回的对象在哪里。For or (equivalently) with you 可以做类似的事情,但必须对如何处理特定于组的术语(即)采取立场。和的参数为您提供了几个选项,但默认情况下,特定于组的术语是有条件的。或者,您可以将(其中一些)设置为零或通过新的组级别传递 a 以集成(一些)组特定的术语。stan_lm
stan_glm
stan_lmer
stan_glmer
family = gaussian(link = "identity")
(1 | group)
re.form
posterior_linpred
posterior_predict
data.frame