arthur-albuquerque

library(marginaleffects)
library(rms)

load(url(
  "https://github.com/vincentarelbundock/modelarchive/raw/main/data-raw/gusto.rda"
))

gusto <- subset(gusto, tx %in% c("tPA", "SK"))
gusto$tx <- factor(gusto$tx, levels = c("tPA", "SK"))

arthur-albuquerque

pacman::p_load(rio, here, brms, dplyr, ggplot2, tidyr)
remotes::install_github("stan-dev/cmdstanr")

d = import(here("/Users/arthur/Downloads/il6_paper.xlsx"))

d = 
  d |> 
  mutate(phecode = as.factor(phecode),
         yi = log(OR),
         # Calculate mean effect size and SE in log scale

arthur-albuquerque

# data from https://cran.microsoft.com/snapshot/2022-07-18/web/packages/ordbetareg/vignettes/package_introduction.html

pacman::p_load(ordbetareg,
              brms)
data("pew")

model_data <- select(pew,therm,age="F_AGECAT_FINAL",
                     sex="F_SEX_FINAL",
                     income="F_INCOME_FINAL",
                     ideology="F_IDEO_FINAL",

arthur-albuquerque

This code was inspired by this post: https://discourse.datamethods.org/t/a-randomized-trial-of-epinephrine-in-out-of-hospital-cardiac-arrest

Install/Load packages

pacman::p_load(brms,
               tidybayes,
               ggplot2,
               dplyr,
 tidyr,

arthur-albuquerque

# Install/Load package

pacman::p_load(brms,
               tidybayes,
               ggplot2,
               dplyr,
               tidyr,
               patchwork)

# Data from

arthur-albuquerque

The objective of this code is to compare outputs from multiple R packages on the “average marginal effect” with a random-effect model fitted with brms.

pacman::p_load(dplyr,
               tidyr,
               tidybayes,
               brms,
               marginaleffects,
               brmsmargins,
               emmeans)

arthur-albuquerque

Make some data

library(MASS)
library(marginaleffects)
library(dplyr)
library(broom)

arthur-albuquerque

Make some data

library(MASS)
library(marginaleffects)
library(dplyr)
library(broom)

arthur-albuquerque

Make some data

library(MASS)
library(marginaleffects)
library(dplyr)
library(broom)

arthur-albuquerque

``` r
# Code based on 
# https://gist.github.com/mattansb/8dd0fb5a0be86e958ef6bdc4c49ac02c?permalink_comment_id=4241354#gistcomment-4241354
# https://gist.github.com/arthur-albuquerque/1a126043fef15b55792392e91f969f83

library(MASS)
library(marginaleffects)
library(dplyr)
#> 
#> Attaching package: 'dplyr'