Skip to main content
Have a personal or library account? Click to login
Prevalence and associated risk factors of frozen shoulder in patients with shoulder musculoskeletal disorders: a cross-sectional analysis Cover

Prevalence and associated risk factors of frozen shoulder in patients with shoulder musculoskeletal disorders: a cross-sectional analysis

Open Access
|Sep 2026

Figures & Tables

Table 1.

Sample structure

QuestionFrequencyPercent
GenderN%
Male2735.5
Female4964.5
Age groupN%
18–3056.6
31–402127.6
41–502938.2
51–601722.4
>6045.3

[i] N - number of participants, % - percentage of participants

Table 2.

Prevalence of FS according to type of shoulder MSD

Type of shoulder disorderFS yesFSTotal
n (%)n (%)n (%)
Rotator cuff injury42 (95.5)2 (4.5)44 (57.9)
Impingement syndrome29 (90.6)3 (9.4)32 (42.1)
Total71 (93.4)5 (6.6)76 (100)
Pearson χ20.703df = 1p = 0.402
Fisher’s exact testp = 0.644

[i] FS - frozen shoulder, MSDs - musculoskeletal disorders

Table 3.

Descriptive statistics and correlation between pain severity and movement limitations

VariableMeanSDMinimumMaximum
Severity of shoulder pain6.532.32110
Movement limitations3.511.3916
Pearson correlation (r)0.685
p-value<0.001
Normality test (Shapiro–Wilk p-value)0.542–0.866 (all subgroups, p > 0.05)

[i] SD - standard deviation, r - Pearson correlation coefficient, p - level of statistical significance (p-value)

Table 4.

Logistic regression analysis of factors associated with FS: standard ML estimates and Firth penalised logistic regression estimates

Predictor (reference category)ML: B (SE)ML p-valueFirth: B (SE)OR (95% CI)Firth p-value
Gender: female vs male (ref = male)18.52 (2.97)0.0241.42 (0.63)4.14 (1.20–14.29)0.025
Age (ordinal, per category increase, 18–30 → >60)−18.23 (1.42)<0.001−0.98 (0.23)0.38 (0.24–0.59)<0.001
History of diabetes: yes vs no (ref = no)−16.15 (5218.97)0.998−0.78 (0.95)0.46 (0.07–2.97)0.412
History of shoulder trauma: yes vs no (ref = no)−3.76 (0.94)<0.001−2.05 (0.54)0.13 (0.045–0.37)<0.001
Model fit−2LL = 48.97Nagelkerke R2 = 0.497Firth − 2LL [nga rianaliza]Firth pseudo-R2 [nga rianaliza]
MulticollinearityTolerance 0.582–0.831VIF 1.203–1.718

[i] CI - confidence interval, ML - maximum-likelihood, OR - odds ratio, SE - standard error, VIF - variance inflation factor

Appendix 1.

Results from regression analysis presented in visual form. VIF - variance inflation factor

Note: The regression coefficients and p-values displayed in this figure correspond to the standard maximum-likelihood (ML) model in Table 4 and are shown here for visualization purposes only. Because the ML estimate for diabetes mellitus was unstable due to quasi-complete separation (p = 0.998), interpretation of this predictor in the text and conclusions is based on the Firth penalized logistic regression estimate (p = 0.412; OR = 0.46, 95% CI 0.07–2.97). The previous “Prob > F” statistic has been removed from the figure, as it is an F-test measure applicable to linear regression and is not an appropriate or applicable statistic for a logistic regression model; see Table 4 for the corresponding model fit statistics.

Language: English
Submitted on: Jun 16, 2026
Accepted on: Sep 7, 2026
Published on: Sep 29, 2026
Published by: Józef Piłsudski University of Physical Education in Warsaw
In partnership with: Paradigm Publishing Services
Publication frequency: 4 issues per year

© 2026 Józef Piłsudski University of Physical Education in Warsaw
This work is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.