Baseline characteristics of adolescents and parents are shown in Table 1. Mean age of adolescents (46.3% boys) was 15.22 ± 1.95 years with no significant difference between the intervention and control groups. The intervention group included more high-risk families (73.4%) than control (55.5%) groups (P < 0.001). Except for mothers’ smoking and fathers’ physical activity, other socio-behavioral variables were significantly different between study groups (p < 0.05). Details of parental risk clusters and their contributing factors have been reported in Fig. 1-Appendix and Table 1-Appendix. Paternal education was the most important factor in the clustering. In the high-risk cluster, most mothers and nearly 80% of fathers were illiterate or had primary education. 42% of fathers were smoker in the high-risk group, compared to 2.5% in low-risk cluster.
Table 1
Baseline characteristics of study participants in intervention and control groups
Adolescent characteristics | Total (n = 797) | Control (n = 605) | Intervention (n = 192) | p-value |
Age (year) | 15.22 ± 1.95 | 15.14 ± 1.95 | 15.50 ± 1.93 | 0.66 |
Sex | | | | 0.75 |
Boys | 369 (46.3%) | 282 (46.6%) | 87 (45.3%) | |
Girls | 428 (53.7%) | 323 (53.4%) | 105 (54.7%) | |
Maternal characteristics | | | | |
Age (year) | 42.12 ± 6.49 | 41.36 ± 6.16 | 44.51 ± 6.92 | 0.024 |
Education | | | | 0.003 |
Illiterate-primary | 351 (44%) | 250 (41.3%) | 101 (52.6%) | |
High school-diploma | 300 (37.6%) | 230 (38%) | 70 (36.5%) | |
Academic degree | 146 (18.3%) | 125 (20.7%) | 21 (10.9%) | |
Employment | | | | 0.002 |
Employed | 62 (7.8%) | 57 (9.4%) | 5 (2.6%) | |
Housewives | 735 (92.2%) | 548 (90.6%) | 187 (97.4%) | |
Smoking status | | | | 0.28 |
Yes | 36 (4.5%) | 30 (5%) | 6 (3.1%) | |
No | 757 (95%) | 571 (95%) | 186 (96.9%) | |
Physical activity* | | | | 0.003 |
Low | 190 (24.1%) | 137 (22.9%) | 53 (27.7%) | |
Moderate | 464 (58.9%) | 343 (57.5%) | 121 (63.4%) | |
High | 134 (17%) | 117 (19.6%) | 17 (8.9%) | |
Paternal characteristics | | | | |
Age (year) | 48.20 ± 7.59 | 47.35 ± 7.13 | 50.90 ± 8.34 | 0.003 |
Education | | | | 0.006 |
Illiterate-primary | 427 (53.6%) | 307 (50.7%) | 120 (62.5%) | |
High school-diploma | 324 (40.7%) | 257 (42.5%) | 67 (34.9%) | |
Academic degree | 46 (5.8%) | 41 (6.8%) | 5 (2.6%) | |
Employment | | | | < 0.001 |
Employed | 645 (80.9%) | 510 (84.3%) | 135 (70.3%) | |
Un-employed | 152 (19.1%) | 95 (15.7%) | 57 (29.7%) | |
Smoking status | | | | 0.008 |
Yes | 208 (26.1%) | 172 (28.4%) | 36 (18.8%) | |
No | 589 (73.9%) | 433 (71.6%) | 156 (81.3%) | |
Physical activity | | | | 0.48 |
Low | 296 (46.1%) | 219 (44.8%) | 77 (50.3%) | |
Moderate | 183 (28.5%) | 142 (29%) | 41 (26.8%) | |
High | 163 (25.4%) | 128 (26.2%) | 35 (22.9%) | |
Parental Cluster | | | | < 0.001 |
Low risk | 320 (40.2%) | 269 (44.5%) | 51 (26.6%) | |
High risk | 477 (59.8%) | 336 (55.5%) | 141 (73.4%) | |
* Physical activity was defined as low (< 600 MET-min/wk), moderate (600–3000 MET-min/wk), and high (≥ 3000 MET-min/wk). |
Figure 2 illustrates the prevalence of current cigarette smoking, passive smoking, and hookah smoking over the follow-ups. Cigarette smoking had no significant difference between the intervention and control, any time. Passive smoking was significantly more prevalent in control group at baseline, 1st and 2nd follow-ups. Although hookah smoking was always less prevalent in the intervention group, the difference was not significant in any follow-up.
The effects of intervention and time on tobacco smoking are shown in Table 2. Based on GEE results, the interaction effect of intervention*time was not significant on cigarette smoking (Wald’s Chi-Square = 0.75, DF = 4, P = 0.94); therefore, it was excluded in the final GEE models. Regardless of sex, the intervention significantly decreased the odds of cigarette smoking by 38% (OR = 0.62, 95% CI = 0.40–0.97; P = 0.04). Compared to baseline, the odds of cigarette smoking showed an increasing trend from the first (OR = 4.02, 95% CI = 2.31–7.02; P < 0.001) to the last (OR = 19.14, 95% CI = 10.89–33.66; P < 0.001) follow-up. High-risk families had higher odds (OR = 1.97, 95% CI = 1.35–2.88; P < 0.001) of cigarette smoking. Girls had lower odds of cigarette smoking, compared to boys (OR = 0.12, 95% CI = 0.08–0.19; P < 0.001). Sex-specific analysis showed that the positive effect of intervention in reducing the odds of cigarette smoking was seen only in boys (OR = 0.60, 95% CI = 0.37–0.99; P = 0.05) not in girls (OR = 0.69, 95% CI = 0.26–1.80; P = 0.45).
Table 2
The total and sex-specific effects of intervention on tobacco smoking in different follow-ups.
| | Cigarette smoking | | Passive smoking | | Hookah smoking* |
| Predictors | OR (95% CI) | P | | OR (95% CI) | P | | OR (95% CI) | P |
Total | Group | | | | | | | | |
Intervention | 0.62 (0.40–0.97) | 0.04 | | 0.43 (0.28–0.64) | < 0.001 | | 0.84 (0.62–1.13) | 0.25 |
Time | | | | | | | | |
1st Follow-up | 4.02 (2.31–7.02) | < 0.001 | | 0.67 (0.56–0.80) | < 0.001 | | Ref a | |
2nd Follow-up | 8.00 (5.25–15.41) | < 0.001 | | 0.46 (0.36–0.58) | < 0.001 | | 2.50 (1.99–3.15) | < 0.001 |
3rd Follow-up | 12.39 (7.06–21.73) | < 0.001 | | 0.32 (0.24–0.42) | < 0.001 | | 2.30 (1.83–2.90) | < 0.001 |
4th Follow-up | 19.14 (10.89–33.66) | < 0.001 | | 0.19 (0.13–0.28) | < 0.001 | | 2.19 (1.72–2.75) | < 0.001 |
parental cluster | | | | | | | . | . |
High-risk | 1.97 (1.35–2.88) | < 0.001 | | 3.50 (2.41–5.09) | < 0.001 | | 1.46 (1.12–1.90) | 0.005 |
Ado. Sex | | | | | | | | |
Girl | 0.12 (0.08–0.19) | < 0.001 | | 1.90 (1.40–2.60) | < 0.001 | | 0.37 (0.29–0.48) | < 0.001 |
Ado. Age | 1.11 (1.01–1.22) | 0.002 | | 0.96 (0.89–1.04) | 0.37 | | 1.01 (0.95–1.08) | 0.71 |
Boys | Group | | | | | | | | |
Intervention | 0.60 (0.37–0.99) | 0.05 | | 0.42 (0.23–0.75) | 0.003 | | 0.88 (0.59–1.29) | 0.51 |
Time | | | | | | | | |
1st Follow-up | 3.66 (2.07–6.45) | < 0.001 | | 0.53 (0.40–0.72) | < 0.001 | | Ref a | |
2nd Follow-up | 7.98 (4.64–13.74) | < 0.001 | | 0.23 (0.13–0.38) | < 0.001 | | 2.22 (1.62–3.06) | < 0.001 |
3rd Follow-up | 11.64 (6.56–20.66) | < 0.001 | | 0.13 (0.07–0.25) | < 0.001 | | 2.15 (1.56–2.94) | < 0.001 |
4th Follow-up | 16.03 (9.00–28.54) | < 0.001 | | 0.10 (0.04–0.22) | < 0.001 | | 2.03 (1.47–2.81) | < 0.001 |
parental cluster | | | | | | | . | . |
High-risk | 2.10 (1.36–3.24) | 0.001 | | 5.27 (2.70–10.31) | < 0.001 | | 1.68 (1.20–2.36) | 0.003 |
Ado. Age | 1.17 (1.04–1.30) | 0.005 | | 0.10 (0.88–1.13) | 0.99 | | 1.09 (1.00–1.18) | 0.04 |
Girls | Group | | | | | | | | |
Intervention | 0.69 (0.26–1.80) | 0.45 | | 0.43 (0.25–0.74) | 0.002 | | 0.75 (0.46–1.22) | 0.25 |
Time | | | | | | | | |
1st Follow-up | Ref | . | | 0.78 (0.63–0.95) | 0.01 | | Ref a | |
2nd Follow-up | 2.71 (1.20–6.10) | 0.02 | | 0.66 (0.51–0.85) | 0.001 | | 2.94 (2.08–4.16) | < 0.001 |
3rd Follow-up | 2.90 (1.16–7.27) | 0.02 | | 0.49 (0.35–0.65) | < 0.001 | | 2.60 (1.84–3.66) | < 0.001 |
4th Follow-up | 7.05 (2.78–17.91) | < 0.001 | | 0.27 (0.17–0.43) | < 0.001 | | 2.50 (1.72–3.62) | < 0.001 |
parental cluster | | | | | | | . | . |
High-risk | 1.45 (0.70–3.04) | 0.31 | | 2.93 (1.86–4.63) | < 0.001 | | 1.17 (0.78–1.75) | 0.44 |
Ado. Age | 0.91 (0.79–1.07) | 0.26 | | 0.95 (0.85–1.05) | 0.28 | | 0.91 (0.83–1.00) | 0.07 |
* evaluating hookah smoking started from the 1st follow-up, therefore 1st follow-up was considered as baseline and three subsequent measurements were conducted over the study duration. |
The interaction effect of intervention*time was not significant on passive smoking (Wald’s Chi-Square = 0.97, DF = 3, P = 0.91), therefore excluded from the final GEE models. Intervention reduced the odds of passive smoking in the total population by 57% (OR = 0.43, 95%CI = 0.28–0.64; P < 0.001). Compared to baseline, the odds of passive smoking in the total population had a decreasing trend from the first (OR = 0.67, 95% CI = 0.56–0.80; P < 0.001) to the last (OR = 0.19, 95% CI = 0.13–0.28; P < 0.001) follow-up. The odds of being a passive smoker was significantly higher in high-risk families (OR = 3.50, 95% CI = 2.41–5.09; P < 0.001). Considering adolescents’ sex and age, while girls had significant higher odds of being a passive smoker (OR = 1.90, 95% CI = 0.1.40–2.60; P < 0.001), the effect of age was not significant (P = 0.37). Sex-specific analysis showed that the intervention reduced the odds of passive smoking in boys by 58% (OR = 0.42, 95% CI = 0.23–0.75; P = 0.003), and in girls by 57% (OR = 0.43, 95% CI = 0.25–0.74; P = 0.002).
In terms of hookah smoking, the interaction effect of intervention*time was not significant (Wald’s Chi-Square = 2.13, DF = 3, P = 0.55) and was excluded from the final GEE models. However, the intervention was not successful in reducing the risk of hookah smoking in adolescents (OR = 0.84, 95% CI = 0.62–1.13; P = 0.25). Living in high-risk families increased the risk of hookah smoking (OR = 1.41, 95% CI = 1.13–1.75; P = 0.002). Girls (OR = 2.30, 95% CI = 1.85–2.58; P < 0.001) had lower odds of hookah smoking; while, age had no significant effect. Moreover, after sex-specific analysis, intervention did not reduce hookah use in boys (OR = 0.88, 95% CI = 0.59–1.29; P = 0.51) and girls (OR = 0.75, 95% CI = 0.46–1.22; P = 0.25).