DIFFERENCE OF QUANTILE

DIFFERENCE OF QUANTILE (LET)

Let Subcommand

Compute the difference between the QUANTILEs for two response variables.

The qth quantile of a data set is defined as that value where a q fraction of the data is below that value and (1-q) fraction of the data is above that value. For example, the 0.5 quantile is the median.

Dataplot supports two methods for computing the quantile.

The first method is based on the order statistic. The formula is:

where

X are the observations sorted in ascending order
NI1 = INT(q*(n+1))
NI2 = NI1 + 1
r = q*(n+1) - INT(q*(n+1))

An alternative method is called the Herrell-Davis estimate. This method attempts to provide a lower standard error for X_q by utilizing all the order statistics rather than a single (or a weighted average of two) order statistic. Note that there are caes where the Herrell-Davis has a substantially smaller standard error than the order statistic method. However, there are also cases where the reverse is true.

To compute the Herrell-Davis estimate, do the following:

1. Sort the X in ascending order.

2. A = (n+1)*q

3. B = (n+1)*(1 - q)

4. W_i = BETCDF(i/n,A,B) - BETCDF((i-1)/n,A,B) where BETCDF is the beta cumulative distribution function with shape parameters A and B.

5. 

For the difference of quantiles, the quantile is computed for each of two samples then their difference is taken.

LET <par> = DIFFERENCE OF QUANTILE <y1> <y2>
                        <SUBSET/EXCEPT/FOR qualification>
where <y1> is the first response variable;
            <y2> is the first response variable;
            <par> is a parameter where the computed difference of the quantiles is stored;
and where the <SUBSET/EXCEPT/FOR qualification> is optional.

LET A = DIFFERENCE OF QUANTILE Y1 Y2
LET A = DIFFERENCE OF QUANTILE Y1 Y2 SUBSET X > 1

The following command is used to determine which method is used to compute the quantile:
SET QUANTILE METHOD <ORDER/HERRELL-DAVIS>

The specific quantile to compute is specified by entering the following command:
LET XQ = <value>

where <value> is a number in the interval (0,1) that specifies the desired quantile.

The difference of QUANTILEs statistic is also supported for the following commands:
DIFFERENCE OF QUANTILE PLOT Y1 Y2 X
CROSS TABULATE DIFFERENCE OF QUANTILE PLOT Y1 Y2 X
BOOTSTRAP DIFFERENCE OF QUANTILE PLOT Y1 Y2
TABULATE DIFFERENCE OF QUANTILE Y1 Y2 X
CROSS TABULATE DIFFERENCE OF QUANTILE Y1 Y2 X
LET Z = CROSS TABULATE DIFFERENCE OF QUANTILE Y1 Y2 X1 X2


Enter

HELP STATISTIC PLOT
HELP CROSS TABULATE PLOT
HELP BOOTSTRAP PLOT
HELP TABULATE
HELP CROSS TABULATE
HELP CROSS TABULATE (LET)


for more information on these plots and commands.

None

None

QUANTILE	= Compute a specified quantile.
MINIMUM	= Compute the minimum.
MAXIMUM	= Compute the maximum.
DIFFERENCE OF MINIMUM	= Compute the difference of the minimums.
DIFFERENCE OF MAXIMUM	= Compute the difference of the maximums.
STATISTICS PLOT	= Generate a statistic versus subset plot.
CROSS TABULATE PLOT	= Generate a statistic versus two subsets plot.
BOOTSTRAP PLOT	= Generate a bootstrap plot.
TABULATE	= Perform a tabulation for a specified statistic.
CROSS TABULATE	= Perform a cross tabulation for a specified statistic.

"Introduction to Robust Estimation and Hypothesis Testing", Rand Wilcox, Academic Press, 1997.

"A New Distribution-Free Quantile Estimator", Frank Herrell and C. E. Davis, Biometrika (1982), 69(3), 635-640.

Data Analysis

3/2003

SKIP 25 
READ IRIS.DAT Y1 TO Y4 X 
. 
LET XQ = 0.9 
. 
LET A = DIFFERENCE OF QUANTILE Y1 Y2 
TABULATE DIFFERENCE OF QUANTILE Y1 Y2 X 
. 
XTIC OFFSET 0.2 0.2 
X1LABEL GROUP ID 
Y1LABEL DIFFERENCE OF (0.9) QUANTILE 
CHAR X 
LINE BLANK 
DIFFERENCE OF QUANTILE PLOT Y1 Y2 X 
CHAR X ALL 
LINE BLANK ALL 
BOOTSTRAP DIFFERENCE OF QUANTILE PLOT Y1 Y2 X  
    

Dataplot generated the following output.