* Supplementary files: Stata codes used in our statistical analysis.


use "rawdata.dta", clear

*Q9. Please select the second to last option from the following choices.  (n=694) 

gen Q9_correct=. 
replace Q9_correct=1 if Q9==4  
replace Q9_correct=0 if Q9==1 | Q9==2 | Q9==3 | Q9==5 

 
*Q89: Respondents who used all drugs （those who selected "almost every day" OR "some days I use sometimes") (n=15) 
gen Q89_allari=0 
replace Q89_allari=1 if /// 
(Q891==4 | Q891==5) & /// 
(Q892==4 | Q892==5) & /// 
(Q893==4 | Q893==5) & /// 
(Q894==4 | Q894==5) & /// 
(Q895==4 | Q895==5) & /// 
(Q896==4 | Q896==5) & /// 
(Q897==4 | Q897==5) & /// 
(Q898==4 | Q898==5) & /// 
(Q899==4 | Q899==5) & /// 
(Q8910==4 | Q8910==5) 


*Q91 Those who selected one of the options 2, 3, 4, 5, 6, 7, for all disease names.(n=18) 
gen Q91_allaru=0 
replace Q91_allaru=1 if /// 
(Q911==2 | Q911==3 | Q911==4 | Q911==5 | Q911==6 | Q911==7) & /// 
(Q912==2 | Q912==3 | Q912==4 | Q912==5 | Q912==6 | Q912==7) & /// 
(Q913==2 | Q913==3 | Q913==4 | Q913==5 | Q913==6 | Q913==7) & /// 
(Q914==2 | Q914==3 | Q914==4 | Q914==5 | Q914==6 | Q914==7) & /// 
(Q915==2 | Q915==3 | Q915==4 | Q915==5 | Q915==6 | Q915==7) & /// 
(Q916==2 | Q916==3 | Q916==4 | Q916==5 | Q916==6 | Q916==7) & /// 
(Q917==2 | Q917==3 | Q917==4 | Q917==5 | Q917==6 | Q917==7) & /// 
(Q918==2 | Q918==3 | Q918==4 | Q918==5 | Q918==6 | Q918==7) & /// 
(Q919==2 | Q919==3 | Q919==4 | Q919==5 | Q919==6 | Q919==7) & /// 
(Q9110==2 | Q9110==3 | Q9110==4 | Q9110==5 | Q9110==6 | Q9110==7) & /// 
(Q9111==2 | Q9111==3 | Q9111==4 | Q9111==5 | Q9111==6 | Q9111==7) & /// 
(Q9112==2 | Q9112==3 | Q9112==4 | Q9112==5 | Q9112==6 | Q9112==7) & /// 
(Q9113==2 | Q9113==3 | Q9113==4 | Q9113==5 | Q9113==6 | Q9113==7) & /// 
(Q9114==2 | Q9114==3 | Q9114==4 | Q9114==5 | Q9114==6 | Q9114==7) & /// 
(Q9115==2 | Q9115==3 | Q9115==4 | Q9115==5 | Q9115==6 | Q9115==7) & /// 
(Q9116==2 | Q9116==3 | Q9116==4 | Q9116==5 | Q9116==6 | Q9116==7) & /// 
(Q9117==2 | Q9117==3 | Q9117==4 | Q9117==5 | Q9117==6 | Q9117==7) & /// 
(Q9118==2 | Q9118==3 | Q9118==4 | Q9118==5 | Q9118==6 | Q9118==7) 

 
*Totalmujun3　they should be excluded from analysis(n=720) 
gen totalmujun3=0 
replace totalmujun3=1 if /// 
Q9_correct==0 | /// 
Q89_allari==1 | /// 
Q91_allaru==1

* pregnancy or postpartum status
gen ninpu=.
replace ninpu=1 if ASC11==1
replace ninpu=1 if ASC12==1
replace ninpu=0 if ASC13==1
replace ninpu=0 if ASC14==1

replace ninpu=1 if BSC1S1==1 | BSC1S1==2 | BSC1S1==3 | BSC1S1==4 | BSC1S1==5 |BSC1S1==6
replace ninpu=0 if BSC1S2==1 | BSC1S2==2

label define ninpu 1"pregnant" 0"postpartum"
label values ninpu ninpu


* age
recode AGE (18/24=1)(25/29=2)(30/34=3)(35/48=4), gen(age4)
label define age4 1"18-24" 2"25-29" 3"30-34" 4"35-48"
label values age4 age4

* education(0: junior high school) Q7 Q8
gen edu=0 if Q7==1
replace edu=0 if Q7==2 & Q8==2
replace edu=0 if Q7==3 & Q8==2

// education=1:high school
replace edu=1 if ///
   Q7==2 & Q8==1  ///	
