Study Characteristics and Measurement Framework
Environmental noise measurements were performed in 15 inpatient departments of a large tertiary teaching hospital between March and August 2025 using a standardized hierarchical sampling protocol (Figure 1). Within each department, four predefined functional zones (ward, patient activity area, corridor, and nurses’ station) were evaluated. Three fixed sampling locations were established within each functional zone, and measurements were obtained during three standardized daytime observation periods (09:00–19:00). At each location and observation period, three consecutive technical replicate measurements were recorded using a calibrated Class 1 sound level meter, generating acoustic data for LAeq, LAFmax, and LCpeak. Operational variables, including bed capacity, occupied beds, staff count, visitor count, alarm events, conversation events, room area, and room volume, were recorded concurrently with each acoustic observation. The planned sampling framework comprised 1,620 acoustic observations (15 departments × 4 functional zones × 3 sampling locations × 3 observation periods × 3 technical replicates). One scheduled observation in the patient activity area of the Intensive Care II Department could not be completed because of temporary emergency-care access restrictions and was handled using complete-case analysis without data imputation, as prespecified in the statistical analysis plan. Additional indoor and outdoor public-area measurements were obtained for descriptive contextual comparison but were excluded from the departmental mixed-effects analyses (Figure 1). The characteristics of the participating departments are summarized in Table 1. Bed capacity ranged from 21 to 72 beds, with corresponding mean occupied beds ranging from 17.8 to 60.6. Mean staff counts were relatively consistent across departments (5.5–5.9 personnel per observation), whereas visitor counts averaged approximately five individuals per observation. Departmental room areas ranged from 145.8 to 500.3 m2, and room volumes ranged from 466.5 to 1,601.0 m3, reflecting substantial variation in the physical characteristics of the participating departments.

