-------------------------------------------------------------------------------
       log:  C:\hochberg.log
  log type:  text
 opened on:  14 Mar 2014, 10:03:58


* input some raw p-values and do a hochberg adjustment

. input raw_p

         raw_p
  1. .09109
  2. .69122
  3. .00177
  4. .57181
  5. .03121
  6. .01413
  7. 
. end


. multproc, pval (raw_p) method(hochberg) crit(CV) reject(rej) puncor(.05)

Method: hochberg
Uncorrected overall critical P-value: .05
Number of P-values: 6
Corrected overall critical P-value: .00833333
Number of rejected P-values: 1

. list

     +--------------------------+
     |  raw_p          CV   rej |
     |--------------------------|
  1. | .09109   .01666667     0 |
  2. | .69122         .05     0 |
  3. | .00177   .00833333     1 |
  4. | .57181        .025     0 |
  5. | .03121       .0125     0 |
     |--------------------------|
  6. | .01413         .01     0 |
     +--------------------------+



. * reverse sort on raw_p for printing 
. gsort -raw_p

. list

     +--------------------------+
     |  raw_p          CV   rej |
     |--------------------------|
  1. | .69122         .05     0 |
  2. | .57181        .025     0 |
  3. | .09109   .01666667     0 |
  4. | .03121       .0125     0 |
  5. | .01413         .01     0 |
     |--------------------------|
  6. | .00177   .00833333     1 |
     +--------------------------+