什么是“linear.predictors”可以从“rstanarm”包中的stan_glm()对象中提取出来的?

时间:2018-03-12 21:31:35

标签: r bayesian rstan hierarchical-bayesian rstanarm

我写信是为了找出返回的"linear.predictors" stan_glm()对象。

显然,"linear.predictors"与用户提供的预测变量不同(文档没有帮助)。

在任何情况下,有没有办法从stan_glm()对象获取预测值?

以下是单个预测变量(即mom_iq)示例:

library(rstanarm) 
data(kidiq)
d <- kidiq  

fit <-stan_glm(kid_score ~ mom_iq,
           data = d,   
           prior = normal(0, 2.5),  
           prior_intercept = normal(0, 10),  
           prior_aux = cauchy(0, 100)) 


 max(fit$linear.predictors) # 110.5605   # As can be seen, these are not the same

 max(d$mom_iq)              # 138.8931

1 个答案:

答案 0 :(得分:0)

stanreg-objects的文档中不清楚哪些内容?

  

...
  fitted.values拟合平均值。对于GLM,线性预测变量通过反向链接函数进行变换   linear.predictors线性拟合链接比例。对于线性模型,这与拟合值相同   ...