Note that
levelsof factor, local(fac_lev)
foreach i of local fac_lev {
gen ydf`i' = cdf if factor==`i'
}
is pretty much equivalent to
separate cdf, by(factor) gen(ydf)
Nick
[email protected]
[email protected]
To plot multiple KM curves on the same graph by levels of some factor,
create a separate variable for each with something like:
*******************************************
sts gen surv = s, by(factor)
gen cdf =100*(1-surv)
****************************************************
* Create a separate CDF for each factor level
****************************************************
sort factor outcome
levelsof factor, local(fac_lev)
foreach i of local fac_lev {
gen ydf`i' = cdf if factor==`i'
}
twoway scatter ydf* outcome if cdf<= 0.5 , // Add remaining options below
On Mon, Aug 31, 2009 at 11:27 AM, <[email protected]> wrote:
> --
> The only solution that I've found to this problem is to generate the
> the failure curve (1 - the survival curve), using -sts gen- and then
> to use -twoway- with an -if- clause to limit the range and a
> "connect" option to connect the points.
> -Steve
>
> On Mon, Aug 31, 2009 at 7:37 AM, Philip Ryan<[email protected]> wrote:
>> sts graph...., .... ylabel(0(.1).5, ang(hor) nogrid).....
>>
>> Phil
>>
>> -----Original Message-----
>> From: [email protected]
>> [mailto:[email protected]] On Behalf Of Fiona Bruinsma
>> Sent: Monday, 31 August 2009 10:00 PM
>> To: [email protected]
>> Subject: st: Decrease range of y axis
>>
>> Hi,
>>
>> I am creating a graph of a kaplan-meier curve. The y axis (cumulative
>> probability) ranges from 0.0 to 1.0 however all my values are less than 0.5
>> making the data difficult to see. I would like to decrease the range of the
>> y axis from 0.0-1.0 to 0.0-0.50. I have tried yscale but notice in the
>> manual that this can only be used to increase the scale. I have not been
>> able to find any reference to decreasing the range.
>>
>> The syntax I have been using is as follows:
>> sts graph, by(tx_cen) failure noorigin ytitle("Cumulative probability of
>> delivery") /*
>> */ xtitle("Gestation (weeks)") ylabel(,angle(horizontal) nogrid) /*
>> */ title(" ") /*
>> */ legend(position(3) ring(0)) legend(col(1)) clpattern("l" "_")
>> legend(region(""))/*
>> */legend(region(lstyle(none))) graphregion(fcolor(white))
>>
*
* 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/
