FourierStuff/R/chisquare.test.R

45 lines
1.6 KiB
R
Raw Normal View History

2014-02-20 01:35:43 +00:00
# Pearson's Chi-Squared test based on the characteristic function
#
# Author: François Pelletier
###############################################################################
#' Pearson's Chi-Squared test based on the characteristic function
#' @param DATA.hist histogram object of the data
2014-02-22 18:27:24 +00:00
#' @param FUN Characteristic function (integral) or
#' Saddlepoint distribution approximation (saddlepoint)
#' @param ... FUN arguments
2014-02-20 01:35:43 +00:00
#' @param alpha tolerance level
2014-02-22 18:27:24 +00:00
#' @param method Method to approximate the distribution function. "integral" or "saddlepoint"
2014-02-20 01:35:43 +00:00
#'
#' @return A list containing the chi-square statistic,
#' degree of freedom, hypothesis reject boolean and p.value
#' @export chisquare.test
2014-02-20 01:35:43 +00:00
#' @author François Pelletier
2014-02-22 18:27:24 +00:00
chisquare.test <- function(DATA.hist,FUN,...,alpha=0.05,method="integral")
2014-02-20 01:35:43 +00:00
{
# Compute expected values for each histogram breaks using the characteristic function
2014-02-22 18:27:24 +00:00
classes <- DATA.hist$breaks
if(method=="integral")
{
expected <- diff(cftocdf(classes,FUN,...)*
sum(observed <- DATA.hist$counts))
}
else if(method=="saddlepoint")
{
expected <- diff(FUN(classes,...)*
sum(observed <- DATA.hist$counts))
}
2014-02-20 01:35:43 +00:00
# Compute the test statistic using chi-square distribution
p.value <- pchisq(chisquare.stat<-sum((observed-expected)^2/expected),
df<-length(classes)-2,lower.tail=FALSE)
# Print output
cat("Chi-Square Test based on CF\n\nTest statistic: ",chisquare.stat,
"\nDegree of freedom: ",df,
"\nP-value: ",p.value,
"\nReject H0 with confidence level ",1-alpha,"?: ",p.value<alpha)
2014-02-20 01:35:43 +00:00
# Create the return list
list(chisquare.stat=chisquare.stat,df=df,reject=p.value<alpha,p.value=p.value)
}