p_mcnemar.test: p-value from McNemar test simulation

In Spower: Power Analyses using Monte Carlo Simulations

p-value from McNemar test simulation

Description

Generates two-dimensional sample data for McNemar test and return a p-value. Uses mcnemar.test.

Usage

p_mcnemar.test(
  n,
  prop,
  OR = NULL,
  prop.disc = NULL,
  two.tailed = TRUE,
  correct = TRUE,
  gen_fun = gen_mcnemar.test,
  return_analysis = FALSE,
  ...
)

gen_mcnemar.test(n, prop, ...)

Arguments

n	total sample size
prop	two-dimensional matrix of proportions/probabilities
OR	instead of supplying the prop table, the odds ratio can be specified instead \pi_{12}/\pi_{21}. Also requires proportion of discordant pairings to be specified
prop.disc	proportion of discordant pairings, \pi_{12} + \pi_{21}
two.tailed	logical; should a two-tailed or one-tailed test be used?
correct	logical; use continuity correction? Only applicable for 2x2 tables
gen_fun	function used to generate the required discrete data. Object returned must be a matrix with k rows and k columns of counts. Default uses gen_mcnemar.test. User defined version of this function must include the argument ...
return_analysis	logical; return the analysis object for further extraction and customization?
...	additional arguments to be passed to gen_fun. Not used unless a customized gen_fun is defined

Value

a single p-value

Author(s)

Phil Chalmers [e-mail]

See Also

gen_mcnemar.test

Examples

# from ?mcnemar.test
Performance <- matrix(c(794, 86, 150, 570),
		   nrow = 2,
		   dimnames = list("1st Survey" = c("Approve", "Disapprove"),
		               "2nd Survey" = c("Approve", "Disapprove")))
(prop <- prop.table(Performance))

# one sample + test and resulting p-value
p_mcnemar.test(n=sum(Performance), prop=prop)

# return analysis model
p_mcnemar.test(n=sum(Performance), prop=prop, return_analysis=TRUE)



# post-hoc power (not recommended)
Spower(p_mcnemar.test(n=sum(Performance), prop=prop))

# odds ratio + discordant proportions supplied instead
OR <- prop[1,2] / prop[2,1]
disc <- prop[1,2] + prop[2,1]
p_mcnemar.test(n=50, OR=.25, prop.disc=disc, two.tailed=FALSE) |>
  Spower(replications=30000)

Spower documentation built on Nov. 28, 2025, 1:06 a.m.