$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Final analytic dataset and participant characteristics
The final analytic dataset included 386 university students. Formal recruitment was conducted from 18 August 2025 to 26 September 2025 through university learning platforms, course communication groups, student email lists, and institutional student channels. A total of 432 submitted responses were screened, and 386 responses were retained, giving a valid-response rate of 89.4%. The locked dataset contained 33 substantive Likert-scale items and one instruction-based attention-check item. The attention-check item was used only for response-quality screening and was excluded from construct scoring, reliability analysis, validity analysis, PLS-SEM estimation, indirect-association analysis, and regression-based robustness analysis. The final sample included 223 female students, 154 male students, and nine participants who selected “prefer not to say” or another response option. There were 318 undergraduate students and 68 postgraduate students. Year-of-study distribution included 110 Year 1 students, 102 Year 2 students, 103 Year 3 students, and 71 students in Year 4 or above. Discipline categories included business and management, social sciences, engineering and technology, medicine and health sciences, humanities and arts, and education. Weekly online learning exposure included less than 2 h per week in 41 students, 2–4 h in 116 students, 5–7 h in 137 students, 8–10 h in 62 students, and more than 10 h in 30 students. Prior online learning experience included none in 38 students, less than 6 months in 68 students, 6–12 months in 190 students, and more than 12 months in 90 students. Weekly online learning hours and prior online learning experience were retained as prespecified control variables.
Expert review, pilot testing, and questionnaire verification
Expert review was completed by five experts in educational psychology, online learning, quantitative methodology, higher education research, and survey design. The review identified 13 wording issues, four redundancy issues, and three construct-alignment issues. Fourteen items were revised for clarity, contextual fit, or construct alignment, and no item was removed before pilot testing. Pilot testing was completed by 42 university students who met the main-study eligibility criteria. The median completion time was 9.6 min, and the minimum completion-time threshold was set at 3.2 min. Preliminary Cronbach’s alpha values ranged from 0.721 to 0.814. No item was removed after pilot testing. The final English-language questionnaire included 12 online learning experience items, nine basic psychological needs items, 12 learning engagement items, and one instruction-based attention-check item. Item-source mapping, variable coding, and scoring rules were locked before formal analysis and are provided in Supplementary File 1.
Response screening, missingness, and data-quality checks
After screening, 386 valid responses were retained from 432 submitted responses. Exclusions included four responses without consent confirmation, 17 eligibility exclusions, eight duplicate responses, six attention-check failures, five below-threshold completion-time responses, four straight-lining responses, and two unresolved eligibility inconsistencies. Six long-string patterns were flagged, but these overlapped with other exclusion categories and were not counted as an additional independent exclusion category. Exclusions were recorded in the de-identified screening log. No duplicate participant codes were retained in the locked dataset. No missing values were observed in the 33 substantive Likert-scale items, background variables, dimension-score variables, or global construct-score variables. No imputation was applied. No retained response met the straight-lining rule, no retained response showed a long-string pattern of 20 or more consecutive identical substantive responses, and no unresolved inconsistency was observed between online learning exposure and weekly online learning hours.
Descriptive distribution and common method bias diagnostics
The global online learning experience score was 3.47 ± 0.45. The global basic psychological needs score was 3.40 ± 0.49. The global learning engagement score was 3.42 ± 0.48. Instructional support had the highest mean among online learning experience dimensions, and interaction quality had the lowest mean. Competence need had the highest mean among basic psychological need dimensions, and relatedness need had the lowest mean. Cognitive engagement had the highest mean among learning engagement dimensions, and emotional engagement had the lowest mean. Item-distribution checks did not show severe floor or ceiling effects. The maximum item-level floor response rate was 1.55% for OLE_PU3. The maximum item-level ceiling response rate was 12.95% for OLE_IS2. Both values were below the prespecified 20% threshold. Common method bias diagnostics did not show severe evidence of a dominant common method factor. Harman’s single-factor test showed that the first unrotated factor explained 24.12% of the total variance. Full collinearity diagnostics produced a maximum variance inflation factor of 1.527.
Measurement-model quality
As shown in Table 3, Cronbach’s alpha values ranged from 0.733 for interaction quality to 0.801 for cognitive engagement at the dimension level. The global online learning experience construct showed an alpha of 0.878, basic psychological needs showed an alpha of 0.846, and learning engagement showed an alpha of 0.867. Standardized outer loadings ranged from 0.704 to 0.842 across the 33 substantive items. Composite reliability values ranged from 0.849 to 0.889. Average variance extracted values ranged from 0.653 to 0.702. The maximum heterotrait-monotrait ratio was 0.742. These values met the prespecified measurement-model criteria.
Table 3: Reliability, descriptive statistics, and measurement-model quality. This table reports Cronbach’s alpha, composite reliability, average variance extracted, standardized loading ranges, means, standard deviations, and retention decisions for the first-order dimensions and higher-order constructs. The attention-check item was excluded from reliability, validity, and model estimation. The maximum HTMT value among the first-order constructs was 0.742, below the prespecified conservative threshold of 0.85. The full outer-loading table and HTMT matrix are provided in Supplementary File 2. Abbreviations: AVE = average variance extracted. Please click here to download this Table.
Online learning experience was correlated with basic psychological needs (r = 0.501) and learning engagement (r = 0.497). Basic psychological needs were correlated with learning engagement (r = 0.510).
Structural pathway and indirect-association results
As shown in Figure 2, online learning experience was positively associated with basic psychological needs (β = 0.501, t = 11.259, p < 0.001, 95% CI: 0.414 to 0.588). Basic psychological needs were positively associated with learning engagement (β = 0.349, t = 7.256, p < 0.001, 95% CI: 0.255 to 0.444). Online learning experience remained directly associated with learning engagement (β = 0.324, t = 6.699, p < 0.001, 95% CI: 0.229 to 0.419). The model explained 25.1% of the variance in basic psychological needs and 33.9% of the variance in learning engagement (Table 4). The f2 effect size for the online learning experience on basic psychological needs was 0.332. The f2 effect sizes for the online learning experience and basic psychological needs on learning engagement were 0.118 and 0.138, respectively.

