Immunization coverage, knowledge, satisfaction, and associated factors of non-National Immunization Program vaccines among migrant and left-behind families in China: evidence from Zhejiang and Henan provinces

Abstract

Background

Migrant and left-behind families are vulnerable in health services utilization, but little is known about their disparities in immunization of non-National Immunization Program (NIP) vaccines. This study aims to evaluate the immunization coverage, knowledge, satisfaction, and associated factors of non-NIP vaccines among local and migrant families in the urban areas and non-left-behind and left-behind families in the rural areas of China.

Methods

A cross-sectional survey was conducted in urban areas of Zhejiang and rural areas of Henan in China. A total of 1648 caregivers of children aged 1–6 years were interviewed face-to-face by a pre-designed online questionnaire, and their families were grouped into four types: local urban, migrant, non-left-behind, and left-behind. Non-NIP vaccines included Hemophilus influenza b (Hib) vaccine, varicella vaccine, rotavirus vaccine, enterovirus 71 vaccine (EV71) and 13-valent pneumonia vaccine (PCV13). Log-binomial regression models were used to calculate prevalence ratios (PRs) and 95% confidence intervals (CIs) for the difference on immunization coverage of children, and knowledge and satisfaction of caregivers among families. The network models were conducted to explore the interplay of immunization coverage, knowledge, and satisfaction. Logistic regression models with odds ratios (ORs) and 95% CIs were used to estimate the associated factors of non-NIP vaccination.

Results

The immunization coverage of all non-NIP vaccines and knowledge of all items of local urban families was the highest, followed by migrant, non-left-behind and left-behind families. Compared with local urban children, the PRs (95% CIs) for getting all vaccinated were 0.65 (0.52–0.81), 0.29 (0.22–0.37) and 0.14 (0.09–0.21) among migrant children, non-left-behind children and left-behind children, respectively. The coverage-knowledge-satisfaction network model showed the core node was the satisfaction of vaccination schedule. Non-NIP vaccination was associated with characteristics of both children and caregivers, including age of children (> 2 years-OR: 1.69, 95% CI: 1.07–2.68 for local urban children; 2.67, 1.39–5.13 for migrant children; 3.09, 1.23–7.76 for non-left-behind children); and below caregivers’ characteristics: family role (parents: 0.37, 0.14–0.99 for non-left-behind children), age (≤ 35 years: 7.27, 1.39–37.94 for non-left-behind children), sex (female: 0.49, 0.30–0.81 for local urban children; 0.31, 0.15–0.62 for non-left-behind children), physical health (more than average: 1.58, 1.07–2.35 for local urban children) and non-NIP vaccines knowledge (good: 0.45, 0.30–0.68 for local urban children; 7.54, 2.64–21.50 for left-behind children).

Conclusions

There were immunization disparities in non-NIP vaccines among migrant and left-behind families compared with their local counterparts. Non-NIP vaccination promotion strategies, including education on caregivers, and optimization of the immunization information system, should be delivered particularly among left-behind and migrant families.

Graphical Abstract

Background

With the development of economic and urbanization, a large proportion of people move from rural to urban areas in order to get a better life in the context of national labor shortages in China over the past few decades [1, 2]. The mass migration of laborers caused a growth of two special groups: left-behind families in rural areas and migrant families in urban areas [3,4,5]. Left-behind families refer to one or both parents migrated into cities for work, leaving their children in the rural communities with other caregivers (e.g., grandparents) for over six months [6]; while migrant families refer to parents migrated into cities for work together with their children for over six months [7]. Till 2022, there was an estimated 6.44 million left-behind children and over 70 million migrant children in China [8]. According to the household registration policy in China, “hukou” gives households access to social benefits in their registration areas but limits access to those outside their registration area, which causes that migrant families often have less access to social resources, including housing, education and medical services [9]. On the other hand, left-behind families in rural areas, characterized by parent–child separation [10], were also socially disadvantaged. Therefore, compared to children of local families, those of migrant and left-behind families usually experienced poorer health, lower quality of life, and lower utilization of health services [11,12,13]. The widening health inequality in health services use (e.g., immunization) in migrant and left-behind children has become an important issue in China [9].

