Thanks to Kit Baum, a new package -qqvalue- is now available for download from SSC. In Stata, use the -ssc- command to do this.
The -qqvalue- package is described as below on my website. It should be useful for users of multiple-test procedures, such as genome scanners, epigenome scanners, or even FFQnome scanners, who measure multiple associations between diseases and items from food frequency questionnaires (FFQs).
Best wishes
Roger
Roger B Newson BSc MSc DPhil
Lecturer in Medical Statistics
Respiratory Epidemiology and Public Health Group
National Heart and Lung Institute
Imperial College London
Royal Brompton Campus
Room 33, Emmanuel Kaye Building
1B Manresa Road
London SW3 6LR
UNITED KINGDOM
Tel: +44 (0)20 7352 8121 ext 3381
Fax: +44 (0)20 7351 8322
Email: [email protected]
Web page: http://www.imperial.ac.uk/nhli/r.newson/
Departmental Web page:
http://www1.imperial.ac.uk/medicine/about/divisions/nhli/respiration/popgenetics/reph/
Opinions expressed are those of the author, not of the institution.
-------------------------------------------------------------------------------
package qqvalue from http://www.imperial.ac.uk/nhli/r.newson/stata10
-------------------------------------------------------------------------------
TITLE
qqvalue: Generate quasi-q-values by inverting multiple-test procedures
DESCRIPTION/AUTHOR(S)
qqvalue is similar to the R package p.adjust. It inputs a single
variable, assumed to contain P-values calculated for multiple
comparisons, in a dataset with 1 observation per comparison. It
outputs a new variable, containing the quasi-q-values
corresponding to these P-values, calculated by inverting a
nultiple-test procedure specified by the user. These
quasi-q-values represent, for each corresponding P-value, the
minimum uncorrected P-value threshold for which that P-value
would be in the discovery set, assuming that the specified
multiple-test procedure was used on the same set of input
P-values to generate a corrected P-value threshold. These
minimum uncorrected P-value thresholds may represent familywise
error rates or false discovery rates, depending on the procedure
used. Optionally, qqvalue may output other variables, containing
the various intermediate results used in calculating the
quasi-q-values. The multiple-test procedures available for
qqvalue are a subset of those available using the multproc module
of the smileplot package, which can be downloaded from SSC.
Author: Roger Newson
Distribution-Date: 08december2009
Stata-Version: 10
INSTALLATION FILES (click here to install)
qqvalue.ado
qqvalue.sthlp
-------------------------------------------------------------------------------
(click here to return to the previous screen)
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
