Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
| From | "Kabaso M.E." <Mushota.Kabaso@soton.ac.uk> |
| To | "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu> |
| Subject | RE: st: RE: Fitting distributions to right-sensored survival data |
| Date | Tue, 26 Nov 2013 15:42:31 +0000 |
Dear Maarten, Thanks for your guidance. I am able to obtain the shape parameter this way, what about the scale parameter? Mushota -----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Maarten Buis Sent: 22 November 2013 13:08 To: statalist@hsphsun2.harvard.edu Subject: Re: st: RE: Fitting distributions to right-sensored survival data On Fri, Nov 22, 2013 at 1:44 PM, Kabaso M.E. wrote: > Essentially, I can use the streg for example to fit a Weibull model to my survival data but how can I get the model's parameters for me to recreate the pdf here fitted elsewhere? The (transformed) parameters are returned in the matrix e(b). You can also use the -[]_b[]- syntax to get at specific parameters, see: -help _variables-. Here is an example: *------------------ begin example ------------------ sysuse cancer streg i.drug age, dist(weibull) nohr matlist e(b) di "p = " exp([ln_p]_b[_cons]) di "1/p = " 1/exp([ln_p]_b[_cons]) *------------------- end example ------------------- * (For more on examples I sent to the Statalist see: * http://www.maartenbuis.nl/example_faq ) Hope this helps, Maarten --------------------------------- Maarten L. Buis WZB Reichpietschufer 50 10785 Berlin Germany http://www.maartenbuis.nl --------------------------------- * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/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/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/