Immunization has been considered as one of the most successful strategies for prevention and control of infectious diseases. Since the National Immunization Program (NIP) was implemented in 1978, China has been working to promote childhood immunization and reduce national burden of vaccine-preventable diseases among children [14]. In China, vaccines were divided into two categories: NIP and non-NIP vaccines [15]. Non-NIP vaccines, an alternative and supplement to NIP vaccines, also contribute to reduce morbidity of some essential infectious diseases (e.g., hepatitis A, and influenza) [16]. Unlike NIP vaccines to be free and mandatory, non-NIP vaccines are voluntarily administered and self-funded, therefore families usually show inadequate knowledge level and low acceptance of non-NIP vaccination for their children, which even widens disparities of immunization among migrant and left-behind families [15, 17].

Early findings have assessed the disparities in immunization between migrant or left-behind children and local children. A systematic review and meta-analysis of 11 studies found rural–urban migrant children in China, India and Nigeria had lower immunization coverage of NIP vaccines than local urban children [2]. A cross-sectional study conducted in Guangdong Province of China also found the immunization coverage of NIP vaccines among migrant children aged 12–59 months was quite low, influenced by the sex and birth place of children, and the caregivers’ occupation, knowledge and attitude towards immunization, as well as the family income [18]. For left-behind children, one study suggested their age-appropriate vaccination coverage for NIP vaccines was significantly lower than that in non-left-behind children [19]. Moreover, the association of immunization coverage, knowledge, practice, and experiences has been reported in several studies. For example, a study among 1820 migrant children aged 12–35 months in Beijing of China suggested that the age-appropriate immunization coverage was significantly influenced by the caregivers’ immunization knowledge but not the satisfaction with vaccination services [20]. Findings from guardians of 0–59-month-old Chinese children reported the coverage of influenza vaccines was negatively related to hesitancy, while vaccine hesitancy could decreased with knowledge and vaccination experience improved [21]. However, previous studies have several limitations: (1) mostly focused on NIP vaccines, (2) only focused on one or a few types of non-NIP vaccines or did not specify the category of vaccines, and (3) did not include children from both the left-behind and migrant families and their local counterparts in one study design. To narrow or eliminate the immunization disparities among migrant and left-behind children, it is necessary to evaluate the immunization coverage, knowledge and satisfaction and their interplays among migrant and left-behind families, particularly in one study design.

Using information from migrant and left-behind families and their local counterparts in China, this study was conducted (1) to estimate the immunization coverage of non-NIP vaccines; (2) to estimate the immunization knowledge and satisfaction of non-NIP vaccines, and their association with non-NIP coverage; and (3) to explore the influence of characteristics of the children and caregivers on non-NIP vaccines coverage.

Methods

Study design and data collection

A cross-sectional survey was conducted during July and October 2022 in Zhejiang and Henan Provinces respectively, either of which is the major labor-importing or labor-exporting province in China and exists a large number of migrant families or left-behind families. The target population of this survey was made up of caregivers of children aged 1–6 years, who resided in urban areas of Xiaoshan District, Hangzhou City, Zhejiang Province and rural areas of Song County, Luoyang City, Henan Province. Families of caregivers included in Xiaoshan district were migrant families and local urban families, and those of caregivers included in Song County were left-behind families and non-left-behind families. Caregivers were excluded if (1) they have been residing locally for < 6 months; or (2) their children could not be vaccinated because of severe illness or disabilities.

Five towns/communities were selected in Song County and Xiaoshan District using the simple random sampling method, and caregivers were recruited in local township health centers and community health service centers. By reviewing the immunization coverage of non-NIP vaccines among urban and rural children (75.8%), using a desired precision of ± 2% with 95% confidence intervals, and assuming a non-response rate of 10.0%, the required number of surveyed caregivers of children was 1937.