Figure 1: Study Design and Hierarchical Environmental Noise Measurement Protocol. Overview of the study design and hierarchical environmental noise measurement protocol. Environmental noise was assessed in 15 inpatient departments of a tertiary teaching hospital using a multistage hierarchical sampling framework. Measurements were obtained from four standardized functional zones (ward, corridor, nurses’ station, and patient activity area), three predefined sampling locations (entrance, middle, and end), three daytime observation periods (morning, midday, and afternoon), and three consecutive technical replicates. Equivalent continuous A-weighted sound pressure level (LAeq) was the primary outcome, and maximum A-weighted sound pressure level (LAFmax) and peak C-weighted sound pressure level (LCpeak) were recorded as secondary acoustic outcomes. Operational variables were recorded concurrently, and all measurements underwent calibration and predefined quality-control procedures before statistical analysis. LAeq, equivalent continuous A-weighted sound pressure level; LAFmax, maximum A-weighted sound pressure level; LCpeak, peak C-weighted sound pressure level. Please click here to view a larger version of this figure.
| Department | Bed Capacity | Mean Occupied Beds | Room Area (m²) | Room Volume (m³) | Mean Staff Count | Mean Visitor Count | Mean Alarm Count | Mean Conversation Count |
| Trauma Orthopedics | 40 | 33.8 | 280.1 | 896.3 | 5.5 | 4.8 | 2.4 | 6.5 |
| Breast and Thyroid Diagnosis and Treatment Center | 60 | 50.6 | 416.8 | 1333.8 | 5.7 | 4.8 | 2.4 | 6.2 |
| Congenital Heart Disease Center | 35 | 29.4 | 243.6 | 779.5 | 5.7 | 5 | 2.4 | 6.5 |
| Otolaryngology–Cervicology Surgery | 50 | 42.2 | 348.1 | 1113.9 | 5.8 | 5 | 2.5 | 6.7 |
| Gastroenterology Diagnosis and Treatment Department I | 38 | 31.9 | 264.8 | 847.3 | 5.8 | 5 | 2.5 | 6.6 |
| Gastroenterology Diagnosis and Treatment Department II | 52 | 43.9 | 361.7 | 1157.4 | 5.7 | 4.9 | 2.4 | 6.5 |
| Hepatobiliary and Pancreatic Surgery | 48 | 40.5 | 334.1 | 1069.2 | 5.8 | 5.1 | 2.5 | 6.8 |
| Intensive Care II Department | 21 | 17.8 | 145.8 | 466.5 | 5.7 | 4.8 | 2.5 | 6.5 |
| Knee Joint Department II | 45 | 37.8 | 312.7 | 1000.6 | 5.8 | 4.9 | 2.5 | 6.5 |
| Parkinson's Disease Department | 55 | 46.2 | 382.2 | 1223 | 5.9 | 5.1 | 2.5 | 6.8 |
| Psychosomatic Medicine Department I | 68 | 57.1 | 472.6 | 1512.3 | 5.8 | 5 | 2.5 | 6.7 |
| Pediatrics Department | 40 | 33.6 | 278.6 | 891.5 | 5.8 | 5.1 | 2.5 | 6.8 |
| Urology Department | 62 | 52.1 | 431.2 | 1380 | 5.8 | 5 | 2.5 | 6.6 |
| Vascular and Nerve Department I | 72 | 60.6 | 500.3 | 1601 | 5.9 | 5.1 | 2.5 | 6.9 |
| Vascular and Nerve Department II | 60 | 50.4 | 417.5 | 1336 | 5.7 | 4.8 | 2.4 | 6.3 |
Table 1: Characteristics of Hospital Departments Included in the Study. Baseline operational characteristics of the 15 inpatient departments included in the environmental noise survey. Department-level variables include bed capacity, mean occupied beds, room area, room volume, and the mean staff count, visitor count, alarm count, and conversation count recorded during standardized daytime observation periods. Operational variables represent average values obtained during the predefined observation periods and were evaluated as exploratory variables in the mixed-effects regression analyses. Room area and room volume are reported in square meters (m2) and cubic meters (m3), respectively.
Descriptive Acoustic Characteristics Across Functional Zones
Across all 1,620 observations, nurses’ stations demonstrated the highest environmental noise levels, followed sequentially by corridors, patient activity areas, and wards (Table 2; Figures 2 and 3). Mean LAeq values were 66.90 ± 2.34 dB(A) at nurses’ stations, 63.97 ± 2.38 dB(A) in corridors, 61.38 ± 2.17 dB(A) in patient activity areas, and 60.38 ± 2.22 dB(A) in wards. Similar spatial gradients were observed for LAFmax and LCpeak, indicating consistent differences across all three-acoustic metrics. Detailed descriptive summaries of the secondary acoustic outcomes (LAFmax and LCpeak) are provided in Supplementary Table 1.
| Functional Zone | Observations (n) | LAeq, dB(A) Mean ± SD | LAFmax, dB(A) Mean ± SD | LCpeak, dB(C) Mean ± SD | Median LAeq (IQR), dB(A) |
| Ward | 405 | 60.38 ± 2.22 | 67.95 ± 2.37 | 79.21 ± 3.31 | 60.4 (58.9–61.9) |
| Corridor | 405 | 63.97 ± 2.38 | 71.46 ± 2.59 | 82.76 ± 3.56 | 64.0 (62.5–65.5) |
| Nurses' Station | 405 | 66.90 ± 2.34 | 74.43 ± 2.56 | 85.72 ± 3.56 | 66.7 (65.4–68.3) |
| Patient Activity Area | 405 | 61.38 ± 2.17 | 69.00 ± 2.42 | 80.20 ± 3.37 | 61.5 (60.0–62.7) |
Table 2: Descriptive Acoustic Characteristics by Hospital Functional Zone. Descriptive acoustic characteristics for the four standardized hospital functional zones. Values are presented as mean ± standard deviation (SD) unless otherwise indicated. The median and interquartile range (IQR) are reported for the equivalent continuous A-weighted sound pressure level (LAeq). Each functional zone contributed 405 observations to the hierarchical dataset (total n = 1,620). LAeq, equivalent continuous A-weighted sound pressure level; LAFmax, maximum A-weighted sound pressure level; LCpeak, peak C-weighted sound pressure level; IQR, interquartile range.

