MCplot: Plot of power function.

In akurz1/tclboot: Conditional Inference in Small Sample Scenarios

Plot of power function.

Description

This function plots the power function of selected test statistics

Usage

MCplot(
  object,
  nstats = c("W", "LR", "RS", "G", "abs", "sqs", "St"),
  alpha = 0.05,
  xlim = 5,
  legend_title = ggplot2::element_blank(),
  legend.direction = "vertical",
  tag_title = ggplot2::element_blank(),
  name_title = ggplot2::element_blank(),
  tcolor = FALSE
)

Arguments

object	An object of class "MCplot"
nstats	Name(s) of test statistics to be plotted. Note a maximum of 4 tests allowed. By default c("W", "LR", "RS", "G") will be plotted. Note "U1" = "sqs" or sum of squared elements of the score function and "U2" = "abs" or sum of absolute values.
alpha	Probability of error of first kind (in plot shown as dotted horizontal line).
xlim	A numeric value, specifying the left/lower limit and the right/upper limit of the x scale. Maximum value is xlim = 10.
legend_title	The text for Legend. By default it is set to "Test".
legend.direction	layout of items in legends ("horizontal" or "vertical")
tag_title	The text for tag label, which will be displayed at the top-left of the plot, by default of Item 1 only. By default it is element_blank().
name_title	The text for title of plot. By default it is element_blank().
tcolor	Logic value. If TRUE then ggplot plot will be in color otherwise black and white.

Value

An object of class "MCplot" containing:

plotlist	A list of length k containing ggplot objects of power plots of each item and selected test statistics.
call	The matched call.

References

Draxler, C., & Dahm, S. (2020). Conditional or Pseudo Exact Tests with an Application in the Context of Modeling Response Times. Psych, 2(4), 198-208. https://doi.org/10.3390/psych2040017

Examples

## Not run: 
# Rasch model with sample size N = 60, k = 8
library(tclboot)

# model parameters
items <- c(0, -2, -1, -0.5, 0, 0.5, 1, 2)
x <- c(rep(0,25), rep(1,35))

# generate initial data matrix
y <- MCsimRasch(N = 60, splitcr = x, items = items)$X

obj <- MCboot(X=y, splitcr = x, n_iter = 1000)

res <- MCplot(object = obj, nstats = c("W", "G"))


## End(Not run)

akurz1/tclboot documentation built on Oct. 23, 2022, 9:07 p.m.