The questionnaire included five sections (characteristics of children, characteristics of caregivers, immunization coverage of children, immunization knowledge of caregivers, and immunization satisfaction of caregivers) and 51 questions. The questionnaire was administered on the online platform Wenjuanxing (https://www.wjx.cn/), the most popular web-based questionnaire platform in China. There were standard instructions for the questionnaire survey, and a total of 35 primary care physicians were invited to administrate the questionnaires. It would take about 10 min to complete the questionnaire, and the caregivers were incentivized to participate in the survey by providing reasonable money. In our survey, only one caregiver per child was face-to-face interviewed, whose answers were entered real-time to the Wenjuanxing platform by primary care physicians. Before the survey, each caregiver of the child enrolled was asked to provide electronic informed consent on the Wenjuanxing online platform. Zhejiang University School of Public Health Medicine Ethics Committees approved the study protocol (ZGL202206-6).

Immunization coverage, knowledge, satisfaction

The five types of non-NIP vaccines included in our survey were Haemophiles influenza b (Hib) vaccine, varicella vaccine, rotavirus vaccine, enterovirus 71 vaccine (EV71) and 13-valent pneumonia vaccine (PCV13). These five vaccines were commonly used in China, but have not been included in the NIP [22]. Caregivers were asked whether their children had been vaccinated the aforementioned non-NIP vaccines, and the coverage of non-NIP vaccines was calculated among local urban, migrant, non-left-behind and left-behind children respectively.

The immunization knowledge of caregivers was measured using seven items, including convenience, category, efficiency, continuity, time, schedule, and adverse events of non-NIP vaccines, which were calculated among local urban, migrant, non-left-behind and left-behind families respectively. The description of each item was provided to caregivers, and they were asked whether they aware of these seven items. For example, the knowledge item of “convenience” was measured by asking caregivers whether they were aware of that children could accept non-NIP vaccination nationwide if they had a vaccination certificate.

The immunization satisfaction of caregivers was investigated using nine items, including convenience, vaccination reminder, vaccination environment, consultation, vaccination skills, service quality, vaccination process, vaccination education and time of vaccination, which were calculated among local urban, migrant, non-left-behind and left-behind families respectively. Similarly, a description of each item was provided to caregivers, and caregivers were asked whether they were satisfied with the aforementioned nine items. For example, the experience item of “schedule” was measured by asking caregivers whether they were satisfied with the non-NIP vaccination schedule of local clinics.

Questions for immunization knowledge and satisfaction were general to all non-NIP vaccines, and the whole questionnaire was presented in the Additional file 1: Questionnaire.

Characteristics of families

The characteristics of children included age, sex (male; female) and birth order (first-born; later-born). The characteristics of caregivers included family role (parents; others), age (≤ 35 years; 35–40 years; ≥ 40 years), sex (male; female), education level (elementary school or lower; middle school; junior college or higher), total household income (less than average; more than average), physical health [assessed by 12-Item Short Form Survey (SF-12)] and mental health (assessed by SF-12). Total household income was measured by the sum of earning income, capital income, pension income, income from government transfers, other income and the total income from other household members during last year.

Statistical analysis

Characteristics of families were described as mean [standard deviation (SD)] for continuous variables and as number (percentage) for categorical variables. The differences of variables across groups by family types were compared using one-way analysis of variance (ANOVA) and Chi-square test.

Log-binomial regression models were used to calculate prevalence ratios (PRs) and 95% confidence intervals (CIs) for the difference on immunization coverage of children, and knowledge and satisfaction of caregivers of different family types, regarding local urban families as the reference. Age and sex of children were adjusted in the log-binomial regression models. The coverage-knowledge-satisfaction network models of non-NIP vaccination were conducted through Mixed Graphical Model (MGM) using the R package “qgraph” [23]. Least Absolute Shrinkage and Selection Operator (LASSO) and Extended Bayesian Information Criterion (EBIC) were utilized to reduce the pseudo-correlation of connections in the network and to find the optimal fitting model [24]. The nodes of the network represent different items in coverage, knowledge, and satisfaction of non-NIP vaccines, while the bridges represent the correlation between nodes, which are represented by calculating the coefficient of partial correlation among the items of vaccines coverage, vaccination knowledge and satisfaction. The closer the nodes are and the thicker the connected line segments are, the stronger the correlations between nodes are. Strength, defined as the absolute value of the shortest distance between all connections of a node, was chosen to represent the centrality of each node, and its magnitude could represent the extent of connection between nodes in the network. Bridge expected influence refers to the sum of the absolute value of the shortest distance between a specific node and all other connected nodes, which was proportionable to the degree of network conduction through the node [25]. In addition, we used stability coefficient (CS) and 95% CI of edge weights to evaluate the stability of the network. CS coefficient represents the maximum sample attenuation ratio when the correlation value between the original network and the regenerative network parameters is 0.70 [26].

The multivariable logistic regression models to conducted to explore the associated factors of non-NIP vaccination among four types of families, with the outcomes of non-NIP vaccination divided into two groups: “having received all five non-NIP vaccines” and “having not received at least one non-NIP vaccines”. Odds ratios (ORs) and 95% CIs were calculated to compare the association of the characteristics of children and caregivers with non-NIP vaccination for each type of families.

Analyses were performed using SAS (Version 9.4, SAS Institute Inc., Cary, NC) and R (Version 3.6.1; R studio, Boston, Massachusetts). The study was conducted and reported in line with the Consensus-Based Checklist for Reporting of Survey Studies (CROSS, Additional file 2: Table).

Results

Basic characteristics

A total of 2366 caregivers were invited to participate in the survey, 2031 agreed to start, and 2018 completed the survey, resulting in response rate of 85.3% and completion rate of 99.4%. Considering the first shot of included non-NIP vaccines were scheduled before one year old, we further excluded respondents with children less than one year old to avoid underestimating the immunization of coverage non-NIP vaccines. Therefore, a number of 1648 caregivers of children aged 1–6 years were included in the analyses. The average age of children was 3.45 (± 1.59) years, and 47.9% of the children were females. The characteristics of four family types are presented in Table 1. The number of local urban children, migrant children, non-left-behind children and left-behind children was 517 (31.4%), 276 (16.8%), 488 (29.6%) and 367 (22.3%), respectively. Migrant children are more likely to be younger and first-born, and their caregivers are more likely to be parents, 35 years or younger and females, and have lower household income. Left-behind children are more likely to be older and first-born, and their caregivers are more likely to be other than parents, ≥ 40 years old and males, and have lower education level and higher household income.

Table 1 Basic characteristics of children and the caregivers according to four types of family in urban areas in Zhejiang Province and rural areas in Henan Province, China 2022 (n = 1648)

Immunization coverage, knowledge, and satisfaction

Immunization coverage

The immunization coverage of non-NIP vaccines among children of different family types is presented in Fig. 1. The coverage for single non-NIP vaccine was highest among local urban children, ranging from 61.1% for EV71 to 85.7% for varicella vaccine. The coverage for varicella vaccine, rotavirus vaccines, EV71 and PCV13 were lowest among left-behind children, which were 55.9%, 16.9%, 41.7%, and 8.2%, respectively. The coverage of all five non-NIP vaccines was highest among local urban children (41.2%), and was lowest among left-behind children (5.7%). After adjustment for age and sex of children, the PRs (95% CIs) for getting all vaccinated were 0.65 (0.52–0.81), 0.29 (0.22–0.37) and 0.14 (0.09–0.21) among migrant children, non-left-behind children and left-behind children, respectively.

Fig. 1

The immunization coverage, prevalence ratio (PRs) and 95% confidence intervals (CIs) of non-National Immunization Program vaccines among children from four types of family in urban areas in Zhejiang Province and rural areas in Henan Province, China 2022. Hib Hemophilus influenza b; EV71 enterovirus 71 vaccine, PCV13 13-valent pneumonia vaccine; PR prevalence ratio

Immunization knowledge

Table 2 shows the immunization knowledge of non-NIP vaccines among the caregivers of different family types. Among the seven knowledge items of non-NIP vaccination, vaccination convenience, category, time and schedule were most aware by the caregivers, while the knowledge of vaccination continuity was lowest. For example, the awareness rate of vaccination time and schedule was 99.2% and 98.9% respectively, while that of vaccination continuity was 27.3% among caregivers from left-behind families. The rate of being aware of all knowledge items was highest in local urban families (41.2%), followed by migrant (23.2%), non-left-behind (19.5%) and left-behind families (11.4%). Compared to caregivers from local urban families, the PRs (95% CIs) for being aware of all knowledge items were 0.55 (0.43–0.70), 0.47 (0.38–0.58) and 0.28 (0.21–0.38) among those from migrant, non-left-behind and left-behind families, respectively.

Table 2 Awareness rate, prevalence ratios (PRs) and 95% confidence intervals (CIs) of immunization knowledge of non-NIP vaccines among the caregivers from four family types in urban areas in Zhejiang Province and rural areas in Henan Province, China 2022 (n = 1648)

Immunization satisfaction

Table 3 shows the immunization satisfaction of non-NIP vaccines among the caregivers of different family types. More than 90% of the caregivers reported satisfied with the nine items of non-NIP vaccination satisfaction. The satisfaction rate was 90.5% (vaccination reminder)–94.6% (vaccination skills), 91.3% (vaccination education)–96.0% (vaccination skills), 90.4% (convenience)–96.1% (vaccination skills), and 95.1% (convenience)–99.2% (vaccination reminder) among those from local urban, migrant, non-left-behind families and left-behind families. Compared to caregivers from local urban families, caregivers from left-behind families had higher satisfaction rate of non-NIP vaccination (PR: 1.12, 1.05–1.19 for being satisfied with all items), while the PRs for being satisfied with all items were not significant among those from migrant and non-left-behind families.

Table 3 Satisfaction rate, prevalence ratios (PRs) and 95% confidence intervals (CIs) of non-NIP vaccination among the caregivers from four family types in urban areas in Zhejiang Province and rural areas in Henan Province, China 2022 (n = 1648)

The interrelationship among immunization coverage, knowledge and satisfaction

Figure 2 displays the coverage-knowledge-satisfaction network model of non-NIP vaccination for children aged 1–6 years. In the types of vaccines, the connection between rotavirus vaccine and PCV13 was strongest, with the edge weight of 0.47. In the knowledge items of non-NIP vaccination, the connection between convenience and vaccines category, as well as the connection between vaccination continuity and adverse events were strongest, with the edge weights of 0.20. In the immunization satisfaction items of non-NIP vaccines, the connection between convenience and vaccination environment was strongest (edge weight: 0.36), followed by the connection between vaccination schedule and vaccination education (edge weight: 0.35). The core nodes in the network models includes vaccination schedule, rotavirus vaccine, and vaccination education due to their high strength. The primary bridge between non-NIP vaccines coverage and knowledge was the rotavirus vaccine node, and the second ones were vaccination continuity and vaccine category.

Fig. 2

Coverage-knowledge-satisfaction network model of non-NIP vaccination among included families in urban areas in Zhejiang Province and rural areas in Henan Province, China 2022. Hib Hemophilus influenza b; EV71 enterovirus 71 vaccine, PCV13 13-valent pneumonia vaccine

Associated factors of non-NIP vaccination

Table 4 shows the association of characteristics of families with immunization coverage of non-NIP vaccines across four family types. Children aged > 2 years among local urban (OR = 1.69, 95% CI: 1.07–2.68), migrant (2.67, 1.39–5.13) and non-left-behind families (3.09, 1.23–7.76) were more likely to get all vaccines. Better physical health of caregivers from local urban families was associated with an increased odds (1.58, 1.07–2.35) of receiving all non-NIP vaccines for their children, while female caregivers (0.49, 0.30–0.81), and caregivers with better immunization knowledge (0.45, 0.30–0.68) showed lower probability to have their children vaccinated all non-NIP vaccines. Among non-left-behind families, younger caregivers were more likely to get their children vaccinated (7.27, 1.39–37.94), while those being parents (0.37, 0.14–0.99) or females (0.31, 0.15–0.62) were less likely to do so. As for left-behind families, good immunization knowledge of caregivers increased their children's odds of have all vaccinated (7.54, 2.64–21.50).

Table 4 Odds ratios (ORs) and 95% confidence interval (CIs) for the association of characteristics of families with non-NIP vaccine coverage among four types of family in urban areas in Zhejiang Province and rural areas in Henan Province, China 2022 (n = 1648)

Discussion

Through this cross-sectional study in Zhejiang and Henan Provinces, we found that the immunization coverage and knowledge of non-NIP vaccines among local urban children and caregivers was highest, followed by migrant, non-left-behind and left-behind children and caregivers. The satisfaction rate of non-NIP vaccination was quite high (> 90%) among the caregivers of four family types, and those from left-behind families were more satisfied with non-NIP vaccination compared to local urban children. The immunization coverage of non-NIP vaccines was associated with characteristics of children and caregivers, including age of children, and family role, age, sex, physical health and non-NIP vaccination knowledge of caregivers. The core node of the interrelationship among non-NIP vaccination coverage, knowledge and satisfaction was the satisfaction of vaccination schedule.

Prior research has provided information about the disparities on immunization coverage between rural and urban children [27, 28], while our study focused on the inequality in non-NIP vaccination among two vulnerable groups of migrant and left-behind children and their counterparts. Our findings were supported by a previous qualitative study of immunization providers in Sichuan, Guangdong and Henan provinces in China, which revealed the low vaccination coverage in migrant and left-behind children [29]. Compared to the caregivers from local urban families, we found those from migrant, non-left-behind and left-behind families showed significantly lower knowledge on non-NIP vaccination. Caregivers from migrant families and rural families may have lower socio-economic status, relatively scarce sources of information [30, 31], and thus have poorer understanding on non-NIP vaccination and less access to healthcare utilization. To improve vaccination service quality and reduce rural-urban health inequalities, children from migrant and left-behind families should be prioritized in the non-NIP vaccines promotion programs. Vaccination education could also have the potential to improve the caregivers’ knowledge and acceptance on non-NIP vaccines, strengthen their connection with local vaccinators, and treat non-NIP vaccination correctly. However, the caregivers from left-behind families in our study showed highest satisfaction on non-NIP vaccination. The caregivers from left-behind families were mainly grandparents, who lived in rural communities for decades, and tended to have a harmonious relationship with community vaccinators [32]. Results from the coverage-knowledge-satisfaction network models of non-NIP vaccination suggesting optimization on vaccination schedule would contribute to caregivers’ satisfaction on non-NIP vaccination for their children, and thus improve immunization coverage of non-NIP vaccines.

The association between socio-demographic characteristics of families and immunization coverage of non-NIP vaccines have been suggested in previous observational studies. An online cross-sectional survey conducted in Jiangsu Province of China suggested parents with younger age, lower education level and health-related occupations had lower acceptance on non-NIP vaccines for their children [33]. As suggested by a recent qualitative study, the child-related determinants of low immunization coverage included sex, siblings and health conditions, while the caregiver-related determinants included socio-economic status, family role, education level and ethnicity [29]. In our study, children aged > 2 years were more likely to receive all five non-NIP vaccines among most types of families. It might be because older children were physically stronger than younger ones, and had less contradictions for vaccination [34]. In addition, some caregivers might have poor access to the vaccination information, especially those from migrant and rural families, and non-NIP vaccinations would be given as the child grows older [34]. Caregivers with better physical health and younger ages would have higher autonomy and energy to get their children received non-NIP vaccines, while caregivers being parents or females, who tended to be more concerned about vaccine adverse events, showed lower probability to take their children for non-NIP vaccination [35]. These findings indicated non-NIP vaccination promotion should be delivered early for families with age-eligible children, and younger family members should be encouraged to participate in the decision-making process of non-NIP vaccination. The opposite association between immunization knowledge level and coverage was observed between local urban families and left-behind families, which might be explained by the difference in sources of information [36]. Caregivers from local urban families had more sources to receive information of vaccination, and the information could be mixed and biased, leading caregivers to express more concerns on non-NIP vaccination [37]. On the contrary, caregivers from left-behind families, mainly grandparents, had more satisfaction with local vaccination services as mentioned above, and therefore were more willing to take their children for non-NIP vaccination. Our findings add weight to the evidence that the effects of child-related factors (e.g., age) and family structure and socio-economic status on the caregiver’s decision-making process on non-NIP vaccination should be stressed. To inform broadly applicable strategies for promoting non-NIP vaccination among migrant and left-behind children, the government should add more vaccines to publicly funded immunization schedule or make policy to decrease the price of non-NIP vaccines to an affordable extent.

Strengths of our study include the consideration of families of four migration types in one study design and the use of network models. In addition, our study had no missing data of variables, since caregivers were incentivized to participate in the survey by providing reasonable money. Limitations to this study also warrant consideration. First, the cross-sectional design precluded us to capture the dynamic change of immunization coverage, knowledge and satisfaction, and to explore the reason of low immunization coverage in migrant and left-behind children. Second, the survey was conducted in local township health centers and community health service centers, which could cause selection bias. For example, the vulnerable groups, migrant families and left-behind families were less accessible to the local vaccine clinics, leading to the overestimate of non-NIP vaccines coverage in these population. Third, recall bias and social desirability bias might exist, and the latter referred when respondents give answers to questions that they believe will make them look good to others, concealing their true opinions or experiences [38]. For example, caregivers might report higher satisfaction with non-NIP vaccination. Fourth, our study only included five non-NIP vaccines (Hib vaccine, varicella vaccine, rotavirus vaccine, EV71 and PCV13), and other common non-NIP vaccines, such us influenza vaccine, were not considered. Fifth, due to the small sample size, the sample of some categories of variables were too small (e.g., nationality of children and caregivers, and the occupation of caregivers) to be included in the models. The small sample size also limited the further network analyses among families of four types. Future large studies are warranted to evaluate more associated factors of non-NIP vaccines coverage and the interplay of immunization coverage, knowledge and satisfaction of non-NIP vaccines among migrant and left-behind families. Sixth, we did not assess NIP vaccination as comparators, and whether the immunization disparities among migrant and left-behind families were specific to non-NIP vaccines or general to both NIP and non-NIP vaccines remained unclear. Finally, the sex of children and caregivers only included two categories of males and females, and “intersex” was not included. In addition, some important information of non-NIP vaccines, including the price and accessibility, was not surveyed. However, we investigated the effects of socioeconomic status of families on non-NIP vaccination, also providing support for lowering the prices of non-NIP vaccines in the future.

Conclusions

The disparities of immunization coverage and knowledge of non-NIP vaccines existed among migrant families in urban areas, and non-left-behind and left behind families in rural areas. Age of children, physical health of caregivers and family structures were associated with the immunization coverage of non-NIP vaccines. When attempting to reduce health inequality in non-NIP vaccination for the children, attention should be paid not only the rural–urban disparity but also the vulnerable groups of left-behind and migrant children.