Figure 2: Structural model with standardized association-based path coefficients. The figure shows the standardized pathways among online learning experience, basic psychological needs, and learning engagement, including the direct path, indirect path, R2 values, and prespecified control variables. All paths are interpreted as association-based estimates. Please click here to view a larger version of this figure.
Table 4: Explained variance of endogenous constructs. This table reports the coefficient of determination (R2), adjusted R2 values, and predictor variables for each endogenous construct in the structural model. All results are interpreted as association-based estimates because the data were cross-sectional. Abbreviations: OLE_TOTAL = global online learning experience score; BPN_TOTAL = global basic psychological needs score; LE_TOTAL = global learning engagement score. Please click here to download this Table.
Table 5 shows that weekly online learning hours were not significantly associated with basic psychological needs (β = 0.001, p = 0.986) or learning engagement (β = −0.020, p = 0.634).
Table 5: Structural pathway estimates and effect sizes. This table reports standardized structural-pathway estimates, corresponding t values, p values, 95% confidence intervals (CI), f2 effect sizes, and pathway decisions for the hypothesized model. Control-variable pathways are also presented. All pathways are interpreted as association-based estimates because the data were cross-sectional. Abbreviations: OLE_TOTAL = global online learning experience score; BPN_TOTAL = global basic psychological needs score; LE_TOTAL = global learning engagement score; WEEKLY_HOURS = weekly online learning hours; PRIOR_EXPERIENCE = prior online learning experience; CI = confidence interval. Please click here to download this Table.
Prior online learning experience was not significantly associated with basic psychological needs (β = −0.005, p = 0.916) or learning engagement (β = 0.013, p = 0.754). The indirect pathway from online learning experience to learning engagement through basic psychological needs was supported. The indirect effect was 0.175, with a 95% bootstrap confidence interval from 0.124 to 0.231. The direct effect of the online learning experience on learning engagement was 0.324, and the total effect was 0.499 (Table 6). This pathway was interpreted as an association-based indirect pathway rather than causal mediation.
Table 6: Direct, indirect, and total association estimates. This table reports the direct, indirect (mediated), and total associations between online learning experience and learning engagement, including bootstrap 95% confidence intervals (CI) and pathway decisions. All effects are interpreted as association-based estimates because the data were cross-sectional. Abbreviations: OLE_TOTAL = global online learning experience score; BPN_TOTAL = global basic psychological needs score; LE_TOTAL = global learning engagement score; CI = confidence interval. Please click here to download this Table.
Psychological-need subdimension sensitivity analysis and robustness analysis
The indirect pathway through autonomy need was supported (indirect effect = 0.087, 95% CI: 0.049 to 0.131). The indirect pathway through competence need was supported and showed the largest indirect effect (indirect effect = 0.148, 95% CI: 0.094 to 0.205). The indirect pathway through relatedness need was also supported (indirect effect = 0.092, 95% CI: 0.056 to 0.132) (Table 7).
Table 7: Sensitivity analysis of basic psychological need subdimensions. This table reports indirect associations through the autonomy, competence, and relatedness subdimensions of basic psychological needs, including bootstrap 95% confidence intervals (CI) and pathway decisions. All effects are interpreted as association-based estimates because the data were cross-sectional. Abbreviations: OLE_TOTAL = global online learning experience score; LE_TOTAL = global learning engagement score; BPN_AU = autonomy need; BPN_CO = competence need; BPN_RE = relatedness need; CI = confidence interval. Please click here to download this Table.
Regression-based robustness analysis produced the same substantive interpretation as the primary PLS-SEM model. Online learning experience remained positively associated with basic psychological needs (β = 0.501, 95% CI: 0.414 to 0.588). Basic psychological needs remained positively associated with learning engagement after adjustment for online learning experience and control variables (β = 0.349, 95% CI: 0.255 to 0.444). The regression-based indirect effect was 0.175, with a 95% bootstrap confidence interval from 0.124 to 0.231. Because the data were cross-sectional, all pathways were interpreted as association-based estimates rather than causal effects (Table 8).
Table 8: Regression-based robustness analysis. This table reports regression-based robustness results for the structural pathways, including adjusted direct and indirect association estimates with bootstrap 95% confidence intervals (CI). All effects are interpreted as association-based estimates because the data were cross-sectional. Abbreviations: OLE_TOTAL = global online learning experience score; BPN_TOTAL = global basic psychological needs score; LE_TOTAL = global learning engagement score; CI = confidence interval. Please click here to download this Table.
Supplementary File 1: Questionnaire, variable coding, and scoring rules. This file provides the English-language questionnaire, eligibility-screening variables, background-variable coding, item wording, item-source mapping, response scale, scoring formulas, variable dictionary, data-handling rules, and consistency checks.Please click here to download this file.
Supplementary File 2: Final analytic dataset verification, analysis settings, and reproducibility records. This file records the locked analytic dataset, fieldwork summary, expert review and pilot-test records, screening summary, missing-data checks, reliability and measurement-model output, structural-model results, sensitivity analysis, robustness analysis, and reproducibility checklist.Please click here to download this file.