Figure 2: Heat Map of Mean Equivalent Continuous A-Weighted Sound Pressure Levels Across Hospital Departments and Functional Zones. Heat map showing mean equivalent continuous A-weighted sound pressure levels (LAeq) for each hospital department and functional zone. Cell values represent mean LAeq in A-weighted decibels [dB(A)], and the color scale indicates increasing sound levels from lower to higher values. Functional zones include wards, corridors, nurses’ stations, and patient activity areas. CHD, congenital heart disease; ICU, intensive care unit; LAeq, equivalent continuous A-weighted sound pressure level. Please click here to view a larger version of this figure.

Figure 3: Distribution of Equivalent Continuous A-Weighted Sound Pressure Levels Across Standardized Hospital Functional Zones. Boxplots showing the distribution of equivalent continuous A-weighted sound pressure levels (LAeq) across four standardized hospital functional zones. Boxes represent the interquartile range (25th–75th percentile), the horizontal line within each box indicates the median, and the × symbol denotes the mean. Whiskers extend to 1.5 times the interquartile range (IQR), and observations beyond the whiskers are displayed as outliers. Each functional zone includes 405 observations (total n = 1,620). Differences among functional zones were evaluated using the primary linear mixed-effects model with Tukey-adjusted pairwise comparisons. LAeq, equivalent continuous A-weighted sound pressure level; IQR, interquartile range. Please click here to view a larger version of this figure.
The heat map presented in Figure 2 demonstrates a consistent spatial pattern across the 15 participating departments. Nurses’ stations consistently exhibited the highest mean LAeq values irrespective of clinical specialty, whereas ward areas demonstrated the lowest equivalent continuous A-weighted sound pressure levels. Departments with larger bed capacities, including Vascular and Nerve Department I, Psychosomatic Medicine Department I, and the Urology Department, generally exhibited higher mean LAeq values across functional zones than departments with smaller bed capacities, such as the Intensive Care II Department and the Congenital Heart Disease Center. Although modest interdepartmental variability was observed, the relative ranking of functional zones remained consistent. Mean LAeq values at nurses’ stations ranged from 64.1 dB(A) in the Intensive Care II Department to 68.9 dB(A) in Vascular and Nerve Department I, whereas ward measurements ranged from 57.9 to 61.6 dB(A). Patient activity areas and corridors exhibited intermediate values across nearly all departments (Figure 2). Figure 3 illustrates the distribution of LAeq measurements within each functional zone. Median LAeq values were highest at nurses’ stations (66.7 dB[A]), followed by corridors (64.0 dB[A]), patient activity areas (61.5 dB[A]), and wards (60.4 dB[A]). Although occasional high-value observations were identified across all functional zones, the interquartile ranges were relatively narrow, indicating limited variability within each functional environment. The boxplots demonstrate clear separation in the distribution of LAeq values among the four functional zones (Figure 3).
Mixed-Effects Comparison of Functional Zones
To account for the hierarchical structure of repeated acoustic measurements within hospital departments, the primary analysis used a linear mixed-effects model with hospital department specified as a random intercept and functional zone as a fixed effect (Table 3). This modeling approach accounted for within-department clustering while estimating adjusted differences in LAeq among the four predefined functional zones.
| Fixed Effects (Outcome: LAeq) |
| Fixed Effect | β | SE | 95% CI | P Value |
| Intercept (Corridor) | 63.972 | 0.266 | 63.451 to 64.493 | <0.001 |
| Nurses' Station | 2.931 | 0.146 | 2.644 to 3.218 | <0.001 |
| Patient Activity Area | −2.594 | 0.146 | −2.881 to −2.307 | <0.001 |
| Ward | −3.589 | 0.146 | −3.876 to −3.302 | <0.001 |
| Estimated Marginal Means |
| Functional Zone | Estimated Mean LAeq, dB(A) |
| Nurses' Station | 66.9 |
| Corridor | 63.97 |
| Patient Activity Area | 61.38 |
| Ward | 60.38 |
| Tukey-Adjusted Pairwise Comparisons |
| Comparison | Mean Difference, dB(A) | Adjusted P Value |
| Nurses' Station vs. Corridor | 2.93 | <0.001 |
| Nurses' Station vs. Patient Activity Area | 5.52 | <0.001 |
| Nurses' Station vs. Ward | 6.52 | <0.001 |
| Corridor vs. Patient Activity Area | 2.59 | <0.001 |
| Corridor vs. Ward | 3.59 | <0.001 |
| Patient Activity Area vs. Ward | 0.99 | <0.001 |
| Model Performance |
| Statistic | Value |
| Intraclass Correlation Coefficient (ICC) | 0.172 |
| Marginal R² | 0.53 |
| Conditional R² | 0.61 |
Table 3: Primary Linear Mixed-Effects Model Comparing Environmental Noise Across Hospital Functional Zones. Results of the primary linear mixed-effects model evaluating differences in equivalent continuous A-weighted sound pressure levels (LAeq) among the four standardized hospital functional zones. Hospital department was specified as a random intercept to account for clustering of repeated observations, and functional zone was included as a fixed effect with corridor as the reference category. Regression coefficients (β), standard errors (SE), Wald 95% confidence intervals (CI), and two-sided P values are presented for the fixed effects. Estimated marginal means were compared using Tukey-adjusted pairwise comparisons. Model performance is summarized using the intraclass correlation coefficient (ICC) and the marginal and conditional coefficients of determination (R2), calculated according to the method of Nakagawa and Schielzeth. LAeq, equivalent continuous A-weighted sound pressure level; β, regression coefficient; SE, standard error; CI, confidence interval; ICC, intraclass correlation coefficient.
The mixed-effects analysis demonstrated a highly significant overall effect of functional zone on environmental noise levels (overall likelihood-ratio test, P < 0.001). Using the corridor as the reference category, nurses’ stations exhibited significantly higher LAeq values (β = 2.931 dB[A], SE = 0.146, 95% CI = 2.644–3.218, P < 0.001), whereas patient activity areas (β = −2.594 dB[A], SE = 0.146, 95% CI = −2.881 to −2.307, P < 0.001) and wards (β = −3.589 dB[A], SE = 0.146, 95% CI = −3.876 to −3.302, P < 0.001) demonstrated significantly lower LAeq values (Table 3). Estimated marginal means derived from the mixed-effects model demonstrated a clear hierarchical pattern of environmental noise across the functional zones. The adjusted mean LAeq was highest at nurses’ stations (66.90 dB[A]), followed by corridors (63.97 dB[A]), patient activity areas (61.38 dB[A]), and wards (60.38 dB[A]). All pairwise comparisons remained statistically significant after Tukey adjustment (all adjusted P < 0.001), indicating statistically significant differences in LAeq among all four functional zones. Mean differences ranged from 0.99 dB(A) between patient activity areas and wards to 6.52 dB(A) between nurses’ stations and wards. Figure 4 summarizes the model-adjusted estimates. Estimated marginal means and their 95% confidence intervals showed minimal overlap among the functional zones, and Tukey grouping assigned each zone to a unique significance group (A–D). Nurses’ stations were assigned to group A with the highest adjusted LAeq, followed by corridors (group B), patient activity areas (group C), and wards (group D). Variance component analysis demonstrated that approximately 17.2% of the total variability in LAeq was attributable to differences between hospital departments (intraclass correlation coefficient = 0.172), supporting inclusion of a department-level random intercept. The fixed effects explained 53% of the variability in environmental noise (marginal R2 = 0.53), whereas the combined fixed- and random-effects model explained 61% of the total variability (conditional R2 = 0.61), indicating good explanatory performance of the hierarchical model (Table 3).

