qvalue [options] <p-values file>
Takes as input a list of p-values and prints a corresponding list of q-values, computed using the method of Benjamini and Hochberg (Journal of the Royal Statistical Society B, 57:289-300, 1995).
The name of a tab separated file containing a column of p-values. Other columns may exist. If the name is - then the input data is read from standard input.
Outputs tab-delimited columns. If the --append option is given, then the line from the input is listed first, followed by the p-value (only when --null is specified) and q-value. Otherwise it outputs the raw score (only when --null is specified). followed by the p-value and q-value. The lines are sorted by p-value.
| Option | Parameter | Description | Default Behaviour |
|---|---|---|---|
| General Options | |||
| --null | file | Treat the primary input file as raw scores, and begin by estimating p-values from the empirical null scores given in this file. The --header and --column arguments apply to this file as well as the primary input file. | The primary input is assumed to actually contain p-values. |
| --good-score | high|low | Choose low if smaller scores are better or high if larger scores are better. | Smaller scores are considered better. |
| --pi-zero | Estimate the fraction (pi-zero) of the empirical distribution that is
drawn according to the null hypothesis. This uses the
method described in:
John Storey, "A direct approach to false discovery rates",
Journal of the Royal Statistical Society B, 64:3, 479-498,
2002.
|
The estimation step is skipped and the program conservatively sets pi-zero to 100%. | |
| --pi-zero-file | file name | Specify a file to store the pi-zero estimate in. This option implies the --pi-zero option. | The results of the pi-zero estimate (if calculated) are not stored in a file. |
| --fdr | Only compute false discovery rate (FDR) estimates. | The FDR values are converted into q-values. | |
| --bootstraps | num | Number of bootstraps to perform when estimating pi-zero. | The program performs 1000 bootstraps. |
| --header | line count | Treat the first line count lines of the p-values file as the header. If the --append option is specified these lines will also be included in the output. | No lines are treated as headers. |
| --column | column number | Read the p-value from the specified column number. The first column is 1, the second column is 2 and so on... | The p-value is read from the first column. |
| --append | Print all the columns in the input file to standard output with additional p-value (raw scores only) and q-value columns. | Print the raw score (only when raw scores are supplied), p-value and q-value in tab separated columns. | |
| --seed | seed | Set the seed value for the pseudo-random number generator. | The seed value is set from the system clock. |