# R/v.chi.sq.R In MOTE: Effect Size and Confidence Interval Calculator

#' V for Chi-Square
#'
#' This function displays V and non-central confidence interval
#' for the specified chi-square statistic.
#'
#' V is calculated by finding the square root of chi-squared divided by the product
#' of the sample size and the degrees of freedom with the lowest value.
#'
#'      v = sqrt(x2 / (n * dfsmall))
#'
#'
#' @param x2 chi-square statistic
#' @param n sample size
#' @param r number of rows in the contingency table
#' @param c number of columns in the contingency table
#' @param a significance level
#' @return Provides V with associated confidence intervals
#' and relevant statistics.
#'
#' \item{v}{v-statistic}
#' \item{vlow}{lower level confidence interval of omega}
#' \item{vhigh}{upper level confidence interval of omega}
#' \item{n}{sample size}
#' \item{df}{degrees of freedom}
#' \item{x2}{significance statistic}
#' \item{p}{p-value}
#' \item{estimate}{the V statistic and confidence interval in
#' APA style for markdown printing}
#' \item{statistic}{the X2-statistic in APA style for markdown printing}
#'
#' @keywords effect size, chi-square
#' @import MBESS
#' @import stats
#' @export
#' @examples
#'
#' #The following example is derived from the "chisq_data" dataset, included
#' #in the MOTE library.
#'
#' #Individuals were polled about their number of friends (low, medium, high)
#' #and their number of kids (1, 2, 3+) to determine if there was a
#' #relationship between friend groups and number of children, as we
#' #might expect that those with more children may have less time for
#' #friendship maintaining activities.
#'
#' chisq.test(chisq_data$kids, chisq_data$friends)
#'
#' v.chi.sq(x2 = 2.0496, n = 60, r = 3, c = 3, a = .05)
#'
#' #Please note, if you see a warning, that implies the lower effect should
#' #be zero, as noted.

v.chi.sq <- function (x2, n, r, c, a = .05) {

if (missing(x2)){
stop("Be sure to include chi-square statistic value.")
}

if (missing(n)){
stop("Be sure to include the sample size.")
}

if (missing(r)){
stop("Be sure to include number of rows.")
}

if (missing(c)){
stop("Be sure to include number of columns")
}

if (a < 0 || a > 1) {
stop("Alpha should be between 0 and 1.")
}

dfsmall <- min(r - 1, c - 1)
v <- sqrt(x2 / (n * dfsmall))
dftotal <- (r - 1) * (c - 1)
ncpboth <- conf.limits.nc.chisq(x2, df = dftotal, conf.level = (1 - a))
vlow <- sqrt((ncpboth$Lower.Limit + dftotal) / (n * dfsmall)) vhigh <- sqrt((ncpboth$Upper.Limit + dftotal) / (n * dfsmall))
p <- pchisq(x2, dftotal, lower.tail = F)

if (p < .001) {reportp = "< .001"} else {reportp = paste("= ", apa(p,3,F), sep = "")}

output = list("v" = v, #v stats
"vlow" = vlow,
"vhigh" = vhigh,
"n" = n, #sample stats
"df" = dftotal,
"x2" = x2, #sig stats,
"p" = p,
"estimate" = paste("$V$ = ", apa(v,2,F), ", ", (1-a)*100, "\\% CI [",
apa(vlow,2,F), ", ", apa(vhigh,2,F), "]", sep = ""),
"statistic" = paste("$\\chi^2$(", dftotal, ") = ", apa(x2,2,T), ", $p$ ",
reportp, sep = ""))

return(output)
}

#' @rdname v.chi.sq
#' @export


## Try the MOTE package in your browser

Any scripts or data that you put into this service are public.

MOTE documentation built on May 2, 2019, 5:51 a.m.