Figure 4: Estimated Marginal Means of Equivalent Continuous A-Weighted Sound Pressure Levels Across Hospital Functional Zones. Estimated marginal means (EMMs) of equivalent continuous A-weighted sound pressure levels (LAeq) across the four hospital functional zones derived from the primary linear mixed-effects model. Points represent model-estimated marginal means, and horizontal error bars indicate 95% confidence intervals (CIs). Hospital department was included as a random effect and functional zone as a fixed effect. Pairwise comparisons were performed using Tukey-adjusted estimated marginal means; functional zones assigned different letters represent statistically significant differences (P < 0.05). Each functional zone contributed 405 observations to the analysis. LAeq, equivalent continuous A-weighted sound pressure level; EMM, estimated marginal mean; CI, confidence interval. Please click here to view a larger version of this figure.
Parsimonious Mixed-Effects Regression of Operational Characteristics Associated with LAeq
To investigate operational characteristics associated with environmental noise while minimizing multicollinearity, a prespecified parsimonious mixed-effects model was fitted with functional zone, visitor count, and departmental bed capacity included as fixed effects and hospital department retained as a random intercept (Table 4). Selection of these variables was informed by the exploratory correlation matrix (Supplementary Figure 1), which demonstrated strong correlations among several structural variables, particularly room area and room volume (Spearman’s ρ = 0.94). Consequently, highly correlated structural variables were excluded from the parsimonious model to improve interpretability and model stability. After adjustment for the operational covariates, functional zone remained the strongest independent factor associated with environmental noise. Relative to corridors, nurses’ stations were associated with an adjusted increase of 2.933 dB(A) (SE = 0.146, 95% CI = 2.647–3.220, P < 0.001), whereas patient activity areas (β = −2.602 dB[A], SE = 0.146, 95% CI = −2.888 to −2.315, P < 0.001) and wards (β = −3.606 dB[A], SE = 0.146, 95% CI = −3.892 to −3.319, P < 0.001) remained significantly quieter than corridors.
| Predictor | β Coefficient | SE | 95% CI | P Value |
| Intercept (Corridor) | 63.97 | 0.27 | 63.45 to 64.49 | <0.001 |
| Nurses' Station (vs. Corridor) | 2.933 | 0.146 | 2.647 to 3.220 | <0.001 |
| Patient Activity Area (vs. Corridor) | −2.602 | 0.146 | −2.888 to −2.315 | <0.001 |
| Ward (vs. Corridor) | −3.606 | 0.146 | −3.892 to −3.319 | <0.001 |
| Visitor Count (per additional visitor) | 0.06 | 0.026 | 0.008 to 0.112 | 0.024 |
| Bed Capacity (per additional bed) | 0.064 | 0.005 | 0.055 to 0.073 | <0.001 |
| Model: LAeq ~ Functional Zone + Visitor Count + Bed Capacity + (1 | Department) |
Table 4: Parsimonious Linear Mixed-Effects Regression of Operational Characteristics Associated with Environmental Noise. Results of the prespecified parsimonious linear mixed-effects regression model evaluating associations between equivalent continuous A-weighted sound pressure level (LAeq) and selected operational characteristics. The model included functional zone, visitor count, and departmental bed capacity as fixed effects and hospital department as a random intercept to account for clustering of repeated observations (model: LAeq ~ Functional Zone + Visitor Count + Bed Capacity + (1 | Department)). Corridor was specified as the reference category for functional zone, and visitor count and bed capacity were modeled as continuous variables. Regression coefficients (β), standard errors (SE), Wald 95% confidence intervals (CI), and two-sided P values are reported. Regression coefficients for visitor count and bed capacity represent the estimated change in LAeq associated with each additional visitor and each additional hospital bed, respectively. The intercept represents the estimated mean LAeq for the reference category (corridor) when continuous predictors are equal to zero. This table presents the prespecified parsimonious model developed to minimize model overfitting; the fully adjusted sensitivity model including all recorded operational characteristics is presented in Supplementary Table S2. LAeq, equivalent continuous A-weighted sound pressure level; β, regression coefficient; SE, standard error; CI, confidence interval.
Among the continuous operational predictors, visitor count demonstrated an independent positive association with LAeq. Each additional visitor present during measurement was associated with an estimated increase of 0.060 dB(A) (SE = 0.026, 95% CI = 0.008–0.112, P = 0.024). Similarly, departmental bed capacity remained independently associated with environmental noise, with each additional hospital bed associated with an estimated 0.064 dB(A) increase in LAeq (SE = 0.005, 95% CI = 0.055–0.073, P < 0.001).
Figure 5 illustrates the fitted predictions generated from the parsimonious mixed-effects model. Predicted LAeq increased approximately linearly across the observed range of departmental bed capacities, and higher visitor counts consistently shifted the predicted sound levels upward across the entire capacity range. The parallel fitted regression lines indicate additive effects of visitor count and bed capacity within the fitted model after adjustment for functional zone and department-level clustering. Model-based predictions demonstrated that, for a department with approximately 20 beds, the predicted LAeq ranged from 61.6 dB(A) at a low visitor count (three visitors) to 62.4 dB(A) at a higher visitor count (seven visitors). At the upper end of the observed bed-capacity range (72 beds), the corresponding predicted LAeq increased to 64.5 dB(A) and 65.3 dB(A), respectively, illustrating the combined contribution of departmental bed capacity and visitor count within the fitted model.

