variables for svyset in Stata (Bangladesh 2011) [message #3566] |
Fri, 09 January 2015 13:31 |
kmorris
Messages: 5 Registered: December 2014
|
Member |
|
|
Hello,
I am currently trying to analyze c-section by wealth quintile (lower 2 and upper 3 as their own groupings).
I am using Stata 12, and when trying to set svyset for the data, I have been encountering an issue with the strata that is affecting my results in analysis.
gen birth5=0
replace birth5=1 if v208>0
label var birth5 "have a live birth in past 5 years"
label define yesno 0 "no" 1 "yes"
label values birth5 yesno
keep if birth5 == 1
gen wt=v005/1000000
gen csect = . if m17 == .
replace csect = 1 if m17 == 1
replace csect = 0 if m17 == 0
label variable csect "C-Section"
label define csect 0 "no" 1 "yes"
label value csect csect
gen psu=v021
gen strata=v023
**note, I generated the weight earlier, above)**
svyset psu [pweight=wt], strata(strata)
**now analyzing using poisson for the lower 2 wealth quintiles**
svy: poisson csect v190 if v190<3
And this is where I run into errors, see results as written below
(running poisson on estimation sample)
Survey: Poisson regression
Number of strata = 19 Number of obs = 3645
Number of PSUs = 499 Population size = 3867.0893
Design df = 480
F( 0, 480) = .
Prob > F = .
------------------------------------------------------------ ------------------
| Linearized
csect | Coef. Std. Err. t P>|t| [95% Conf. Interval]
-------------+---------------------------------------------- ------------------
v190 | 1.011505 . . . . .
_cons | -4.721347 . . . . .
------------------------------------------------------------ ------------------
Note: missing standard errors because of stratum with single sampling unit.
Did I do something wrong with the sampling unit? Is this error unique to this dataset or will I have to change my syntax?
Thank you!
Kate Morris
|
|
|