0

RでbayestestRandを使用してベイジアン ロジスティック回帰を実行したいと思います。出力はlog(odds ratio)にあると思います。中心性、不確実性、存在、有意性の指標など、すべてを代わりにオッズ比に変換する方法を知っていますか? パッケージの関数には引数があり、すべてをORで返すことを知っています。rstanarmtbl_summarygtsummaryexponentiate = TRUE

コード：

library(rstanarm)
library(bayestestR)

data <- iris %>%
  filter(Species != "setosa") %>%
  droplevels()

model <- stan_glm(Species ~ Sepal.Width, data = data, family = "binomial", refresh = 0)

describe_posterior(model)
# Summary of Posterior Distribution 
# 
# Parameter   | Median |          95% CI |     pd |          ROPE | % in ROPE |  Rhat |     ESS
# ---------------------------------------------------------------------------------------------
# (Intercept) |  -6.16 | [-10.46, -2.30] | 99.95% | [-0.18, 0.18] |        0% | 1.001 | 2842.00
# Sepal.Width |   2.15 | [  0.83,  3.64] | 99.92% | [-0.18, 0.18] |        0% | 1.001 | 2799.00
4

1 に答える 1