Statalist


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

Re: st: Re: limiting observations within a program


From   "Sergiy Radyakin" <serjradyakin@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: Re: limiting observations within a program
Date   Thu, 2 Oct 2008 20:06:34 -0400

And where is varlist X defined?


On Thu, Oct 2, 2008 at 7:56 PM, Michael McCulloch <mm@pinest.org> wrote:
> Thanks Martin. This is my program. It first estimates probability of
> receiving treatment (called "tcm"), and then weights the data for -stcox-.
> There was one observation in which the weight was missing, I'm attempting to
> include a line in my program to drop that one observation.
>
> capture program drop msmcox_stageALL_includeIPTWd
>        program define msmcox_stageALL_includeIPTWd, rclass
>        tempname p_tcm_ALL              p_notcm_ALL     p_ALL        p_noALL
>         wt_stab_ALL
>
>        * fit the full model used in Traditional Cox
>        logit tcm dxage sex celltype txsurg txrad txchemo
>
>        *estimate probability of treatment
>        predict `p_tcm_ALL'
>        gen `p_notcm_ALL' = 1-`p_tcm_ALL' if tcm==0     //want p(tcm) for
> both users and nonusers
>        replace `p_tcm_ALL' =`p_notcm_ALL' if tcm==0
>
>        *create stabilized weight==P(A)/P(A|W)
>        logit tcm        predict `p_ALL'
>        gen `p_noALL'           =       1-`p_ALL' if tcm==0
>        replace `p_ALL' =       `p_noALL' if tcm==0
>        gen `wt_stab_ALL'       =       `p_ALL'/`p_tcm_ALL'
>
>        * stset the data
>        stset lastcontact [`iweight=wt_stab_ALL'], failure(failed)
> origin(datedx) scale(30.4375)   /*wt: stabilized IPTW*/
>
>        * cox
>                stcox tcm if `wt_stab_ALL'!=.
>        indeplist, local
>        foreach var of varlist `X' {
>                return scalar `var' = exp(_b[`var'])
>        }
> end set seed 12358
> eststo MSM_t_I: bootstrap tcm=r(tcm), eform reps(10):
> msmcox_stageALL_includeIPTWd
>
>
>
>> That is too little info. How does `var1' come into existence? A
>> -bootstrap- would cycle through the -program- several times, so that is an
>> important issue. Also -set trace on- and let us know the output...
>>
>> HTH
>> Martin
>> _______________________
>> ----- Original Message ----- From: "Michael McCulloch" <mm@pinest.org>
>> To: "Statalist" <statalist@hsphsun2.harvard.edu>
>> Sent: Friday, October 03, 2008 1:10 AM
>> Subject: st: limiting observations within a program
>>
>>> Hello,
>>> I've written a short program which estimates a probability of censoring,
>>> and have a question. The program works as written, but if I add the
>>> following line:
>>>
>>> drop if `var1'==.
>>>
>>> the program returns error message:
>>>
>>> "varlist required an error occurred when bootstrap executed program1
>>> r(100);"
>>>
>>> Any assistance would be appreciated!
>>>
>>>
>>> --
>>>
>>> Best wishes,
>>> Michael McCulloch
>>>
>>>
>>>
>>> Pine Street Foundation
>>> 124 Pine St., San Anselmo, CA 94960-2674
>>> Tel: (415) 407-1357
>>> Fax: (415) 485-1065
>>> mcculloch@pinestreetfoundation.org
>>> www.pinestreetfoundation.org
>>> *
>>> *   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/
>
> *
> *   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/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index