* 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 |
+--------------------------+
