Bookmark and Share

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]

st: RE: RE: Summarising multiple visit data for each individual participant


From   Matthew Hurley <Matthew.Hurley@nottingham.ac.uk>
To   "statalist@hsphsun2.harvard.edu" <statalist@hsphsun2.harvard.edu>
Subject   st: RE: RE: Summarising multiple visit data for each individual participant
Date   Fri, 2 Mar 2012 16:31:32 +0000

Thank you Nick,
This was exactly what I was trying (and miserably failing) to do.
I'm very grateful for your help,
Matt

-----Original Message-----
From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Nick Cox
Sent: 02 March 2012 15:02
To: 'statalist@hsphsun2.harvard.edu'
Subject: st: RE: Summarising multiple visit data for each individual participant

... if chrpo[1]== chrpo[2]== chrpo[3]

is legal but it doesn't do what you seem to expect. I think what Stata does is evaluate from left to right, so 

chrpo[1]== chrpo[2]

is evaluated as true (1) or false (0) -- call this <result> 

and then

<result> == chrpo[3] 

is evaluated. Only exceptionally will this coincide with what you want. However, even with the expression 

chrpo if chrpo[1] == chrpo[2] & chrpo[2] == chrpo[3] 

this just copies chrpo for all observations if, and only if, values for the first three observations are identical. 

I think you want either the constant dose, when there is one, 

by ID ( visitno), sort: gen sumabx = chrpo[1] if chrpo[1]== chrpo[2] & chrpo[2] == chrpo[3] 

although clearly [2] or [3] would work as well, or an indicator 

by ID ( visitno), sort: gen constantdose = chrpo[1]== chrpo[2] & chrpo[2] == chrpo[3]

Nick 
n.j.cox@durham.ac.uk 

Matthew Hurley

I have patient registry data for which visit data is grouped by patient ID.  At each visit individuals are administered an antibiotic.  I wish to determine if specific antibiotic usage in this dataset is associated with infection in a survival analysis.

I have completed the survival analysis assuming that the antibiotic taken on the first visit represents all antibiotic taken throughout the period, however this assumption is incorrect as antibiotic usage changes on some occasions.

Is there a way I can 'summarise' each individual patients antibiotic use please?  The most pragmatic solution I can think of is to perform a sensitivity analysis for those patients whose antibiotic usage does not change (e.g. for the first 3 visits) to see if this changes my results.

I have tried by generating an indicator variable 'visitno'. 'chrpo' is my categorical variable that reports the antibiotic taken at each visit - generating a 'sumabx' variable that reports the consistently administered antibiotic for each patient.

by ID ( visitno), sort: gen sumabx = chrpo if chrpo[1]== chrpo[2]== chrpo[3] ....however this doesn't work.
I have a separate indicator variable that selects the first visit per patient.


*
*   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/
This message and any attachment are intended solely for the addressee and may contain confidential information. If you have received this message in error, please send it back to me, and immediately delete it.   Please do not use, copy or disclose the information contained in this message or in any attachment.  Any views or opinions expressed by the author of this email do not necessarily reflect the views of the University of Nottingham.

This message has been checked for viruses but the contents of an attachment
may still contain software viruses which could damage your computer system:
you are advised to perform your own checks. Email communications with the
University of Nottingham may be monitored as permitted by UK legislation.
*
*   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   |   Site index