In jeksterslabds/jeksterslabRmedsimple: Jekster's Lab Simple Mediation Model Utilities

knitr::opts_chunk$set(
  error = TRUE,
  collapse = TRUE,
  comment = "#>",
  out.width = "100%"
)

library(jeksterslabRmedsimple)
data(results_mvn_mnar_mc.mvn_ci, package = "jeksterslabRmedsimple")
data(results_mvn_mnar_nb.fiml_ci, package = "jeksterslabRmedsimple")
results_mnar_ci <- rbind(
  results_mvn_mnar_mc.mvn_ci,
  results_mvn_mnar_nb.fiml_ci
)
results_mnar_10_ci <- results_mnar_ci[which(results_mnar_ci$missing == "MNAR.10"), ]
results_mnar_20_ci <- results_mnar_ci[which(results_mnar_ci$missing == "MNAR.20"), ]
results_mnar_30_ci <- results_mnar_ci[which(results_mnar_ci$missing == "MNAR.30"), ]

Plots

alpha 0.05

mnar_10

plot_type1(results_mnar_10_ci)
plot_power(results_mnar_10_ci)
plot_miss(results_mnar_10_ci)

mnar_20

plot_type1(results_mnar_20_ci)
plot_power(results_mnar_20_ci)
plot_miss(results_mnar_20_ci)

mnar_30

plot_type1(results_mnar_30_ci)
plot_power(results_mnar_30_ci)
plot_miss(results_mnar_30_ci)

jeksterslabds/jeksterslabRmedsimple documentation built on Oct. 16, 2020, 11:30 a.m.