r - 使用 brms 包中的 ranef() 提取随机效应

Question

我正在使用 brms 包中的 ranef() 管理随机效应的结果。

bmodel<- brm(pop ~  RDB2000pop + Temperature2003 + Population2003 + 
                     (1+RDB2000pop+Temperature2003+Population2003 |species_id),
                   data      = dfpop_chenv,
                   family    = cumulative(link = "logit", threshold = "flexible"),
                   prior   = c(set_prior("normal(0,10)", class = "b")), 
                   warmup    = 1000,
                   iter      = 10000,
                   chains    = 4,
                   save_all_pars = TRUE) 

ranef(bmodel)[[1]]

我想提取“Population2003”的结果。

ranef(bmodel)$Population2003
NULL

如何提取模型中特定变量的结果？