Figure 5: Model-Based Predictions of Equivalent Continuous A-Weighted Sound Pressure Levels According to Departmental Bed Capacity. Predicted equivalent continuous A-weighted sound pressure levels (LAeq) across the observed range of departmental bed capacity derived from the parsimonious linear mixed-effects regression model. Colored lines represent model-predicted LAeq values at low (3 visitors), average (5 visitors), and high (7 visitors) visitor counts while holding all other model terms constant. Shaded regions indicate 95% confidence intervals for the predicted values, and points represent observed measurements. Predictions are shown only within the observed range of departmental bed capacity included in the study. LAeq, equivalent continuous A-weighted sound pressure level; CI, confidence interval. Please click here to view a larger version of this figure.
Sensitivity analyses using the fully adjusted mixed-effects model (Supplementary Table 2), which additionally included occupied beds, room area, room volume, staff count, alarm events, conversation events, and measurement period, produced highly consistent estimates for functional zone and visitor count. Nurses’ stations remained the noisiest functional zone, and visitor count remained independently associated with higher LAeq (β = 0.060, P = 0.023), whereas occupied beds, room area, room volume, staff count, alarm events, conversation events, and measurement period were not independently associated with environmental noise after mutual adjustment. These findings support the robustness of the parsimonious model and indicate that functional zone and visitor count were the principal variables associated with environmental noise in this study.
Descriptive Comparison with Reference Acoustic Criteria
To place the observed environmental noise levels into context, measured LAeq values were descriptively compared with three widely cited reference acoustic criteria: the WHO daytime guideline value for hospital patient-care environments (35 dB[A]), the EN ISO 16032 reference value for general indoor environments (40 dB[A]), and the Chinese environmental quality standard GB3096-2022 Class 1 reference value (45 dB[A]). Because these reference values differ in scope, intended application, measurement methodology, and averaging period, the comparisons were performed for descriptive purposes only and do not represent formal regulatory compliance assessments (Figure 6). Across all 1,620 acoustic observations, environmental noise levels exceeded each of the selected reference values. Overall, 1,615 observations (99.7%) exceeded the WHO daytime guideline value of 35 dB(A), 1,586 observations (97.9%) exceeded the EN ISO 16032 reference value of 40 dB(A), and 1,513 observations (93.2%) exceeded the Chinese GB3096-2022 Class 1 reference value of 45 dB(A).

