Chapter Contents Previous Next
 Functions

STDMED Function

computes the standard deviation of a sample median.

Syntax

STDMED(n)

where n is the sample size.

Description

The STDMED function gives the standard deviation of the median of a normally distributed sample with a mean of 0 and a variance of 1. This function gives the standard error used to determine the width of the control limits for charts produced by the MCHART and MRCHART statements in PROC SHEWHART.

Let n represent the sample size and X(i) represent the ith order statistic. Then, when n is odd, the STDMED function calculates , where

where Band is the gamma function, is the standard normal cumulative distribution function, and is the corresponding density function.

If n is even, the function calculates the square root of the following:

Var [ [( X(n/2) + X((n/2)+1) )/2] ] =
(1/4) [ E(X(n/2)2) + E(X((n/2)+1)2) + 2 E(X(n/2)X((n/2)+1)) ]

where

For more details, refer to David (1981), Kendall and Stuart (1977, 252), and Sarhan and Greenberg (1962).

Examples

These statements use a loop to calculate the standard deviation of the median for sample sizes from 6 to 12:

   data;
do n=6 to 12;
s=stdmed(n);
put s;
output;
end;
run;


The statements produce these values:

   0.4634033519
0.4587448763
0.410098592
0.4075552495
0.3719226208
0.3703544701
0.3428063408


 Chapter Contents Previous Next Top