Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down at the end of May, and its replacement, **statalist.org** is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
Nick Cox <njcoxstata@gmail.com> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: RE: RE: Plotting survival curves after multiple imputation |

Date |
Wed, 20 Feb 2013 11:59:14 +0000 |

It's a very small side-issue, and I can't illuminate your main question, but foreach i in 5 10 15 20 25 30 35 40 45 { replace age = . in `i' } could be re-written as forval i = 5(5)45 { replace age = . in `i' } or as replace age = . if mod(_n, 5) == 0 or as replace age = . if !mod(_n, 5) The shortest code here is also, unfortunately, the most devious and obscure. Nick On Wed, Feb 20, 2013 at 11:06 AM, Robert Long <W.R.Long@leeds.ac.uk> wrote: > Excuse me, there was a line missing in the code. Here is the correct listing: > > webuse drugtr, clear > > foreach i in 5 10 15 20 25 30 35 40 45 { > > replace age = . in `i' > > } > > > mi set wide > mi register impute age > > ice age studytime died drug, m(5) saving(drugtr_imp, replace) > > use drugtr_imp.dta, clear > > * no problem obtaining pooled estimates, and here I use the storebv option > * in an attempt to give -stcurve- what it needs mim , storebv: stcox age drug > mim, storebv: stcox age drug > > > * but it doesn't work: > stcurve, survival at1(drug=0) at2(drug=1) > > > > > -----Original Message----- > From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Robert Long > Sent: 20 February 2013 09:47 > To: 'statalist@hsphsun2.harvard.edu' > Subject: st: RE: Plotting survival curves after multiple imputation > > I have attempted to solve this using the following code, which is reproducible, but relies on the -ice- and -mim- commands > > webuse drugtr, clear > > * create some missing value > foreach i in 5 10 15 20 25 30 35 40 45 { > replace age = . in `i' > } > > mi set wide > mi register impute age > > ice age studytime died drug, m(5) saving(drugtr_imp, replace) > > use drugtr_imp.dta, clear > > * no problem obtaining pooled estimates, and here I use the storebv option > * in an attempt to give -stcurve- what it needs mim , storebv: stcox age drug > > * but it doesn't work: > stcurve, survival at1(drug=0) at2(drug=1) > > The error is > last estimates not found > r(301); > > This error seems to be produced in -stcurve- which is looking for e(cmd2) to be set to "stcox", but e(cmd2) is not set at all - e(cmd) is set to "stcox". However, running stcox normally (ie without -mim- ) , e(cmd2) /is/ set to "stcox".. > > Can anyone help ? > > Thanks > Robert Long > > > > -----Original Message----- > From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Robert Long > Sent: 15 February 2013 11:55 > To: 'statalist@hsphsun2.harvard.edu' > Subject: st: Plotting survival curves after multiple imputation > > Hello Statalist > > This sounds like something very simple, or at least it should be something that is commonly required, but so far I haven't found an answer. > > I want to use multiple imputation for missing covariate data. No problems with that. > > Having obtained m different datasets, I want to run cox regressions, obtaining pooled results. No problems with that. > > But I also want to plot a "pooled" survival curve for different levels of a covariate. > > How do I obtain such plots ? > > Thanks > Robert Long > > > * > * 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/ > > * > * 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/

**References**:**st: test of coefficients of the same regression equation***From:*Nahla Betelmal <nahlaib@gmail.com>

**Re: st: test of coefficients of the same regression equation***From:*Florian Wendelspiess Chávez Juárez <florian@chavezjuarez.com>

**Re: st: test of coefficients of the same regression equation***From:*Nahla Betelmal <nahlaib@gmail.com>

**Re: st: test of coefficients of the same regression equation***From:*Florian Wendelspiess Chávez Juárez <florian@chavezjuarez.com>

**Re: st: test of coefficients of the same regression equation***From:*Nahla Betelmal <nahlaib@gmail.com>

**st: Plotting survival curves after multiple imputation***From:*Robert Long <W.R.Long@leeds.ac.uk>

**st: RE: Plotting survival curves after multiple imputation***From:*Robert Long <W.R.Long@leeds.ac.uk>

**st: RE: RE: Plotting survival curves after multiple imputation***From:*Robert Long <W.R.Long@leeds.ac.uk>

- Prev by Date:
**Re: st: comparing equality of coefficients from two subsamples** - Next by Date:
**RE: st: Capturing unique actors meeting conditions by focal actor and year (excluding focal actor)** - Previous by thread:
**st: RE: RE: Plotting survival curves after multiple imputation** - Next by thread:
**st: trouble with traveltime command** - Index(es):