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

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

From |
Adam Olszewski <adam.olszewski@gmail.com> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: KM curve adjusted by the Inverse Probability of Treatment Weights(IPTW) |

Date |
Thu, 20 Dec 2012 12:31:35 -0500 |

Hi Marcos, In Stata, you just need to stset your dataset indicating appropriate weights: stset time [pw=iptw], fail(remission) Then you can graph anything you want, and the weight adjustment is done automatically: sts graph, by(tx) I'm not going into the intricacies of what such graphs might or might not represent, since I trust you understand them. Curve comparison can be only done with a Cox regression (or another regression model if you prefer, such as Royston-Parmar). You will not be able to use the efron or exact tie adjustment. sts test tx which should be the same as stcox test I am not a statistician and do not understand very well why/if non-parametric tests can or cannot be done with pweighted data, but simple logrank certainly cannot be done and Stata will refuse to perform it (it will default to Cox). On the other hand, the -stpepemori- (available through SSC) command from Enzo Coviello works in a pweighted dataset and I had at least one expert declare that it should be valid to use. If any of our statisticians on the list can (p)weigh in on the subject I would be delighted to hear! Great reference: Cole, S. R. and M. A. Hernan (2004). "Adjusted survival curves with inverse probability weights." Comput Methods Programs Biomed 75(1): 45-49. And of course the famous: Hernan, M. A., et al. (2000). "Marginal structural models to estimate the causal effect of zidovudine on the survival of HIV-positive men." Epidemiology 11(5): 561-570. Adam Olszewski On Thu, Dec 20, 2012 at 11:55 AM, Marcos Vinicius <bioestatistico@inbox.com> wrote: > Hello, > I would like to generate a Kaplan Meier curve adjusted(controlled/corrected) by the Inverse Probability of Treatment Weights(IPTW). > Supose I have a dataset with columns: tx=treatment (0/1) , ps=propensity score weig= Inverse Probability of Treatment Weights (iptw),time (=survival time)remission(1/0) (0=censored), > In IPTW, individuals are weighted by the inverse probability of receiving the treatment that they actually received. > In SAS I tried this code but it is not working although it is indicated in the SAS/STAT > 12.1 User’s Guide The LIFETEST Procedure (Chapter) > Someone could help me with a Stata code to produte and compare 2 Km curves according to treatment group but adjusted by iptw(weig) > proc lifetest method=km data=ps plots = (survival ); > > time time*remission(0); > weight weig; > survival out=surv_data conftype=loglog ; > /*weight weig;*/ > ODS EXCLUDE ProductLimitEstimates Means LogrankHomCov WILHOMCOV HomStats TaroneHomCov PetoHomCov Rank HomTests ; > /*weight weig;*/ > strata tx; > /*weight weig;*/ > run; > [ > Error message from SAS log > ERROR 180-322: Statement is not valid or it is used out of > proper order. > > > After program the propensity score weight , also referred to as the inverse probability oftreatment weight (IPTW)( weig in the SAS code) I tried to use the weight statatement from proc Lifetest to produce a Kaplan-Meier curve adjusted by iptw. It is not working. > Someone could help me in the code implementation or guide another way to produce the KM curve adjusted by propensity score?: > Many thanks, > Vinicius. > > ____________________________________________________________ > GET FREE SMILEYS FOR YOUR IM & EMAIL - Learn more at http://www.inbox.com/smileys > Works with AIM®, MSN® Messenger, Yahoo!® Messenger, ICQ®, Google Talk™ and most webmails > > > > * > * 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/

**Follow-Ups**:**Re: st: KM curve adjusted by the Inverse Probability of Treatment Weights(IPTW)***From:*Marcos Vinicius <bioestatistico@inbox.com>

**References**:**st: KM curve adjusted by the Inverse Probability of Treatment Weights(IPTW)***From:*Marcos Vinicius <bioestatistico@inbox.com>

- Prev by Date:
**Re: st: weight in tab1** - Next by Date:
**RE: st: split a sample** - Previous by thread:
**st: KM curve adjusted by the Inverse Probability of Treatment Weights(IPTW)** - Next by thread:
**Re: st: KM curve adjusted by the Inverse Probability of Treatment Weights(IPTW)** - Index(es):