Nothing
## ----results='hide', message=FALSE, warning=FALSE-----------------------------
library(bayesian)
## ----results='hide', message=FALSE, warning=FALSE-----------------------------
library(recipes)
library(workflows)
## -----------------------------------------------------------------------------
epi_recipe <- epilepsy |>
recipe() |>
update_role(count, new_role = "outcome") |>
update_role(Trt, Age, Base, patient, new_role = "predictor") |>
add_role(patient, new_role = "group") |>
step_normalize(Age, Base)
## -----------------------------------------------------------------------------
print(epi_recipe)
## -----------------------------------------------------------------------------
epi_model <- bayesian(
family = poisson()
) |>
set_engine("brms") |>
set_mode("regression")
## -----------------------------------------------------------------------------
print(epi_model)
## -----------------------------------------------------------------------------
epi_model <- epi_model |>
update(family = poisson())
## -----------------------------------------------------------------------------
epi_workflow <- workflow() |>
add_recipe(epi_recipe) |>
add_model(
spec = epi_model,
formula = count ~ Trt + Base + Age + (1 | patient)
)
## -----------------------------------------------------------------------------
print(epi_workflow)
## ----results='hide', echo = FALSE---------------------------------------------
run_on_linux <- grepl("linux", R.Version()$os, ignore.case = TRUE)
## ----results='hide', eval = run_on_linux--------------------------------------
epi_workflow_fit <- epi_workflow |>
fit(data = epilepsy)
## ----eval = run_on_linux------------------------------------------------------
print(epi_workflow_fit)
## ----eval = run_on_linux------------------------------------------------------
epi_fit <- epi_workflow_fit |>
extract_fit_parsnip()
## ----eval = run_on_linux------------------------------------------------------
epi_brmsfit <- epi_workflow_fit |>
extract_fit_engine()
## ----eval = run_on_linux------------------------------------------------------
class(epi_brmsfit)
## -----------------------------------------------------------------------------
newdata <- epilepsy[1:5, ]
## ----eval = run_on_linux------------------------------------------------------
epi_workflow_fit |>
predict(
new_data = newdata,
type = "conf_int",
level = 0.95
)
## ----eval = run_on_linux------------------------------------------------------
epi_workflow_fit |>
predict(
new_data = newdata,
type = "conf_int",
level = 0.95,
std_error = TRUE
)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.