Descriptive Results
Descriptive statistics and correlations between study variables are reported in Table
2. Participants who were romantically active in their adolescence reported having had their first relationship in middle adolescence, had more than one romantic relationship on average, and spent around 24 total months of their adolescence in romantic relationships. These three variables were significantly correlated with each other: The younger participants were at their first relationship, the more partners and the longer the total length of romantic involvement they reported by age 20.
Table 2
Descriptive characteristics and pearson correlations among the study variables
Romantic involvement |
1. Age of first relationship | – | – | – | – | – | – | – |
2. Number of partners | −0.54*** | – | – | – | – | – | – |
3. Months in relationships | −0.53*** | 0.36*** | – | – | – | – | – |
Psychosocial adjustment |
4. Life satisfaction | 0.03 | −0.04 | 0.03 | – | – | – | – |
5. Self-esteem | 0.03 | −0.03 | −0.04 | 0.61*** | – | – | – |
6. Loneliness | −0.06* | 0.00 | −0.05* | −0.50*** | −0.61*** | – | – |
7. Depressive symptoms | −0.02 | 0.03 | 0.02 | −0.68*** | −0.76*** | 0.65*** | – |
M | 16.25 | 1.93 | 24.31 | 7.75 | 3.91 | 2.09 | 1.75 |
SD | 1.85 | 1.12 | 16.48 | 1.02 | 0.65 | 0.85 | 0.37 |
The psychosocial adjustment variables were also all significantly correlated with each other: Both the correlations between life satisfaction and self-esteem and between loneliness and depressive symptoms were positive. In evaluating the correlations between romantic relationship indicators and psychosocial adjustment, loneliness was found to be related to two of the indicators: The later participants started dating and the more time they spent in relationships, the less lonely they felt.
Latent Profile Analysis
Table
3 presents the fit statistics for the LPA models ranging from two to six classes. The model with six classes could not be properly identified, as the best log likelihood values in the model estimation could not be replicated and estimates were unreliable. Out of the remaining models, the three-class solution was chosen for the final model for the following four reasons: First, although each fit statistic decreased across the two- to the five-class solution, the smallest decrease was found when moving from the three- to the four-class solution, suggesting minimal improvement when a fourth class was included. Second, the LMR comparing the three- to the four-class model was not significant, again suggesting that a model with four classes did not fit the data better than the model with three classes. Third, beyond the solution of three classes, the sample size of the smallest group lay below the 5% minimum recommended by Nylund-Gibson and Choi (
2018), as this would limit statistical power for further analyses. Fourth, when comparing the distribution of romantic relationship indicators in the three- and the four-class solutions, the additional fourth class was found to be conceptually redundant to one of the other three classes.
Table 3
Fit statistics for LPA models ranging from two to six classes
2 classes | 28632 | 28697 | 28687 | 28656 | p < 0.001 | p < 0.001 | 0.69 | 42.19% |
3 classes | 28140 | 28232 | 28218 | 28174 | p = 0.015 | p < 0.001 | 0.77 | 9.14% |
4 classes | 27854 | 27972 | 27954 | 27897 | p = 0.404 | p < 0.001 | 0.83 | 3.08% |
5 classes | 27380 | 27525 | 27503 | 27433 | p = 0.022 | p < 0.001 | 0.95 | 3.08% |
6 classes | Model was not well identifieda |
After deciding on the final model, individuals were assigned to classes based on the highest affiliation probability. The entropy score for the final model indicated good classification accuracy. In addition to the three classes covering romantic involvement during adolescence, a fourth class was included for those participants who remained single during their adolescence. Descriptive information on the four classes is shown in Table
4. The final number of classes was in line with the first hypothesis.
Table 4
Descriptive information of the classes
Sample size | 907 | 833 | 175 | 542 |
% females | 0.48 | 0.58 | 0.50 | 0.40 |
Romantic involvement |
Age of first rel. M (SD) | 17.67 (1.29)a | 15.12 (1.18)b | 14.32 (1.37)c | – |
No. of partners M (SD) | 1.19 (0.40)a | 2.20 (0.69)b | 4.5 (0.88)c | 0 (0) |
Months in rel. M (SD) | 13.83 (10.25)a | 33.56 (15.43)b | 34.99 (14.37)b | 0 (0) |
A MANOVA comparing the three classes from the LPA with regard to romantic experiences proved to be significant, Wilks’ Λ = 0.16, F(2, 1903) = 969.81, p < 0.001. Tukey’s post hoc test revealed significant differences between the three groups on all of the romantic relationship variables: age of first relationship, F(2, 1912) = 1138, p < 0.001, number of romantic partners, F(2, 1912) = 2425, p < 0.001, and romantic involvement in months, F(2, 1903) = 552.5, p < 0.001.
The first and largest class was labeled late starters (36.91%) because this group was characterized by a later age of entering the first romantic relationship than all the other groups (all p < 0.001). In addition, this group showed the lowest number of partners and lowest overall romantic involvement duration (all p < 0.001).
The smallest group was called frequent changers (7.12%). This group presented the highest number of romantic partners among all groups (all p < 0.001) and was also characterized by the earliest age of entering one’s first relationship (all p < 0.001).
The second-largest group comprised 33.90% of the sample and was named moderate daters, as this group lay in between the late starters and frequent changers. Although this group spent as many months in romantic relationships as the frequent changers (p = 0.390), participants in this group started dating at a later age (p < 0.001) and had significantly fewer partners (p < 0.001) than the frequent changers. However, they started at an earlier age (p < 0.001) and reported more partners (p < 0.001) than the late starters.
The group of those who reported not having had a romantic relationship by age of 20 was termed continuous singles and comprised 22.06% of participants.
Although these four groups differed with regard to sex, χ2(3) = 33.19, p < 0.001—More female participants were classified into the moderate daters group and more male participants into the singles group—the effect size was small (Cramer’s V = 0.13).
First-Order Latent Basis Growth Curve Models
Results of the latent basis growth curve models for each outcome separated by group can be found in Table
5.
Table 5
Parameter estimates of the latent growth curve models for each outcome separated by group
Life satisfaction |
Intercept W1 | 7.95 | 0.04 | <0.001 | 8.06a | 0.04 | <0.001 | 7.79 | 0.09 | <0.001 | 7.77b | 0.05 | <0.001 |
Slope | −0.27 | 0.05 | <0.001 | −0.36 | 0.05 | <0.001 | −0.39 | 0.11 | <0.001 | −0.38 | 0.06 | <0.001 |
Intercept W10 | 7.68a | 0.04 | <0.001 | 7.70a | 0.05 | <0.001 | 7.40 | 0.09 | <0.001 | 7.40b | 0.05 | <0.001 |
Self-esteem |
Intercept W2 | 3.99 | 0.03 | <0.001 | 3.99 | 0.03 | <0.001 | 3.93 | 0.06 | <0.001 | 3.93 | 0.03 | <0.001 |
Slope | −0.12 | 0.03 | <0.001 | −0.11 | 0.03 | <0.001 | −0.20 | 0.07 | 0.003 | −0.13 | 0.04 | 0.001 |
Intercept W10 | 3.87 | 0.03 | <0.001 | 3.88 | 0.03 | <0.001 | 3.73 | 0.06 | <0.001 | 3.80 | 0.03 | <0.001 |
Loneliness |
Intercept W4 | 2.05a | 0.03 | <0.001 | 1.93a | 0.04 | <0.001 | 2.06 | 0.08 | <0.001 | 2.24b | 0.04 | <0.001 |
Slope | 0.01 | 0.04 | 0.792 | 0.16 | 0.05 | <0.001 | 0.14 | 0.10 | 0.160 | 0.15 | 0.06 | 0.008 |
Intercept W10 | 2.06a | 0.04 | <0.001 | 2.09a | 0.04 | <0.001 | 2.20 | 0.08 | <0.001 | 2.39b | 0.05 | <0.001 |
Depressive symptoms |
Intercept W2 | 1.68 | 0.01 | <0.001 | 1.65 | 0.01 | <0.001 | 1.72 | 0.03 | <0.001 | 1.68 | 0.02 | <0.001 |
Slope | 0.14 | 0.02 | <0.001 | 0.16 | 0.02 | <0.001 | 0.17 | 0.03 | <0.001 | 0.16 | 0.02 | <0.001 |
Intercept W10 | 1.82 | 0.02 | <0.001 | 1.80 | 0.02 | <0.001 | 1.89 | 0.03 | <0.001 | 1.84 | 0.02 | <0.001 |
Life satisfaction
The latent growth curve model for life satisfaction provided a good fit with RMSEA = 0.04, CFI = 0.95, and TLI = 0.95. With regard to initial group differences in life satisfaction, continuous singles tended to show lower levels of initial life satisfaction compared to those in the moderate daters group (d = −0.29, p < 0.001). There were no other significant differences in the intercepts between the groups. The slope across measurements was significantly negative for each group, suggesting a general decrease in life satisfaction over time. In young adulthood, the significant difference between singles and moderate daters was still present (d = −0.29, p < 0.001) and in addition, singles also reported lower life satisfaction than the late starters in young adulthood (d = −0.27, p < 0.001).
Self-esteem
The latent growth curve model for self-esteem suggested a good fit (RMSEA = 0.04, CFI = 0.98, TLI = 0.97). Results indicated no significant group differences in either the intercepts or the slopes. Independent of group affiliation, self-esteem decreased significantly across the years.
Loneliness
The model fit for loneliness was good with RMSEA = 0.03, CFI = 0.98, and TLI = 0.98. Group comparisons indicated higher initial loneliness levels among the continuous singles compared to both the moderate daters (d = 0.43, p < 0.001) and late starters (d = 0.26, p < 0.001). These differences were also present in young adulthood with d = 0.38, p < 0.001 for continuous singles vs. moderate daters, and d = 0.41, p < 0.001 for continuous singles vs. late starters, respectively. The magnitude of change in loneliness was equal in all groups and indicated a general increase in loneliness over time.
Depressive symptoms
The RMSEA = 0.05, CFI = 0.96, and TLI = 0.95 suggested a good fit for the latent growth curve model for depressive symptoms. Adolescents’ depressive symptoms and their change over time did not differ across groups. In all groups, depressive symptoms increased over time through young adulthood.
Together, these results suggest that part of the hypothesis regarding the link between romantic involvement and psychosocial adjustment was supported. In particular, those with no romantic experiences in adolescence indicated lower life satisfaction and more loneliness than those who dated moderately. However, the group characterized by frequent dating did not differ from the others in these aspects and no group differences were found in self-esteem or depressive symptoms. Regarding the exploratory question, although romantic involvement had no impact on the rate of change of psychosocial adjustment, lasting associations were found for life satisfaction and loneliness in young adulthood.
Sensitivity Analyses
Two sensitivity analyses were conducted to test the robustness of the results. First, parametric bootstrapping with 1000 replications was applied to the growth curve models. The results of these analyses fully replicated those of the main analyses. Second, to ensure that findings were not sensitive to missing data, all models were further estimated utilizing solely participants with complete data. Again, the results of the main analyses were largely confirmed. All statistically significant effects found in the main analyses were fully replicated. However, two additional group differences emerged: Late starters reported significantly higher life satisfaction than continuous singles (d = 0.36, p < 0.001) at the first wave, and continuous singles reported more depressive symptoms than the moderate daters (d = 0.30, p = 0.001) at the last wave.
Overall, the results of these sensitivity analyses largely replicated the main findings. In the following section, only the robust results of the main analyses will be discussed. Finally, it should be noted that all data analyses were conducted as planned and no participants or variables were excluded due to lack of significance.