|
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: sts test: displaying log-rank test of equality within a Kaplan Meier graph
From |
Michael Norman Mitchell <[email protected]> |
To |
[email protected] |
Subject |
Re: sts test: displaying log-rank test of equality within a Kaplan Meier graph |
Date |
Thu, 11 Feb 2010 20:48:59 -0800 |
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: [email protected]
[mailto:[email protected]] 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/