## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(persuasio)

## -----------------------------------------------------------------------------
persuasio(
   est = "apr",
   y = "voteddem_all",
   t = "readsome",
   z = "post",
   level = 0.80,
   method = "normal",
   data = GKB
 )

## -----------------------------------------------------------------------------
persuasio(
   est = "lpr",
   y = "voteddem_all",
   t = "readsome",
   z = "post",
   level = 0.80,
   method = "normal",
   data = GKB
 )

## -----------------------------------------------------------------------------
persuasio(
  est = "yz",
  y = "voteddem_all",
  z = "post",
  x = "MZwave2",
  data = GKB,
  level = 0.80,
  model = "interaction",
  method = "bootstrap",
  nboot = 100
)

## -----------------------------------------------------------------------------
voteddem_all_0 <- mean(GKB$voteddem_all[GKB$post == 0], na.rm = TRUE)
voteddem_all_1 <- mean(GKB$voteddem_all[GKB$post == 1], na.rm = TRUE)

readsome_0 <- mean(GKB$readsome[GKB$post == 0], na.rm = TRUE)
readsome_1 <- mean(GKB$readsome[GKB$post == 1], na.rm = TRUE)

## -----------------------------------------------------------------------------
calc4persuasio(
  y1 = voteddem_all_1,
  y0 = voteddem_all_0,
  e1 = readsome_1,
  e0 = readsome_0
)

## -----------------------------------------------------------------------------
aprlb(y = "voteddem_all", z = "post", data = GKB)

## -----------------------------------------------------------------------------
aprlb(y = "voteddem_all", z = "post", x = "MZwave2", data = GKB)

## -----------------------------------------------------------------------------
gkb_groups <- split(GKB, GKB$MZwave2)
lapply(gkb_groups, function(sub_data) {
  aprlb(y = "voteddem_all", z = "post", data = sub_data)
})

## -----------------------------------------------------------------------------
aprub(y = "voteddem_all", t = "readsome", z = "post", data = GKB)

## -----------------------------------------------------------------------------
aprub(y = "voteddem_all", t = "readsome", z = "post", x = "MZwave2", data = GKB)

## -----------------------------------------------------------------------------
gkb_groups <- split(GKB, GKB$MZwave2)
lapply(gkb_groups, function(sub_data) {
  aprub(y = "voteddem_all", t = "readsome", z = "post", data = sub_data)
})

## -----------------------------------------------------------------------------
lpr4ytz(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  data = GKB
)

## -----------------------------------------------------------------------------
lpr4ytz(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  x = "MZwave2",
  data = GKB
)

## -----------------------------------------------------------------------------
lpr4ytz(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  x = "MZwave2",
  model = "interaction",
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4ytz(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  method = "normal",
  level = 0.80,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4ytz(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  method = "bootstrap",
  level = 0.80,
  nboot = 100,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4ytz(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  x = "MZwave2",
  model = "interaction",
  method = "bootstrap",
  level = 0.80,
  nboot = 100,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4ytz2lpr(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  method = "normal",
  level = 0.80,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4ytz2lpr(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  method = "bootstrap",
  level = 0.80,
  nboot = 100,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4ytz2lpr(
  y = "voteddem_all",
  t = "readsome",
  z = "post",
  x = "MZwave2",
  model = "interaction",
  method = "bootstrap",
  level = 0.80,
  nboot = 100,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4yz(
  y = "voteddem_all",
  z = "post",
  method = "normal",
  level = 0.80,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4yz(
  y = "voteddem_all",
  z = "post",
  method = "bootstrap",
  level = 0.80,
  nboot = 100,
  data = GKB
)

## -----------------------------------------------------------------------------
persuasio4yz(
  y = "voteddem_all",
  z = "post",
  x = "MZwave2",
  model = "interaction",
  method = "bootstrap",
  level = 0.80,
  nboot = 100,
  data = GKB
)