Study sample
We used data from the Longitudinal Survey of Newborns in the 21st Century, a population-based survey conducted by the Ministry of Health, Labor, and Welfare in Japan. The study sample included all babies born in Japan between January 10 and 17 or July 10 and 17, 2001, using birth records from national vital statistics. The baseline survey was mailed to parents when their infants were six months old (n = 53,575). A total of 47,015 caregivers responded to the baseline questionnaire (response rate: 87.7%). After the baseline survey, annual surveys were conducted by sending questionnaires to participants by mail. We used data from 2001, 2002, and 2006 survey waves. A total of 43,925 (93.4%) caregivers responded in the 2002 survey wave, and 38,540 (82.0%) responded in the 2006 survey wave. We included children who lived with both parents and children whose mothers answered the questionnaire to maintain consistency in the assessment of childcare (excluded samples n = 3,458, 9.0%). We excluded responses with missing data for the variables used in the analysis: child problem behavior (n = 1,443, 3.7%), paternal involvement in childcare and domestic chores (n = 2,533, 6.6%), maternal involvement in childcare and domestic chores (n = 1,277, 3.3%), parental education (n = 578, 1.5%), parental employment (n = 807, 2.1%), annual household income in 2002 (n = 2,705, 7.0%), gestational age (n = 14, 0.04%), weight at birth (n = 7, 0.02%), and child temperament (n = 374, 1.0%). Finally, 27,870 newborns were included in the analyses (72.3% of the respondents in the 2006 survey wave).
Exposure: paternal childcare
Paternal childcare was assessed at 18 months of age using the following two measurements: paternal childcare hours on weekdays or weekends, and frequency of each type of caregiving (feeding, changing diapers, bathing, putting the child to sleep, playing with the child at home, and taking the child outside).
Paternal childcare hours were ascertained by asking the following questions: How much time does the father spend with his child on average in one day, except for sleeping hours? Respondents answered in respect of both weekdays and weekends by choosing from the following options: none, < 0.5, 0.5-0.9 hours, 1-1.9 hours, 2-3.9 hours, 4-5.9 hours, and ≥ 6 hours. Then, we defined the categories, making sure not to highly skew the distribution of each one, into “less than 0.5 hours,” “0.5 to 1.9 hours,” and “2 or more hours” for weekdays, and “less than 4 hours,” “4 to 5.9 hours,” and “6 or more hours” for weekends.
The frequency of each type of paternal caregiving was retrieved using six items: 1) feeding, 2) changing diapers, 3) bathing, 4) putting the child to sleep, 5) playing with the child at home, and 6) taking the child outside. Responses for each question included “not at all,” “rarely,” “sometimes,” and “always” For these categories, we classified “always” and “sometimes” as a high degree of paternal caregiving, and “rarely” and “not at all” as a low degree of paternal caregiving. To calculate the total paternal caregiving scores, each response was scored from 0 to 3 (i.e., “not at all” = 0 and “always” = 3). Thus, the total paternal caregiving scores ranged from 0 to 18 (Cronbach's alpha = 0.77). Paternal caregiving in one type of activity was not strongly correlated with another type of activity (Spearman’s correlation coefficients ranging from 0.27-0.50, Table 1). The measurement and the method of specifying the caregiving variable were in line with previous studies [9, 22]. We divided the total caregiving scores into three groups, namely the “high degree of paternal caregiving” group with more than 1 SD above the mean, the “low degree of paternal caregiving” group with less than 1 SD under the mean, and the “middle degree of paternal caregiving” group.
Table 1. Spearman correlation coefficients between specific types of paternal involvement in caregiving
|
1)
|
2)
|
3)
|
4)
|
5)
|
6)
|
1) Feeding
|
1.00
|
|
|
|
|
|
2) Changing diaper
|
0.50
|
1.00
|
|
|
|
|
3) Bathing
|
0.35
|
0.31
|
1.00
|
|
|
|
4) Putting the child to sleep
|
0.42
|
0.42
|
0.40
|
1.00
|
|
|
5) Playing with the child at home
|
0.35
|
0.31
|
0.36
|
0.30
|
1.00
|
|
6) Taking the child outside
|
0.32
|
0.30
|
0.27
|
0.28
|
0.44
|
1.00
|
Outcome: child behavioral problems
Child behavioral problems were assessed when the children were 5.5 years old by asking the following six questions which required a yes/no response: 1) Is your child able to listen without fidgeting? 2) Is your child able to focus on a specific task? 3) Is your child patient? 4) Is your child able to express their emotions appropriately? 5) Is your child able to behave in a group situation? 6) Is your child able to keep promises? These variables were developed to identify early signs of behavioral and developmental problems and have been previously used as a set of measurements for behavioral problems [23, 24]. In addition, we defined children with any of these problem behaviors as “children with behavioral problems.”
Covariates
We selected the following variables as potential confounders: the number of siblings, living with grandparent(s), paternal and maternal age at childbirth (< 25, 25-29, 30-34, 35-39, ≥ 40), paternal and maternal education (lower than high school degree, high school degree, 2-year college or vocational school degree, and a college degree or more), annual household income in 2002 (< JPY 4 million, 4-5.9 million, 6-7.9 million, 8-9.9 million, ≥ 10 million), maternal childcare hours on weekdays or weekends at 18 months old, and total scores of maternal involvement in caregiving. Maternal involvement in caregiving was scored using the same items used to assess fathers’ involvement in caregiving, and the summed score of the six items was used in the analysis.
How fathers are involved with their children is affected by factors relating to the children themselves, such as disabilities and temperament [25, 26]. Therefore, the baseline characteristics of children should also be considered when estimating the effect of paternal childcare: gestational age (< 37 weeks, 37-41 weeks, ≥ 42 weeks), child’s history of hospital admission or visits for congenital diseases, and child’s temperament (self-reported dichotomous answers to questions about 18 traits: active, shy of strangers, playful, short-tempered, careful, aggressive, timid, competitive, obedient, independent, stubborn, full of curiosity, fickle, restless, shy, spoiled, easygoing, impatient).
Statistical analysis
We developed a logistic regression model with odds ratios (ORs) and 95% confidence intervals (95% CIs) calculated by adjusting the covariates shown above. In addition to the crude models, two adjusted models were used: Model 1 adjusted for covariates, and Model 2 additionally adjusted for other variables of paternal childcare to examine how each of the paternal childcare variables, that is, paternal childcare hours on weekdays, weekends, and total caregiving scores, affected a child's behavioral problems. We also examined whether the frequency of each type of paternal caregiving was associated with each type of child behavioral problem after adjusting for paternal childcare hours.
All analyses were conducted separately for boys and girls because the effects of paternal childcare may differ due to biologically established sex differences in child development [27]. Analyses were performed using Stata version 15.0 (Stata Corp., College Station, TX, USA).
Approval from the ethics committee of the Tokyo Medical and Dental University was waived because the data were anonymous and available from the Ministry of Health, Labour, and Welfare in Japan upon request. Questionnaire responses from caregivers were considered as informed consent to participate in the study.