| Svyset in Stata for HIV results and ART use data [message #637] |
Thu, 18 July 2013 10:20  |
petermacpherson
Messages: 1
Registered: July 2013
Location: UK
|
Member |
|
|
Hello,
I am conducting a study investigating the prevalence of HIV-negative status on HIV testing among individuals who report that they are currently taking antiretrovirals. I plan to conduct this analysis using data from DHS surveys in sub-Saharan Africa, but have started with the Malawi 2010 DHS to get a feel for the data. I am using Stata 12.1 for analysis.
So far, with the Malawi 2010 DHS, I have appended the men's individual records to the women's individual records and merged with the HIV testing records, following instructions in a previous forum post.
However, I am not completely sure about how to procedure from here. I am stuck on three main areas:
1) I am not sure about how to correctly "svyset" the data in Stata. This is my code:
gen psu=v021
egen strata = group(v024 v025), label
gen weight=hiv05/1000000
svyset psu [pweight=weight], strata(strata)
Is this correct? I am relatively unfamiliar with analysing survey data.
2) Estimating the prevalence of HIV-negativity among individuals who report taking ART. Again, I am not sure if this is correct syntax
svy: tab hiv03 s1318
3) Finally, I want to construct a logistic regression model to investigate potential risk factors for HIV-negative status in individuals who report taking ART (e.g. sex, age group etc)
My intention is to generate a new variable (hivneg_art) with 1=HIV negative and taking ART, and 0=HIV positive and taking ART.
My code for analysis would then be:
xi:svy: logistic hivneg_art i.sex
Is this correct?
Thanks so much for the assistance.
|
|
|
|
Members
Search
Help
Register
Login
Home
