mr.dimensionality.test: Multiple-response dimensionality test

View source: R/mr.dimensionality.test.R

mr.dimensionality.testR Documentation

Multiple-response dimensionality test

Description

Performs a multiple-response dimensionality test as defined in Mahieu, Schlich, Visalli, and Cardot (2021) using random permutations to estimate the null distribution

Usage

mr.dimensionality.test(data, nperm = 2000, alpha = 0.05)

Arguments

data

A data.frame of observations in rows whose first column is a factor (the categories) and subsequent columns are binary numeric or integer, each column being a response option

nperm

Number of permuted datasets to estimate the distribution of the statistic under the null hypothesis. See details

alpha

The alpha risk of the test

Details

  • nperm: The distribution of the statistic under the null hypothesis of no associations between categories and response options is estimated using nperm datasets generated thanks to random permutations of the response vectors along observations.

Value

A list with the following elements:

dim.sig

The number of significant dimensions

statistics

Observed multiple-response chi-square statistic of each dimension

p.values

P-value of the test of each dimension adjusted for closed testing procedure

References

Loughin, T. M., & Scherer, P. N. (1998). Testing for Association in Contingency Tables with Multiple Column Responses. Biometrics, 54(2), 630-637.

Mahieu, B., Schlich, P., Visalli, M., & Cardot, H. (2021). A multiple-response chi-square framework for the analysis of Free-Comment and Check-All-That-Apply data. Food Quality and Preference, 93.

Examples

nb.obs=200
nb.response=5
nb.category=5
vec.category=paste("C",1:nb.category,sep="")
right=matrix(rbinom(nb.response*nb.obs,1,0.25),nb.obs,nb.response)
category=sample(vec.category,nb.obs,replace = TRUE)
dset=cbind.data.frame(category,right)
dset$category=as.factor(dset$category)


mr.dimensionality.test(dset)


MahieuB/MultiResponseR documentation built on June 22, 2024, 8:08 a.m.