Figure 6: Percentage of Acoustic Observations Exceeding Selected International and National Reference Acoustic Criteria Across Hospital Functional Zones. Percentage of equivalent continuous A-weighted sound pressure level (LAeq) observations exceeding selected international and national reference acoustic criteria across the four standardized hospital functional zones. Bars represent the proportion of observations exceeding the reference thresholds of the World Health Organization (WHO) guidance value (>35 dB[A]), EN ISO 16032 (>40 dB[A]), and the Chinese National Standard GB3096-2022 (>45 dB[A]). Percentages were calculated as the number of observations exceeding each reference threshold divided by the total number of observations within each functional zone (n = 405). Because these reference documents differ in their intended application, measurement methodology, environmental setting, and averaging period, the comparisons are presented for descriptive purposes only and do not represent formal compliance assessments. LAeq, equivalent continuous A-weighted sound pressure level; WHO, World Health Organization. Please click here to view a larger version of this figure.
Marked differences were observed among the functional zones. Nurses’ stations demonstrated the highest proportion of observations exceeding each reference value, with 100.0% exceeding both the WHO and EN ISO reference values and 99.0% exceeding the Chinese reference value. Corridors likewise demonstrated high proportions of exceedance (100.0%, 99.3%, and 96.0%, respectively). In comparison, wards exhibited the lowest—but still substantial—proportions of observations exceeding the three reference values (99.3%, 95.3%, and 87.9%, respectively), whereas patient activity areas demonstrated intermediate proportions (99.5%, 97.0%, and 90.6%, respectively) These descriptive findings were consistent with the mixed-effects analyses, with nurses’ stations demonstrating the highest environmental noise levels and wards the lowest across the four functional zones. The high proportions of observations exceeding all three reference values indicate that elevated daytime environmental noise was observed throughout the hospital rather than being limited to individual departments or functional zones.
Sensitivity Analyses and Model Diagnostics
The robustness of the primary findings was evaluated using a fully adjusted sensitivity mixed-effects model that additionally included occupied beds, room area, room volume, staff count, alarm events, conversation events, and measurement period (Supplementary Table 2). Functional-zone estimates remained highly consistent with those observed in the prespecified parsimonious model. Nurses’ stations continued to demonstrate significantly higher LAeq than corridors, whereas patient activity areas and wards remained significantly quieter (all P < 0.001). Visitor count also remained independently associated with higher LAeq (β = 0.060, P = 0.023), whereas occupied beds, room area, room volume, staff count, alarm events, conversation events, and measurement period were not independently associated with environmental noise after simultaneous adjustment (Supplementary Table 2). These findings indicate that the principal findings were robust to additional adjustment for correlated operational characteristics. Model assumptions were evaluated using residual diagnostic procedures (Supplementary Figure 1). The normal quantile–quantile (Q–Q) plot demonstrated close agreement between the observed and theoretical residual quantiles, indicating no substantial departure from normality. The residual-versus-fitted plot showed random scatter around zero without discernible trends or evidence of heteroscedasticity, supporting the assumptions of linearity and constant residual variance. Consistent with these visual findings, the Shapiro–Wilk test did not demonstrate a statistically significant deviation from normality (W = 0.992, P = 0.084). Additional model-performance statistics are summarized in Supplementary Table 3. The residual standard deviation was 1.92 dB(A), and the department-level random-intercept standard deviation was 0.90 dB(A). The intraclass correlation coefficient was 0.172, indicating moderate clustering of observations within hospital departments. The mixed-effects model explained 53% of the variability through the fixed effects alone (marginal R2 = 0.53) and 61% when the department-level random effects were incorporated (conditional R2 = 0.61). Variance inflation factors were uniformly low (maximum VIF = 2.14), indicating no evidence of problematic multicollinearity among the retained predictors. Furthermore, all Cook’s distance values were <1.0 (maximum = 0.18), suggesting that no individual observation exerted undue influence on the fitted model. The exploratory correlation matrix (Supplementary Figure 2) demonstrated expected positive associations among several operational characteristics. Strong correlations were observed between room area and room volume (Spearman’s ρ = 0.94), bed capacity and occupied beds (ρ = 0.86), and bed capacity and room volume (ρ = 0.81). Moderate correlations were observed between visitor count and occupied beds (ρ = 0.71) and between staff count and conversation events (ρ = 0.62). These findings supported the prespecified variable-selection strategy for the parsimonious mixed-effects model, whereby highly correlated structural variables were excluded to minimize multicollinearity while preserving model interpretability.
Overall, the primary mixed-effects analyses, sensitivity analyses, residual diagnostics, and exploratory correlation assessment demonstrated internally consistent findings. Across all analyses, functional zone demonstrated the strongest independent association with environmental noise, whereas visitor count and departmental bed capacity also remained independently associated with LAeq. Most other recorded operational characteristics were not independently associated with environmental noise after mutual adjustment.
Data Availability:
The data supporting the findings of this study are available within the article and its Supplementary Materials. The complete raw acoustic measurement dataset, calibration records, operational-variable dataset, figure source data, and R statistical analysis code have been deposited in the Zenodo repository and are publicly available at https://doi.org/10.5281/zenodo.21404807. Supplementary Tables 1–3 and Supplementary Figures 1 and 2 provide additional processed datasets, sensitivity analyses, model diagnostics, and supporting analyses.
Supplementary Figure 1: Residual Diagnostic Plots for the Parsimonious Linear Mixed-Effects Model. Diagnostic plots evaluating the assumptions of the parsimonious linear mixed-effects regression model. (A) Normal quantile–quantile (Q–Q) plot of the conditional residuals demonstrating approximate normality of the residual distribution. (B) Plot of conditional residuals versus fitted values showing no substantial deviation from homoscedasticity or systematic pattern across the fitted range. Hospital department was included as a random intercept, and functional zone, visitor count, and bed capacity were included as fixed effects. Residual diagnostics support the adequacy of the model assumptions for the primary regression analysis. Diagnostic plots were generated using the DHARMa package in R. LAeq, equivalent continuous A-weighted sound pressure level.Please click here to download this file.
Supplementary Figure 2: Spearman Correlation Matrix of Operational Variables Included in the Exploratory Mixed-Effects Analyses. Spearman correlation matrix showing pairwise associations among departmental operational variables evaluated during exploratory mixed-effects regression analyses. Cell values represent Spearman’s rank correlation coefficients (ρ), and the color scale indicates the strength and direction of the correlation, with blue representing positive correlations and red representing negative correlations. Correlation analysis was performed to evaluate potential multicollinearity among candidate explanatory variables before model development. Variables demonstrating strong collinearity were considered during model selection to minimize multicollinearity in the parsimonious regression model. Operational variables included departmental bed capacity, occupied beds, room area, room volume, staff count, visitor count, alarm count, and conversation count. Spearman’s ρ, Spearman’s rank correlation coefficient.Please click here to download this file.
Supplementary Table 1: Secondary Acoustic Outcomes According to Hospital Functional Zone. Descriptive summary of the secondary acoustic outcomes measured across the four standardized hospital functional zones. Values are presented as mean ± standard deviation (SD) and represent pooled observations from all 15 inpatient departments. Equivalent continuous A-weighted sound pressure level (LAeq) was the primary acoustic outcome, whereas maximum A-weighted sound pressure level (LAFmax) and peak C-weighted sound pressure level (LCpeak) were recorded as secondary acoustic measures to characterize maximum and peak sound levels during routine daytime hospital operations. LAeq, equivalent continuous A-weighted sound pressure level; LAFmax, maximum A-weighted sound pressure level; LCpeak, peak C-weighted sound pressure level; SD, standard deviation.Please click here to download this file.
Supplementary Table 2: Fully Adjusted Sensitivity Linear Mixed-Effects Regression Model. Results of the fully adjusted sensitivity linear mixed-effects regression model evaluating associations between equivalent continuous A-weighted sound pressure level (LAeq) and all recorded operational characteristics. Hospital department was included as a random intercept to account for clustering of repeated observations. Functional zone, time period, departmental bed capacity, occupied beds, room area, room volume, staff count, visitor count, alarm count, and conversation count were included as fixed effects. Corridor served as the reference category for functional zone, and afternoon served as the reference category for time period. Regression coefficients (β), standard errors (SE), Wald 95% confidence intervals (CI), and two-sided P values are reported. This model was performed as a sensitivity analysis to evaluate the robustness of the parsimonious model presented in Table 4 after adjustment for all recorded operational characteristics. LAeq, equivalent continuous A-weighted sound pressure level; β, regression coefficient; SE, standard error; CI, confidence interval.Please click here to download this file.
Supplementary Table 3: Model Performance and Diagnostic Statistics for the Primary Linear Mixed-Effects Model. Summary of model performance and diagnostic statistics for the primary linear mixed-effects model evaluating equivalent continuous A-weighted sound pressure levels (LAeq) across hospital functional zones. Residual standard deviation (σ) represents the unexplained within-model variability, whereas the random-intercept standard deviation quantifies between-department variability. The intraclass correlation coefficient (ICC) estimates the proportion of total variance attributable to clustering within hospital departments. Marginal R2 represents the proportion of variance explained by the fixed effects alone, whereas conditional R2 represents the proportion of variance explained by both the fixed and random effects. Akaike information criterion (AIC), Bayesian information criterion (BIC), and log-likelihood are provided as measures of model fit. Variance inflation factors (VIFs) were all <5, indicating no evidence of problematic multicollinearity, and Cook’s distance values were <1.0 for all observations, indicating that no individual observation exerted undue influence on the model estimates. LAeq, equivalent continuous A-weighted sound pressure level; ICC, intraclass correlation coefficient; VIF, variance inflation factor.Please click here to download this file.