| |
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]
st: graph comparing results of several models after -parmby-
From |
Michael McCulloch <[email protected]> |
To |
Statalist <[email protected]> |
Subject |
st: graph comparing results of several models after -parmby- |
Date |
Tue, 17 Oct 2006 11:30:00 -0700 |
Hello,
MY PROJECT:
I am comparing several different regression models whose results I have
saved using -parmby-, and then concatenated into one data set using dsconcat.
WHAT I HAVE DONE:
*run models, adjusted and unadjusted (thanks to Roger Newson's paper)
tempfile tf1 tf2
parmby "stcox tcm", lab saving (tf1, replace) idn(1)
ids(Unadjusted) eform
parmby "stcox tcm stage", lab saving (tf2, replace) idn(2)
ids(Adjusted) eform
dsconcat tf1 tf2
format estimate min* max* %8.2f
format p %8.2e
and when* forest plot comparing results (thanks to Scott Merryman)
encode ids, gen(Model)
twoway rcap min95 max95 Model if parmseq==1, horizontal
msize(huge) ///
xtitle(Hazard Ratio and 95% Confidence Interval, margin(0 0 0 5)) ///
ytitle(Analysis Method) ylabel( 1 "Adjusted" 2 "Unadjusted",
angle(h)) ///
|| scatter Model estimate if parmseq==1, msize(medlarge) legend(off)
WHAT I WOULD LIKE TO ACHIEVE:
How would I modify the rcap/scatter command so that I could adjust the
vertical distance between plotted whisker bars?
How would I modify the ylabel specification so I could automate the
labeling without manually typing "adjusted", "unadjusted", etc. (I'm
planning to compare many different models.)
Best wishes,
Michael McCulloch
____________________________________
Michael McCulloch
Pine Street Foundation
124 Pine Street, San Anselmo, CA 94960-2674
tel 415.407.1357
fax 415.485.1065
email: [email protected]
web: www.pinestreetfoundation.org
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/