NAG Library Routine Document
G08AGF performs the Wilcoxon signed rank test on a single sample of size .
|SUBROUTINE G08AGF (
||N, X, XME, TAIL, ZER, W, WNOR, P, N1, WRK, IFAIL)
||N, N1, IFAIL
||X(N), XME, W, WNOR, P, WRK(3*N)
The Wilcoxon one-sample signed rank test may be used to test whether a particular sample came from a population with a specified median. It is assumed that the population distribution is symmetric. The data consists of a single sample of observations denoted by . This sample may arise from the difference between pairs of observations from two matched samples of equal size taken from two populations, in which case the test may be used to test whether the median of the first population is the same as that of the second population.
The hypothesis under test,
, often called the null hypothesis, is that the median is equal to some given value
, and this is to be tested against an alternative hypothesis
- : population median ; or
- : population median ; or
- : population median ,
using a two tailed, upper tailed or lower tailed probability respectively. You select the alternative hypothesis by choosing the appropriate tail probability to be computed (see the description of parameter TAIL
in Section 5
The Wilcoxon test differs from the Sign test (see G08AAF
) in that the magnitude of the scores is taken into account, rather than simply the direction of such scores.
The test procedure is as follows
||For each , for , the signed difference is found, where is a given test value for the median of the sample.
||The absolute differences are ranked with rank and any tied values of are assigned the average of the tied ranks. You may choose whether or not to ignore any cases where by removing them before or after ranking (see the description of the parameter ZER in Section 5).
||The number of nonzero is found.
||To each rank is affixed the sign of the to which it corresponds. Let .
||The sum of the positive-signed ranks, , is calculated.
||the test statistic ;
||the number of nonzero ;
||the approximate Normal test statistic , where
||the tail probability, , corresponding to , depending on the choice of the alternative hypothesis, .
If , is computed exactly; otherwise, an approximation to is returned based on an approximate Normal statistic corrected for continuity according to the tail specified.
The value of can be used to perform a significance test on the median against the alternative hypothesis. Let be the size of the significance test (that is, is the probability of rejecting when is true). If then the null hypothesis is rejected. Typically might be or .
Conover W J (1980) Practical Nonparametric Statistics Wiley
Neumann N (1988) Some procedures for calculating the distributions of elementary nonparametric teststatistics Statistical Software Newsletter 14(3) 120–126
Siegel S (1956) Non-parametric Statistics for the Behavioral Sciences McGraw–Hill
- 1: N – INTEGERInput
On entry: , the size of the sample.
- 2: X(N) – REAL (KIND=nag_wp) arrayInput
On entry: the sample observations, .
- 3: XME – REAL (KIND=nag_wp)Input
On entry: the median test value, .
- 4: TAIL – CHARACTER(1)Input
: indicates the choice of tail probability, and hence the alternative hypothesis.
- A two tailed probability is calculated and the alternative hypothesis is : population median .
- An upper tailed probability is calculated and the alternative hypothesis is : population median .
- A lower tailed probability is calculated and the alternative hypothesis is : population median .
, or .
- 5: ZER – CHARACTER(1)Input
: indicates whether or not to include the cases where
in the ranking of the
- All are included when ranking.
- All , are ignored, that is all cases where are removed before ranking.
- 6: W – REAL (KIND=nag_wp)Output
On exit: the Wilcoxon rank sum statistic, , being the sum of the positive ranks.
- 7: WNOR – REAL (KIND=nag_wp)Output
: the approximate Normal test statistic,
, as described in Section 3
- 8: P – REAL (KIND=nag_wp)Output
: the tail probability,
, as specified by the parameter TAIL
- 9: N1 – INTEGEROutput
On exit: the number of nonzero 's, .
- 10: WRK() – REAL (KIND=nag_wp) arrayWorkspace
- 11: IFAIL – INTEGERInput/Output
must be set to
. If you are unfamiliar with this parameter you should refer to Section 3.3
in the Essential Introduction for details.
For environments where it might be inappropriate to halt program execution when an error is detected, the value
is recommended. If the output of error messages is undesirable, then the value
is recommended. Otherwise, if you are not familiar with this parameter, the recommended value is
. When the value is used it is essential to test the value of IFAIL on exit.
unless the routine detects an error or a warning has been flagged (see Section 6
6 Error Indicators and Warnings
If on entry
, explanatory error messages are output on the current error message unit (as defined by X04AAF
Errors or warnings detected by the routine:
|On entry,||, or .|
|or|| or .|
The whole sample is identical to the given median test value.
The approximation used to calculate when will return a value with a relative error of less than for most cases. The error may increase for cases where there are a large number of ties in the sample.
The time taken by G08AGF increases with , until , from which point on the approximation is used. The time decreases significantly at this point and increases again modestly with for .
This example performs the Wilcoxon signed rank test on two matched samples of size , taken from two populations. The distribution of the differences between pairs of observations from the two populations is assumed to be symmetric. The test is used to test whether the medians of the two distributions of the populations are equal or not. The test statistic, the approximate Normal statistic and the two tailed probability are computed and printed.
9.1 Program Text
Program Text (g08agfe.f90)
9.2 Program Data
Program Data (g08agfe.d)
9.3 Program Results
Program Results (g08agfe.r)