math/

bootstrap_median.pro

top source bootstrap_median

Math

result = bootstrap_median(values, nboot=nboot, conflimit=conflimit, uniqlist=uniqlist)

Calculates the median and a confidence limit on the median based on bootstrap resampling.

Return value

Returns a 3-element vector containing the lower limit, median, and upper limit.

Parameters

values in: A vector of values whose median and error is to be calculated.

Keywords

nboot in: Number of bootstrap resamplings. Default: 1000.
conflimit in: Confidence limit. Default: 0.68 (equivalent to 1sigma for a normal distribution).
uniqlist in: If independent points are associated with more than one element of Values, then they should all be included or excluded together in the bootstrap resampling. In this case, set UNIQLIST to the result of running UNIQ on a list with the same length as Values containing the unique identifier associated with each. Note that for this to work, Values must be sorted in order of the identifier.

Examples

    Calculates the error in the median of 5000 values distributed normally:

    IDL> vals = 2.5*RANDOMN(seed,5000)
    IDL> PRINT, BOOTSTRAP_MEDIAN(vals)
         -0.25968859     -0.15505694     0.095240064

Author information

History:: Written by: Jeremy Bailin 12 June 2008 Public release in JBIU

File attributes

Modification date:	Wed Apr 15 16:07:53 2009
Lines:	82

Documentation for jbiu/

Generated by IDLdoc