Thanks to Kit Baum, there is now a new version of the -smileplot- package
available for download from SSC. In Stata, use the -ssc- command to do this.
The -smileplot- package is defined as below on my website. In the new
version, the default value of the -refopts()- option is now
-refopts(lstyle(xyline))-, and not -refopts(lstyle(refline))- as in the
previous version. Therefore, unless the user specifies otherwise, the X-
and Y-axis reference lines, representing the null hypothesis and the
critical P-values, will be drawn using the same style as other X- and
Y-axis reference lines, specified by the -xline()- and -yline()- options of
-graph twoway-.
I hope this change doesn't cause too much inconvenience to users. It seems
that StataCorp understand the term "reference lines" to mean something
other than X- and Y-axis reference lines, at least if you read their own
scheme files. To find more about line styles in Stata graphics, type
whelp linestyle
Best wishes
Roger
------------------------------------------------------------------------------------
package smileplot from http://www.kcl-phs.org.uk/rogernewson/stata8
------------------------------------------------------------------------------------
TITLE
smileplot: Multiple test procedures and smile plots
DESCRIPTION/AUTHOR(S)
This package contains the programs multproc, smileplot and smileplot7.
multproc inputs a data set with 1 observation for each of a set of
multiple
significance tests and data on the P-values, and carries out a
multiple test
procedure chosen by the user to define a corrected overall critical
P-value
for accepting or rejecting the null hypotheses tested. These procedures
may be one-step, step-up or step-down, and may control the
familywise error
rate (eg the Bonferroni, Sidak, Holm, Holland-Copenhaver, Hochberg
and Rom
procedures) or the false discovery rate (eg the Simes, Benjamini-Liu,
Benjamini-Yekutieli and Benjamini-Krieger-Yekutieli procedures).
smileplot,
and its Stata 7 version smileplot7, work by calling multproc and then
creating a smile plot, with data points corresponding to multiple
estimated
parameters, the P-values (on a reverse log scale) on the Y-axis, and the
corresponding parameter estimates (or another variable) on the
X-axis. There
are Y-axis reference lines at the uncorrected and corrected overall
critical
P-values. The reference line at the corrected critical P-value,
known as the
parapet line, is interpreted informally as a boundary between data
mining and
data dredging. multproc, smileplot and smileplot7 are used on data
sets with
one observation per estimated parameter and data on estimates and their
P-values, which may be created using parmby, parmest, statsby or
postfile.
Author: Roger Newson
Distribution-Date: 22november2005
Stata-Version: 8
INSTALLATION FILES (click here to install)
multproc.ado
multproc.hlp
smileplot.ado
smileplot.hlp
smileplot7.ado
smileplot7.hlp
------------------------------------------------------------------------------------
(click here to return to the previous screen)
--
Roger Newson
Lecturer in Medical Statistics
Department of Public Health Sciences
Division of Asthma, Allergy and Lung Biology
King's College London
5th Floor, Capital House
42 Weston Street
London SE1 3QD
United Kingdom
Tel: 020 7848 6648 International +44 20 7848 6648
Fax: 020 7848 6620 International +44 20 7848 6620
or 020 7848 6605 International +44 20 7848 6605
Email: [email protected]
Website: http://phs.kcl.ac.uk/rogernewson/
Opinions expressed are those of the author, not the institution.
