Re: sts test: displaying log-rank test of equality within a Kaplan Meier graph

[Michael McCulloch <mm@pinestreetfoundation.org>]

Thank you kindly Michael!

On Feb 11, 2010, at 8:48 PM, Michael Norman Mitchell wrote:

> Dear Michael
>
>  If I understand your last message, you want to be able to change  
> the placement of the "log rank P= 3.67" text. If so, you can control  
> that by the coordinates specified in the "text()" option. For  
> example, if you do this...
>
> sysuse cancer.dta, clear
> stset studytime
> drop if drug==3
> sts test drug
> local chi_res  logrank P = `=round(r(chi2), .01)'
> sts graph, by(drug) risktable text(.95 28 "`chi_res'", size(small)  
> col(red))
>
>  The text will be displayed in the upper right corner, because the  
> coordinates are y=.95 x=28. You can fiddle with those (Y,X) values  
> to place the text whereever yo like. I hope that helps,
>
> Michael N. Mitchell
> See the Stata tidbit of the week at...
> http://www.MichaelNormanMitchell.com
> Visit me on Facebook at...
> http://www.facebook.com/MichaelNormanMitchell
>
> On 2010-02-11 8.28 PM, Michael McCulloch wrote:
> > Thanks Phillip! This gives me the idea. May I ask:
> >    * regardless of what place option I use: ne, se, nw, sw, etc.,  
> > they all put the text box in the lower left corner. I've read -help  
> > compassdirstyle- and
> >    * bexpand also doesn't expand the size of the box.
> > Michael
> >
> >
> > On Feb 11, 2010, at 8:05 PM, Philip Ryan wrote:
> >
> > > Michael
> > >
> > > I note you wanted a "window".  Here is a box instead. I also  
> > > correct my
> > > misuse of the return scalars :-(
> > >
> > >
> > > sysuse cancer.dta, clear
> > > stset studytime
> > > drop if drug==3
> > > sts test drug
> > > local chi_P  logrank P = `=round(chi2tail(r(df),r(chi2)), .01)'
> > > sts graph, by(drug) risktable text(.05 5 "`chi_P'", place(ne) box
> > > size(small) col(red))
> > >
> > >
> > >
> > > Philip Ryan
> > > University of Adelaide
> > >
> > >
> > > -----Original Message-----
> > > From: owner-statalist@hsphsun2.harvard.edu
> > > [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Michael  
> > > McCulloch
> > > Sent: Friday, 12 February 2010 1:23 PM
> > > To: Statalist Statalist
> > > Subject: sts test: displaying log-rank test of equality within a  
> > > Kaplan
> > > Meier graph
> > >
> > > Dear Statalist members,
> > > I'm creating a Kaplan-Meier graph, and also performing the log-rank
> > > test of equality, for example:
> > >
> > > sysuse cancer.dta, clear
> > > stset studytime
> > > drop if drug==3
> > > sts graph, by(drug) risktable
> > > sts test drug
> > >
> > > My question: is it possible to display the results of the log-rank
> > > test in a small window within the area of the Kaplan-Meier survival
> > > curve?
> > >
> > >
> > > Best wishes,
> > >
> > > Michael McCulloch
> > > Pine Street Foundation
> > > 124 Pine Street
> > > San Anselmo, CA 94960-2674
> > > tel:    415-407-1357
> > > fax:     206-338-2391
> > >
> > > *
> > > *   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/
> > >
> > > *
> > > *   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/
> >
> >
> > Best wishes,
> >
> > Michael McCulloch
> > Pine Street Foundation
> > 124 Pine Street
> > San Anselmo, CA 94960-2674
> > tel:    415-407-1357
> > fax:     206-338-2391
> >
> > *
> > *   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/
> *
> *   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/


Best wishes,

Michael McCulloch
Pine Street Foundation
124 Pine Street
San Anselmo, CA 94960-2674
tel:	415-407-1357
fax: 	206-338-2391

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