0

我写信是为了找出对象"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
4

1 回答 1

0

从文档到什么不清楚stanreg-objects

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

于 2018-03-12T22:07:16.633 回答