NSM3 (version 1.1)

pBohnWolfe: Function to compute the P-value for the observed Bohn-Wolfe U statistic.

Description

This function computes the U statistic and then uses Monte Carlo sampling to compute the corresponding P-value. The Monte Carlo samples are simulated based on the order statistics of a uniform(0,1) distribution.

Usage

pBohnWolfe(x,y,k,q,c,d,method="Monte Carlo",n.mc=10000)

Value

Returns a list with "NSM3Ch5p" class containing the following components:mnumber of observations in RSS for the first data group (X)nnumber of observations in RSS for the second data group (Y)obs.statthe observed U statisticp.valupper tail P-value

References

Bohn, Lora L., and Douglas A. Wolfe. "Nonparametric two-sample procedures for ranked-set samples data." Journal of the American Statistical Association 87.418 (1992): 552-561

Examples

##Hollander, Wolfe, Chicken Example 15.4 Body Mass Index:
male<-c(18.0, 20.5, 21.3, 21.3, 22.3, 23.8, 23.8, 24.6, 25.0, 25.2, 25.3, 25.9, 26.1, 27.0,
27.4, 27.4, 28.4, 29.4, 29.6, 32.8)
female<-c(17.2, 17.8, 19.9, 20.0, 21.7, 22.0, 22.3, 23.1, 23.9, 25.8, 27.1, 29.6, 30.1, 30.3,
30.7, 31.1, 35.2, 35.6, 38.1, 42.5)

pBohnWolfe(male,female,4,4,5,5)
##To use more Monte Carlo samples:
#pBohnWolfe(male,female,4,4,5,5,n.mc=100000)