| (Q7==2 & Q8==3) ///  
| (Q7==3 & Q8==1) ///	
| (Q7==3 & Q8==3) ///	
| (Q7==4 & Q8==2) ///	
| (Q7==5 & Q8==2) ///	
| (Q7==6 & Q8==2) ///	
| (Q7==7 & Q8==2) ///	
| (Q7==8 & Q8==2)     
	
//education=2: college
replace edu=2 if ///
   Q7==4 & Q8==1 /// 
| (Q7==4 & Q8==3) /// 
| (Q7==5 & Q8==1) /// 
| (Q7==5 & Q8==3)  
	
//education==3: University
replace edu=3 if ///
   Q7==6 & Q8==1 ///
| (Q7==6 & Q8==3) /// 
| (Q7==7 & Q8==1) ///
| (Q7==7 & Q8==3) /// 
| (Q7==8 & Q8==1) /// 
| (Q7==8 & Q8==3) /// 
| (Q7==9 & Q8==2)	  
	
//education==4: Postgraduate
replace edu=4 if ///
   Q7==9 & Q8==1 /// 
| (Q7==9 & Q8==3)   
	

//Other
replace edu=0 if Q7==10

label define edu 0"no" 1"high school" 2"college" 3"university" 4"graduate school", replace
label values edu edu
		
//Education level
gen edu2=.
replace edu2=0 if edu>=0 & edu<2
replace edu2=1 if edu>=2 & edu<=4

recode edu2 (0=1)(1=2)
label define edu2 1"high school" 2"college"
label values edu2 edu2

* marital status Q5
recode Q5 (1=3) (2=1) (3/4=2) (5/6=3), gen (marital)
label define marital 1"married_cohabit" 2"married_separate" 3"single/divorced/widowed", replace
label values marital marital

* Income
gen hhincome=Q641 + Q642
replace hhincome=. if Q641==8888 | Q641==9999
replace hhincome=. if Q642==8888 | Q642==9999
replace hhincome=. if Q641==88888

centile hhincome, centile(25, 50, 75)  /*500, 670, 850*/

gen hhincome4=.
replace hhincome4=1 if hhincome>=0 & hhincome<500
replace hhincome4=2 if hhincome>=500 & hhincome<670
replace hhincome4=3 if hhincome>=670 & hhincome<850
replace hhincome4=4 if hhincome>=850 & hhincome<.
replace hhincome4=5 if hhincome==.
replace hhincome4=5 if Q641==8888 | Q641==9999 | Q641==88888 
replace hhincome4=5 if Q642==8888 | Q642==9999

label define hhincome4 1"Q1" 2"Q2" 3"Q3" 4"Q4" 5"declined", replace
label values hhincome4 hhincome4


* Parity: number of live birth Q202
recode Q202 (0=0)(1=1)(2/7=2)(39=.), gen(parity)
label define parity 0"0" 1"1" 2"2+", replace
label values parity parity

gen parity5=.
replace parity5=0 if ninpu==1 & parity==0
replace parity5=1 if ninpu==0 & parity==1
replace parity5=2 if ninpu==1 & parity==1
replace parity5=3 if ninpu==0 & parity==2
replace parity5=4 if ninpu==1 & parity==2

*Health literacy
gen hl=(Q801+Q802+Q803+Q804+Q805)/5
centile hl
gen hl2=1 if hl>=1 & hl<3.6
replace hl2=2 if hl>=3.6 & hl<=5
label define hl2 1"low" 2"high"
label values hl2 hl2

* Partner support
recode Q486 (1/2=1)(3/4=2), gen(pasupport)
replace pasupport=2 if marital==3

label define pasupport 1"supportive" 2"not supportive"
label values pasupport pasupport

* Current depression or mental illness
gen curdeppsy=0
replace curdeppsy=1 if Q9110==3 | Q9110==5  | Q9110==7 
replace curdeppsy=1 if Q9111==3 | Q9111==5  | Q9111==7

label define curdeppsy 1"yes" 0"no"
label values curdeppsy curdeppsy
  
* Current disability
gen disability=0
replace disability=1 if Q9112==3 | Q9112==5  | Q9112==7 
replace disability=1 if Q9113==3 | Q9113==5  | Q9113==7 
replace disability=1 if Q9114==3 | Q9114==5  | Q9114==7 
replace disability=1 if Q9115==3 | Q9115==5  | Q9115==7 
replace disability=1 if Q9116==3 | Q9116==5  | Q9116==7 

label define disability 1"yes" 0"no"
label values disability disability

*Current concerns at the time of survey
label variable Q511 "Q51.1 pregnancy delivery"
label variable Q512 "Q51.2 economy"
label variable Q513 "Q51.3 self"
label variable Q514 "Q51.4 partner"
label variable Q515 "Q51.5 family"
label variable Q516 "Q51.6 childcare"
label variable Q517 "Q51.7 work"
label variable Q518 "Q51.8 other"

