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_nb_ci, package = "jeksterslabRmedsimple")
data(results_mvn_sem_mc.mvn_ci, package = "jeksterslabRmedsimple")
results_mvn_complete_ci <- rbind(
  results_mvn_sem_mc.mvn_ci,
  results_mvn_nb_ci
)
results_mvn_complete_ci <- results_mvn_complete_ci[which(results_mvn_complete_ci$n > 20), ]
results_mvn_complete_power_ci <- results_mvn_complete_ci[which(results_mvn_complete_ci$n > 20 & results_mvn_complete_ci$n < 150), ]

Plots

alpha 0.05

width <- 11.7
height <- 8.3
cairo_pdf("complete-type1.pdf", width = width, height = height)
plot_type1(results_mvn_complete_ci)
dev.off()
cairo_pdf("complete-power.pdf", width = width, height = height)
plot_power(results_mvn_complete_power_ci)
dev.off()
cairo_pdf("complete-miss.pdf", width = width, height = height)
plot_miss(results_mvn_complete_ci)
dev.off()

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