*Study outcome: wanted to consult but could not
gen yakusyo=0
replace yakusyo=1 if Q521==3 | Q522==3

gen byoin=0
replace byoin=1 if Q526==3

gen Q531y=0
replace Q531y=1 if Q5311==1 | Q5321==1
replace Q531y=. if yakusyo==0

gen Q532y=0
replace Q532y=1 if Q5312==1 | Q5322==1
replace Q532y=. if yakusyo==0

gen Q533y=0
replace Q533y=1 if Q5313==1 | Q5323==1
replace Q533y=. if yakusyo==0

gen Q534y=0
replace Q534y=1 if Q5314==1 | Q5324==1
replace Q534y=. if yakusyo==0

gen Q535y=0
replace Q535y=1 if Q5315==1 | Q5325==1
replace Q535y=. if yakusyo==0

gen Q536y=0
replace Q536y=1 if Q5316==1 | Q5326==1
replace Q536y=. if yakusyo==0

gen Q537y=0
replace Q537y=1 if Q5317==1 | Q5327==1
replace Q537y=. if yakusyo==0

gen Q538y=0
replace Q538y=1 if Q5318==1 | Q5328==1
replace Q538y=. if yakusyo==0

gen Q539y=0
replace Q539y=1 if Q5319==1 | Q5329==1
replace Q539y=. if yakusyo==0

gen Q5310y=0
replace Q5310y=1 if Q53110==1 | Q53210==1
replace Q5310y=. if yakusyo==0

gen Q5311y=0
replace Q5311y=1 if Q53111==1 | Q53211==1
replace Q5311y=. if yakusyo==0

gen Q5312y=0
replace Q5312y=1 if Q53112==1 | Q53212==1
replace Q5312y=. if yakusyo==0

gen Q5313y=0
replace Q5313y=1 if Q53113==1 | Q53213==1
replace Q5313y=. if yakusyo==0

gen Q531b=0
replace Q531b=1 if Q5361==1
replace Q531b=. if byoin==0

gen Q532b=0
replace Q532b=1 if Q5362==1
replace Q532b=. if byoin==0

gen Q533b=0
replace Q533b=1 if Q5363==1
replace Q533b=. if byoin==0

gen Q534b=0
replace Q534b=1 if Q5364==1
replace Q534b=. if byoin==0

gen Q535b=0
replace Q535b=1 if Q5365==1
replace Q535b=. if byoin==0

gen Q536b=0
replace Q536b=1 if Q5366==1
replace Q536b=. if byoin==0

gen Q537b=0
replace Q537b=1 if Q5367==1
replace Q537b=. if byoin==0

gen Q538b=0
replace Q538b=1 if Q5368==1
replace Q538b=. if byoin==0

gen Q539b=0
replace Q539b=1 if Q5369==1
replace Q539b=. if byoin==0

gen Q5310b=0
replace Q5310b=1 if Q53610==1
replace Q5310b=. if byoin==0

gen Q5311b=0
replace Q5311b=1 if Q53611==1
replace Q5311b=. if byoin==0

gen Q5312b=0
replace Q5312b=1 if Q53612==1
replace Q5312b=. if byoin==0

gen Q5313b=0
replace Q5313b=1 if Q53613==1
replace Q5313b=. if byoin==0

** Analysis
drop if totalmujun3==1
drop if parity==.

* Table 1　basic characteristics
ta age4
ta edu2
ta marital
ta hhincome4
ta parity
ta hl2
ta pasupport
ta curdeppsy
ta disability

ta Q511 
ta Q512 
ta Q513 
ta Q514 
ta Q515 
ta Q516 
ta Q517 
ta Q518

ta yakusyo
ta byoin


*Figure 1 the reasons that she wanted to consult with public service but could not
ta Q531y 
ta Q532y 
ta Q533y 
ta Q534y 
ta Q535y
ta Q536y 
ta Q537y 
ta Q538y 
ta Q539y 
ta Q5310y
ta Q5311y 
ta Q5312y 
ta Q5313y


**Figure2 the reasons that she wanted to consult with hospital but could not
ta Q531b 
ta Q532b 
ta Q533b 
ta Q534b 
ta Q535b 
ta Q536b 
ta Q537b 
ta Q538b 
ta Q539b 
ta Q5310b 
ta Q5311b 
ta Q5312b 
ta Q5313b


*Table 3 Associated factors: wanted to consult with public service
logistic yakusyo b4.i.age4 b2.i.edu2 b4.i.hhincome4 b0.i.parity5 b2.hl2 b1.pasupport b0.curdeppsy b0.disability

* Table 4 Associated factors: wanted to consult with hospital
logistic byoin b4.i.age4 b2.i.edu2 b4.i.hhincome4 b0.i.parity5 b2.hl2 b1.pasupport b0.curdeppsy b0.disability