Caratteristiche dello studio e quadro di riferimento per le misurazioni
Le misurazioni del rumore ambientale sono state effettuate in 15 reparti ospedalieri di un grande ospedale universitario tra marzo e agosto 2025, utilizzando un protocollo di campionamento gerarchico standardizzato (Figura 1). All'interno di ciascun reparto, sono state valutate quattro zone funzionali predefinite (reparto, area di attività del paziente, corridoio e postazione infermieristica). In ogni zona funzionale sono stati stabiliti tre punti fissi di campionamento, e le misurazioni sono state effettuate durante tre periodi standardizzati di osservazione diurna (09:00–19:00). In ciascun punto e periodo di osservazione, sono state registrate tre misurazioni tecniche replicate consecutive mediante un fonometro di Classe 1 calibrato, generando dati acustici per LAeq, LAFmax e LCpeak. Variabili operative, tra cui capacità letti, letti occupati, numero di personale, numero di visitatori, eventi di allarme, eventi di conversazione, superficie delle stanze e volume delle stanze, sono state registrate contemporaneamente a ciascuna osservazione acustica. Il piano di campionamento prevedeva 1.620 osservazioni acustiche (15 reparti × 4 zone funzionali × 3 punti di campionamento × 3 periodi di osservazione × 3 replicate tecniche). Un'osservazione programmata nell'area di attività del paziente del Reparto di Terapia Intensiva II non è potuta essere completata a causa di restrizioni temporanee di accesso per cure d'emergenza ed è stata gestita mediante analisi dei casi completi senza imputazione dei dati, come specificato in anticipo nel piano di analisi statistica. Sono state effettuate ulteriori misurazioni in aree pubbliche interne ed esterne a scopo descrittivo comparativo, ma sono state escluse dalle analisi multilivello sui reparti (Figura 1). Le caratteristiche dei reparti partecipanti sono riassunte in Tabella 1. La capacità letti variava da 21 a 72 letti, con una media di letti occupati compresa tra 17,8 e 60,6. I conteggi medi del personale erano relativamente costanti tra i reparti (5,5–5,9 persone per osservazione), mentre il numero di visitatori medi era di circa cinque persone per osservazione. Le superfici delle stanze nei reparti variavano da 145,8 a 500,3 m2, e i volumi delle stanze da 466,5 a 1.601,0 m3, riflettendo una sostanziale variabilità nelle caratteristiche fisiche dei reparti partecipanti.

Figura 1: Progetto dello Studio e Protocollo Gerarchico per la Misurazione del Rumore Ambientale. Panoramica del progetto dello studio e del protocollo gerarchico per la misurazione del rumore ambientale. Il rumore ambientale è stato valutato in 15 reparti ospedalieri di un ospedale universitario mediante un framework di campionamento gerarchico multistadio. Le misurazioni sono state effettuate in quattro zone funzionali standardizzate (reparto, corridoio, postazione infermieristica e area di attività del paziente), in tre posizioni di campionamento predefinite (ingresso, centro e fondo), in tre periodi di osservazione diurni (mattina, mezzogiorno e pomeriggio) e con tre repliche tecniche consecutive. Il livello equivalente continuo di pressione sonora ponderato A (LAeq) è stato l'outcome primario, mentre il livello massimo di pressione sonora ponderato A (LAFmax) e il livello di picco di pressione sonora ponderato C (LCpeak) sono stati registrati come outcome acustici secondari. Le variabili operative sono state registrate contemporaneamente e tutte le misurazioni hanno subito procedure di calibrazione e di controllo qualità predefinite prima dell'analisi statistica. LAeq, livello equivalente continuo di pressione sonora ponderato A; LAFmax, livello massimo di pressione sonora ponderato A; LCpeak, livello di picco di pressione sonora ponderato C. Cliccare qui per visualizzare una versione ingrandita di questa figura.
| Dipartimento | Capacità Letti | Media Letti Occupati | Superficie Stanza (m²) | Volume Stanza (m³) | Media Personale | Media Visitatori | Media Allarmi | Media Conversazioni |
| Traumatologia Ortopedia | 40 | 33,8 | 280,1 | 896,3 | 5,5 | 4,8 | 2,4 | 6,5 |
| Centro Diagnosi e Trattamento Seno e Tiroide | 60 | 50,6 | 416,8 | 1333,8 | 5,7 | 4,8 | 2,4 | 6,2 |
| Centro Malattie Congenite del Cuore | 35 | 29,4 | 243,6 | 779,5 | 5,7 | 5 | 2,4 | 6,5 |
| Chirurgia Otorinolaringoiatria–Cervicologia | 50 | 42,2 | 348,1 | 1113,9 | 5,8 | 5 | 2,5 | 6,7 |
| Dipartimento Diagnosi e Trattamento Gastroenterologia I | 38 | 31,9 | 264,8 | 847,3 | 5,8 | 5 | 2,5 | 6,6 |
| Dipartimento Diagnosi e Trattamento Gastroenterologia II | 52 | 43,9 | 361,7 | 1157,4 | 5,7 | 4,9 | 2,4 | 6,5 |
| Chirurgia Epatobiliare e Pancreatica | 48 | 40,5 | 334,1 | 1069,2 | 5,8 | 5,1 | 2,5 | 6,8 |
| Dipartimento di Terapia Intensiva II | 21 | 17,8 | 145,8 | 466,5 | 5,7 | 4,8 | 2,5 | 6,5 |
| Dipartimento Articolazione del Ginocchio II | 45 | 37,8 | 312,7 | 1000,6 | 5,8 | 4,9 | 2,5 | 6,5 |
| Dipartimento Malattia di Parkinson | 55 | 46,2 | 382,2 | 1223 | 5,9 | 5,1 | 2,5 | 6,8 |
| Dipartimento Medicina Psicosomatica I | 68 | 57,1 | 472,6 | 1512,3 | 5,8 | 5 | 2,5 | 6,7 |
| Dipartimento di Pediatria | 40 | 33,6 | 278,6 | 891,5 | 5,8 | 5,1 | 2,5 | 6,8 |
| Dipartimento di Urologia | 62 | 52,1 | 431,2 | 1380 | 5,8 | 5 | 2,5 | 6,6 |
| Dipartimento Vascolare e Neurologico I | 72 | 60,6 | 500,3 | 1601 | 5,9 | 5,1 | 2,5 | 6,9 |
| Dipartimento Vascolare e Neurologico II | 60 | 50,4 | 417,5 | 1336 | 5,7 | 4,8 | 2,4 | 6,3 |
Tabella 1: Caratteristiche dei reparti ospedalieri inclusi nello studio. Caratteristiche operative di base dei 15 reparti di degenza inclusi nel rilevamento del rumore ambientale. Le variabili a livello di reparto comprendono la capacità di posti letto, il numero medio di letti occupati, la superficie delle stanze, il volume delle stanze e il numero medio di personale, visitatori, allarmi e conversazioni registrati durante periodi standardizzati di osservazione diurna. Le variabili operative rappresentano valori medi ottenuti durante i periodi di osservazione predeterminati e sono state valutate come variabili esplorative nelle analisi di regressione a effetti misti. La superficie delle stanze e il volume delle stanze sono riportati rispettivamente in metri quadrati (m2) e metri cubi (m3).
Caratteristiche acustiche descrittive nelle diverse zone funzionali
Tra tutte le 1.620 osservazioni, le postazioni infermieristiche hanno mostrato i livelli di rumore ambientale più elevati, seguite in sequenza da corridoi, aree di attività dei pazienti e stanze comuni (Tabella 2; Figure 2 e 3). I valori medi di LAeq sono stati di 66,90 ± 2,34 dB(A) presso le postazioni infermieristiche, 63,97 ± 2,38 dB(A) nei corridoi, 61,38 ± 2,17 dB(A) nelle aree di attività dei pazienti e 60,38 ± 2,22 dB(A) nelle stanze comuni. Gradienti spaziali simili sono stati osservati per LAFmax e LCpeak, indicando differenze costanti in tutte e tre le metriche acustiche. Sono riportati in Tabella supplementare 1 dettagliati riassunti descrittivi degli esiti acustici secondari (LAFmax e LCpeak).
| Zona Funzionale | Osservazioni (n) | LAeq, dB(A) Media ± DS | LAFmax, dB(A) Media ± DS | LCpeak, dB(C) Media ± DS | Mediana LAeq (IQR), dB(A) |
| Reparto | 405 | 60.38 ± 2.22 | 67.95 ± 2.37 | 79.21 ± 3.31 | 60.4 (58.9–61.9) |
| Corridoio | 405 | 63.97 ± 2.38 | 71.46 ± 2.59 | 82.76 ± 3.56 | 64.0 (62.5–65.5) |
| Posto di Guardia Infermieri | 405 | 66.90 ± 2.34 | 74.43 ± 2.56 | 85.72 ± 3.56 | 66.7 (65.4–68.3) |
| Area di Attività del Paziente | 405 | 61.38 ± 2.17 | 69.00 ± 2.42 | 80.20 ± 3.37 | 61.5 (60.0–62.7) |
Tabella 2: Caratteristiche acustiche descrittive per zona funzionale ospedaliera. Caratteristiche acustiche descrittive per le quattro zone funzionali ospedaliere standardizzate. I valori sono espressi come media ± deviazione standard (DS), salvo diversamente indicato. Sono riportati la mediana e il range interquartile (IQR) per il livello equivalente continuo di pressione sonora ponderato A (LAeq). Ogni zona funzionale ha contribuito con 405 osservazioni al dataset gerarchico (totale n = 1.620). LAeq, livello equivalente continuo di pressione sonora ponderato A; LAFmax, livello massimo di pressione sonora ponderato A; LCpeak, livello di picco di pressione sonora ponderato C; IQR, range interquartile.

Figura 2: Mappa Termica dei Livelli Equivalente Continui di Pressione Sonora Ponderati A Nei Reparti Ospedalieri e Nelle Zone Funzionali. Mappa termica che mostra i livelli equivalenti continui di pressione sonora ponderati A (LAeq) per ciascun reparto ospedaliero e zona funzionale. I valori delle celle rappresentano la media di LAeq in decibel ponderati A [dB(A)], e la scala cromatica indica un aumento dei livelli sonori da valori più bassi a valori più alti. Le zone funzionali includono reparti, corridoi, postazioni infermieristiche e aree di attività dei pazienti. CHD, malformazione cardiaca congenita; ICU, unità di terapia intensiva; LAeq, livello equivalente continuo di pressione sonora ponderato A. Cliccare qui per visualizzare una versione ingrandita di questa figura.

Figura 3: Distribuzione dei livelli equivalenti continui di pressione sonora ponderati A nelle zone funzionali ospedaliere standardizzate. Boxplot che mostrano la distribuzione dei livelli equivalenti continui di pressione sonora ponderati A (LAeq) in quattro zone funzionali ospedaliere standardizzate. I box rappresentano l'intervallo interquartile (25°–75° percentile), la linea orizzontale all'interno di ciascun box indica la mediana e il simbolo × indica la media. I baffi si estendono fino a 1,5 volte l'intervallo interquartile (IQR) e le osservazioni al di là dei baffi sono visualizzate come valori anomali. Ogni zona funzionale comprende 405 osservazioni (totale n = 1.620). Le differenze tra le zone funzionali sono state valutate utilizzando il modello primario a effetti misti lineari con confronti a coppie aggiustati secondo Tukey. LAeq, livello equivalente continuo di pressione sonora ponderato A; IQR, intervallo interquartile. Cliccare qui per visualizzare una versione ingrandita di questa figura.
La mappa termica presentata nella Figura 2 mostra un modello spaziale costante nei 15 reparti partecipanti. Le postazioni infermieristiche hanno mostrato costantemente i valori medi più elevati di LAeq, indipendentemente dalla specialità clinica, mentre le aree dei reparti hanno registrato i livelli più bassi di pressione sonora continua equivalente ponderata A. I reparti con maggiore capacità di posti letto, tra cui il Reparto di Vascolare e Neurologia I, il Reparto di Medicina Psicosomatica I e il Reparto di Urologia, hanno generalmente mostrato valori medi di LAeq più elevati nelle diverse zone funzionali rispetto ai reparti con minor numero di posti letto, come il Reparto di Terapia Intensiva II e il Centro per le Malattie Congenite del Cuore. Sebbene sia stata osservata una modesta variabilità tra i reparti, il posizionamento relativo delle zone funzionali è rimasto costante. I valori medi di LAeq presso le postazioni infermieristiche sono variati da 64,1 dB(A) nel Reparto di Terapia Intensiva II a 68,9 dB(A) nel Reparto di Vascolare e Neurologia I, mentre le misurazioni nei reparti sono comprese tra 57,9 e 61,6 dB(A). Le aree per l’attività dei pazienti e i corridoi hanno mostrato valori intermedi in quasi tutti i reparti (Figura 2). La Figura 3 illustra la distribuzione delle misurazioni di LAeq all’interno di ciascuna zona funzionale. I valori medi di LAeq sono risultati più elevati presso le postazioni infermieristiche (66,7 dB[A]), seguiti dai corridoi (64,0 dB[A]), dalle aree per l’attività dei pazienti (61,5 dB[A]) e dai reparti (60,4 dB[A]). Sebbene siano state rilevate occasionalmente misurazioni con valori elevati in tutte le zone funzionali, gli intervalli interquartili erano relativamente ristretti, indicando una variabilità limitata all’interno di ciascun ambiente funzionale. I boxplot mostrano una chiara separazione nella distribuzione dei valori di LAeq tra le quattro zone funzionali (Figura 3).
Confronto a Effetti Misti delle Zone Funzionali
Per tenere conto della struttura gerarchica delle misurazioni acustiche ripetute all'interno dei reparti ospedalieri, l'analisi principale ha utilizzato un modello lineare a effetti misti con il reparto ospedaliero specificato come intercetta casuale e la zona funzionale come effetto fisso (Tabella 3). Questo approccio modellistico ha considerato il raggruppamento all'interno dei singoli reparti, stimando al contempo le differenze aggiustate nei livelli di LAeq tra le quattro zone funzionali predefinite.
| Effetti Fissi (Esito: LAeq) |
| Effetto Fisso | β | SE | IC 95% | Valore P |
| Intercetta (Corridoio) | 63,972 | 0,266 | 63,451 a 64,493 | <0,001 |
| Posto Infermieri | 2,931 | 0,146 | 2,644 a 3,218 | <0,001 |
| Area Attività Paziente | −2,594 | 0,146 | −2,881 a −2,307 | <0,001 |
| Reparto | −3,589 | 0,146 | −3,876 a −3,302 | <0,001 |
| Medie Marginali Stimate |
| Zona Funzionale | Media Stimata LAeq, dB(A) |
| Posto Infermieri | 66,9 |
| Corridoio | 63,97 |
| Area Attività Paziente | 61,38 |
| Reparto | 60,38 |
| Confronti a Coppie con Aggiustamento di Tukey |
| Confronto | Differenza Media, dB(A) | Valore P Aggiustato |
| Posto Infermieri vs. Corridoio | 2,93 | <0,001 |
| Posto Infermieri vs. Area Attività Paziente | 5,52 | <0,001 |
| Posto Infermieri vs. Reparto | 6,52 | <0,001 |
| Corridoio vs. Area Attività Paziente | 2,59 | <0,001 |
| Corridoio vs. Reparto | 3,59 | <0,001 |
| Area Attività Paziente vs. Reparto | 0,99 | <0,001 |
| Prestazioni del Modello |
| Statistica | Valore |
| Coefficiente di Correlazione Intraclasse (ICC) | 0,172 |
| R² Marginale | 0,53 |
| R² Condizionale | 0,61 |
Tabella 3: Modello primario a effetti misti lineari per il confronto del rumore ambientale tra le diverse zone funzionali dell'ospedale. Risultati del modello primario a effetti misti lineari volto a valutare le differenze nei livelli equivalenti continui di pressione sonora ponderata A (LAeq) tra le quattro zone funzionali standardizzate dell'ospedale. Il reparto ospedaliero è stato definito come intercetta casuale per tenere conto del raggruppamento delle osservazioni ripetute, mentre la zona funzionale è stata inclusa come effetto fisso, con il corridoio come categoria di riferimento. Sono riportati i coefficienti di regressione (β), gli errori standard (SE), gli intervalli di confidenza al 95% calcolati con il metodo di Wald (IC) e i valori P bilaterali per gli effetti fissi. Le medie marginali stimate sono state confrontate mediante confronti multipli a coppie corretti con il metodo di Tukey. Le prestazioni del modello sono riassunte tramite il coefficiente di correlazione intraclasse (ICC) e i coefficienti marginali e condizionali di determinazione (R2), calcolati secondo il metodo di Nakagawa e Schielzeth. LAeq, livello equivalente continuo di pressione sonora ponderata A; β, coefficiente di regressione; SE, errore standard; IC, intervallo di confidenza; ICC, coefficiente di correlazione intraclasse.
L'analisi a effetti misti ha dimostrato un effetto complessivo altamente significativo della zona funzionale sui livelli di rumore ambientale (test del rapporto di verosimiglianza complessivo, P < 0,001). Utilizzando il corridoio come categoria di riferimento, le postazioni infermieristiche hanno mostrato valori di LAeq significativamente più elevati (β = 2,931 dB[A], SE = 0,146, IC 95% = 2,644–3,218, P < 0,001), mentre le aree di attività del paziente (β = −2,594 dB[A], SE = 0,146, IC 95% = −2,881 a −2,307, P < 0,001) e i reparti (β = −3,589 dB[A], SE = 0,146, IC 95% = −3,876 a −3,302, P < 0,001) hanno mostrato valori di LAeq significativamente più bassi (Tabella 3). Le medie marginali stimate derivate dal modello a effetti misti hanno evidenziato un chiaro schema gerarchico del rumore ambientale nelle diverse zone funzionali. La media aggiustata di LAeq è risultata più alta nelle postazioni infermieristiche (66,90 dB[A]), seguita dai corridoi (63,97 dB[A]), dalle aree di attività del paziente (61,38 dB[A]) e dai reparti (60,38 dB[A]). Tutti i confronti a coppie sono rimasti statisticamente significativi dopo la correzione di Tukey (tutti i P aggiustati < 0,001), indicando differenze statisticamente significative nei valori di LAeq tra tutte e quattro le zone funzionali. Le differenze medie variavano da 0,99 dB(A) tra le aree di attività del paziente e i reparti a 6,52 dB(A) tra le postazioni infermieristiche e i reparti. Figura 4 riassume le stime aggiustate dal modello. Le medie marginali stimate e i rispettivi intervalli di confidenza al 95% hanno mostrato un minimo sovrapposizione tra le zone funzionali, e il raggruppamento secondo Tukey ha assegnato ciascuna zona a un gruppo di significatività unico (A–D). Le postazioni infermieristiche sono state assegnate al gruppo A, con il LAeq aggiustato più alto, seguite dai corridoi (gruppo B), dalle aree di attività del paziente (gruppo C) e dai reparti (gruppo D). L'analisi delle componenti della varianza ha mostrato che circa il 17,2% della variabilità totale in LAeq è attribuibile alle differenze tra i reparti ospedalieri (coefficiente di correlazione intraclasse = 0,172), supportando l'inclusione di un intercetta casuale a livello di reparto. Gli effetti fissi hanno spiegato il 53% della variabilità del rumore ambientale (R2 marginale = 0,53), mentre il modello combinato con effetti fissi e casuali ha spiegato il 61% della variabilità totale (R2 condizionale = 0,61), indicando un'elevata capacità esplicativa del modello gerarchico (Tabella 3).

Figura 4: Medie Marginali Stimate dei Livelli Equivalenti Continui di Pressione Sonora Ponderata A nei Diversi Settori Funzionali dell'Ospedale. Medie marginali stimate (EMM) dei livelli equivalenti continui di pressione sonora ponderata A (LAeq) nei quattro settori funzionali ospedalieri, ottenuti dal modello lineare primario a effetti misti. I punti rappresentano le medie marginali stimate dal modello, mentre le barre d'errore orizzontali indicano gli intervalli di confidenza (IC) al 95%. Il reparto ospedaliero è stato incluso come effetto casuale e il settore funzionale come effetto fisso. I confronti a coppie sono stati effettuati mediante medie marginali stimate aggiustate con il metodo di Tukey; i settori funzionali contrassegnati con lettere diverse presentano differenze statisticamente significative (P < 0,05). Ogni settore funzionale ha contribuito con 405 osservazioni all'analisi. LAeq, livello equivalente continuo di pressione sonora ponderata A; EMM, media marginale stimata; IC, intervallo di confidenza. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Regressione parsimoniosa a effetti misti delle caratteristiche operative associate al LAeq
Per analizzare le caratteristiche operative associate al rumore ambientale riducendo al contempo la multicollinearità, è stato adattato un modello prespecificato parsimonioso a effetti misti, includendo come effetti fissi la zona funzionale, il numero di visitatori e la capacità letti del reparto, e mantenendo il reparto ospedaliero come intercetta casuale (Tabella 4). La selezione di queste variabili è stata guidata dalla matrice esplorativa di correlazione (Figura Supplementare 1), che ha evidenziato forti correlazioni tra diverse variabili strutturali, in particolare tra superficie della stanza e volume della stanza (rho di Spearman = 0,94). Di conseguenza, le variabili strutturali altamente correlate sono state escluse dal modello parsimonioso per migliorarne l'interpretabilità e la stabilità. Dopo aggiustamento per i covarianti operativi, la zona funzionale è rimasta il fattore indipendente più forte associato al rumore ambientale. Rispetto ai corridoi, le postazioni infermieristiche erano associate a un aumento aggiustato di 2,933 dB(A) (SE = 0,146, IC 95% = 2,647–3,220, P < 0,001), mentre le aree di attività del paziente (β = −2,602 dB[A], SE = 0,146, IC 95% = −2,888 a −2,315, P < 0,001) e i reparti (β = −3,606 dB[A], SE = 0,146, IC 95% = −3,892 a −3,319, P < 0,001) risultavano significativamente più silenziosi rispetto ai corridoi.
| Predittore | Coefficiente β | SE | IC 95% | Valore P |
| Intersezione (Corridoio) | 63,97 | 0,27 | 63,45 a 64,49 | <0,001 |
| Posto degli infermieri (rispetto al Corridoio) | 2,933 | 0,146 | 2,647 a 3,220 | <0,001 |
| Area attività paziente (rispetto al Corridoio) | −2,602 | 0,146 | −2,888 a −2,315 | <0,001 |
| Reparto (rispetto al Corridoio) | −3,606 | 0,146 | −3,892 a −3,319 | <0,001 |
| Numero visitatori (per visitatore aggiuntivo) | 0,06 | 0,026 | 0,008 a 0,112 | 0,024 |
| Capacità letti (per letto aggiuntivo) | 0,064 | 0,005 | 0,055 a 0,073 | <0,001 |
| Modello: LAeq ~ Zona Funzionale + Numero Visitatori + Capacità Letti + (1 | Dipartimento) |
Tabella 4: Regressione lineare a effetti misti parsimoniosa delle caratteristiche operative associate al rumore ambientale. Risultati del modello di regressione lineare a effetti misti parsimonioso predeterminato, utilizzato per valutare le associazioni tra il livello equivalente continuo di pressione sonora ponderato in scala A (LAeq) e alcune caratteristiche operative selezionate. Il modello includeva come effetti fissi la zona funzionale, il numero di visitatori e la capacità letti del reparto, e come intercetta casuale il reparto ospedaliero, per tenere conto del raggruppamento delle osservazioni ripetute (modello: LAeq ~ Zona Funzionale + Numero Visitatori + Capacità Letti + (1 | Reparto)). La categoria di riferimento per la zona funzionale era il corridoio, mentre il numero di visitatori e la capacità letti erano modellati come variabili continue. Sono riportati i coefficienti di regressione (β), gli errori standard (SE), gli intervalli di confidenza al 95% di Wald (IC) e i valori P bilaterali. I coefficienti di regressione per il numero di visitatori e la capacità letti rappresentano la variazione stimata del LAeq associata a ciascun visitatore aggiuntivo e a ciascun letto ospedaliero aggiuntivo, rispettivamente. L'intercetta rappresenta la media stimata del LAeq per la categoria di riferimento (corridoio) quando i predittori continui sono pari a zero. Questa tabella presenta il modello parsimonioso predeterminato sviluppato per minimizzare l'eccesso di adattamento del modello; il modello di sensibilità completamente aggiustato, che include tutte le caratteristiche operative registrate, è riportato nella Tabella S2 Supplementare. LAeq, livello equivalente continuo di pressione sonora ponderato in scala A; β, coefficiente di regressione; SE, errore standard; IC, intervallo di confidenza.
Tra i predittori operativi continui, il numero di visitatori ha dimostrato un'associazione positiva indipendente con LAeq. Ogni visitatore aggiuntivo presente durante la misurazione era associato a un aumento stimato di 0,060 dB(A) (SE = 0,026, IC 95% = 0,008–0,112, P = 0,024). Analogamente, la capacità di posti letto del reparto è rimasta indipendentemente associata al rumore ambientale, con ogni letto ospedaliero aggiuntivo associato a un aumento stimato di 0,064 dB(A) in LAeq (SE = 0,005, IC 95% = 0,055–0,073, P < 0,001).
Figura 5 illustra le previsioni adattate generate dal modello a effetti misti parsimonioso. L'LAeq previsto è aumentato approssimativamente in modo lineare nell'intervallo osservato di capacità di posti letto del reparto, e un numero maggiore di visitatori ha spostato costantemente verso l'alto i livelli sonori previsti lungo l'intero intervallo di capacità. Le rette di regressione adattate parallele indicano effetti additivi del numero di visitatori e della capacità di posti letto all'interno del modello adattato, dopo la correzione per zona funzionale e raggruppamento a livello di reparto. Le previsioni basate sul modello hanno dimostrato che, per un reparto con circa 20 posti letto, l'LAeq previsto variava da 61,6 dB(A) con un numero ridotto di visitatori (tre visitatori) a 62,4 dB(A) con un numero maggiore di visitatori (sette visitatori). All'estremo superiore dell'intervallo osservato di capacità di posti letto (72 posti letto), l'LAeq previsto corrispondente è aumentato rispettivamente a 64,5 dB(A) e 65,3 dB(A), evidenziando il contributo combinato della capacità di posti letto del reparto e del numero di visitatori all'interno del modello adattato.

Figura 5: Previsioni basate sul modello dei livelli equivalenti continui di pressione sonora ponderati A in funzione della capacità di letti del reparto. Livelli equivalenti continui di pressione sonora ponderati A (LAeq) previsti nell'intervallo osservato di capacità di letti del reparto, ottenuti dal modello semplificato di regressione lineare a effetti misti. Le linee colorate rappresentano i valori di LAeq previsti dal modello a un numero basso (3 visitatori), medio (5 visitatori) e alto (7 visitatori) di visitatori, mantenendo costanti tutti gli altri termini del modello. Le aree ombreggiate indicano gli intervalli di confidenza al 95% per i valori previsti, mentre i punti rappresentano le misurazioni osservate. Le previsioni sono mostrate soltanto all'interno dell'intervallo osservato di capacità di letti del reparto incluso nello studio. LAeq, livello equivalente continuo di pressione sonora ponderato A; CI, intervallo di confidenza. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Le analisi di sensibilità condotte utilizzando il modello a effetti misti completamente aggiustato (Tabella Supplementare 2), che includeva inoltre letti occupati, superficie della stanza, volume della stanza, numero di personale, eventi di allarme, eventi di conversazione e periodo di misurazione, hanno prodotto stime altamente coerenti per la zona funzionale e il numero di visitatori. Le postazioni infermieristiche sono rimaste la zona funzionale più rumorosa, e il numero di visitatori è risultato indipendentemente associato a un maggiore LAeq (β = 0,060, P = 0,023), mentre letti occupati, superficie della stanza, volume della stanza, numero di personale, eventi di allarme, eventi di conversazione e periodo di misurazione non sono risultati indipendentemente associati al rumore ambientale dopo l’aggiustamento reciproco. Questi risultati confermano la robustezza del modello parsimonioso e indicano che la zona funzionale e il numero di visitatori sono state le variabili principali associate al rumore ambientale in questo studio.
Confronto descrittivo con criteri acustici di riferimento
Per contestualizzare i livelli di rumore ambientale osservati, i valori misurati di LAeq sono stati confrontati descrittivamente con tre criteri acustici di riferimento ampiamente citati: il valore guida dell'OMS per le ore diurne negli ambienti ospedalieri di cura dei pazienti (35 dB[A]), il valore di riferimento della norma EN ISO 16032 per gli ambienti interni generici (40 dB[A]) e il valore di riferimento della classe 1 dello standard cinese sulla qualità ambientale GB3096-2022 (45 dB[A]). Poiché questi valori di riferimento differiscono per ambito, applicazione prevista, metodologia di misurazione e periodo di media, i confronti sono stati effettuati esclusivamente a scopo descrittivo e non rappresentano valutazioni formali di conformità normativa (Figura 6). Su un totale di 1.620 osservazioni acustiche, i livelli di rumore ambientale hanno superato ciascuno dei valori di riferimento selezionati. Nel complesso, 1.615 osservazioni (99,7%) hanno superato il valore guida dell'OMS di 35 dB(A), 1.586 osservazioni (97,9%) hanno superato il valore di riferimento EN ISO 16032 di 40 dB(A) e 1.513 osservazioni (93,2%) hanno superato il valore di riferimento cinese GB3096-2022 classe 1 di 45 dB(A).

Figura 6: Percentuale di Osservazioni Acustiche che Superano i Criteri Acustici di Riferimento Internazionali e Nazionali Selezionati nelle Zone Funzionali dell'Ospedale. Percentuale di osservazioni del livello equivalente continuo di pressione sonora ponderato A (LAeq) che superano i criteri acustici di riferimento internazionali e nazionali selezionati nelle quattro zone funzionali ospedaliere standardizzate. Le barre rappresentano la proporzione di osservazioni che superano le soglie di riferimento delle indicazioni dell'Organizzazione Mondiale della Sanità (OMS) (>35 dB[A]), della norma EN ISO 16032 (>40 dB[A]) e della norma nazionale cinese GB3096-2022 (>45 dB[A]). Le percentuali sono state calcolate come numero di osservazioni superiori a ciascuna soglia di riferimento diviso per il numero totale di osservazioni all'interno di ciascuna zona funzionale (n = 405). Poiché questi documenti di riferimento differiscono per applicazione prevista, metodologia di misurazione, contesto ambientale e periodo di mediazione, i confronti sono presentati esclusivamente a scopo descrittivo e non rappresentano valutazioni formali di conformità. LAeq, livello equivalente continuo di pressione sonora ponderato A; OMS, Organizzazione Mondiale della Sanità. Cliccare qui per visualizzare una versione ingrandita di questa figura.
Sono state osservate differenze significative tra le zone funzionali. Le postazioni infermieristiche hanno mostrato la percentuale più alta di osservazioni superiori a ciascun valore di riferimento, con il 100,0% che superava sia i valori di riferimento dell'OMS che quelli della norma EN ISO, e il 99,0% che superava il valore di riferimento cinese. Anche i corridoi hanno mostrato proporzioni elevate di superamento (rispettivamente 100,0%, 99,3% e 96,0%). In confronto, i reparti hanno mostrato le percentuali più basse — sebbene ancora sostanziali — di osservazioni superiori ai tre valori di riferimento (rispettivamente 99,3%, 95,3% e 87,9%), mentre le aree di attività dei pazienti hanno mostrato proporzioni intermedie (rispettivamente 99,5%, 97,0% e 90,6%). Questi risultati descrittivi sono coerenti con le analisi a effetti misti, in cui le postazioni infermieristiche presentavano i livelli più elevati di rumore ambientale e i reparti i livelli più bassi tra le quattro zone funzionali. Le elevate proporzioni di osservazioni che superano tutti e tre i valori di riferimento indicano che livelli elevati di rumore ambientale diurno sono stati rilevati in tutta l'ospedale, piuttosto che essere limitati a singoli reparti o zone funzionali.
Analisi di Sensibilità e Diagnostica del Modello
La robustezza dei risultati principali è stata valutata utilizzando un modello misto a effetti misti completamente aggiustato, che includeva inoltre letti occupati, superficie della stanza, volume della stanza, numero di personale, eventi di allarme, eventi di conversazione e periodo di misurazione (Tabella Supplementare 2). Le stime delle zone funzionali sono rimaste molto coerenti con quelle osservate nel modello parsimonioso prefissato. Le postazioni degli infermieri hanno continuato a mostrare livelli di LAeq significativamente più elevati rispetto ai corridoi, mentre le aree di attività dei pazienti e i reparti sono rimasti significativamente più silenziosi (tutti P < 0,001). Il numero di visitatori è rimasto indipendentemente associato a un LAeq più alto (β = 0,060, P = 0,023), mentre letti occupati, superficie della stanza, volume della stanza, numero di personale, eventi di allarme, eventi di conversazione e periodo di misurazione non sono risultati indipendentemente associati al rumore ambientale dopo aggiustamento simultaneo (Tabella Supplementare 2). Questi risultati indicano che i principali risultati erano robusti anche dopo ulteriore aggiustamento per caratteristiche operative correlate. Le assunzioni del modello sono state valutate mediante procedure diagnostiche sui residui (Figura Supplementare 1). Il grafico quantile-normale (Q–Q) ha mostrato un accordo molto stretto tra i quantili osservati e quelli teorici, indicando alcuna deviazione sostanziale dalla normalità. Il grafico dei residui rispetto ai valori previsti ha mostrato una distribuzione casuale intorno allo zero, senza tendenze evidenti o segni di eteroschedasticità, sostenendo le assunzioni di linearità e varianza costante dei residui. In linea con questi risultati visivi, il test di Shapiro–Wilk non ha mostrato una deviazione statisticamente significativa dalla normalità (W = 0,992, P = 0,084). Ulteriori statistiche sulle prestazioni del modello sono riassunte nella Tabella Supplementare 3. La deviazione standard dei residui era di 1,92 dB(A), e la deviazione standard dell'intercetta casuale a livello di reparto era di 0,90 dB(A). Il coefficiente di correlazione intraclasse era 0,172, indicando un raggruppamento moderato delle osservazioni all'interno dei reparti ospedalieri. Il modello a effetti misti ha spiegato il 53% della variabilità attraverso gli effetti fissi soltanto (R2 marginale = 0,53) e il 61% quando sono stati inclusi gli effetti casuali a livello di reparto (R2 condizionale = 0,61). I fattori di inflazione della varianza erano uniformemente bassi (VIF massimo = 2,14), indicando assenza di multicollinearità problematica tra i predittori mantenuti. Inoltre, tutti i valori della distanza di Cook erano <1,0 (massimo = 0,18), suggerendo che nessuna singola osservazione esercitava un'influenza eccessiva sul modello stimato. La matrice esplorativa delle correlazioni (Figura Supplementare 2) ha mostrato associazioni positive attese tra diverse caratteristiche operative. Correlazioni forti sono state osservate tra superficie della stanza e volume della stanza (ρ di Spearman = 0,94), capacità letti e letti occupati (ρ = 0,86), e capacità letti e volume della stanza (ρ = 0,81). Correlazioni moderate sono state osservate tra numero di visitatori e letti occupati (ρ = 0,71) e tra numero di personale ed eventi di conversazione (ρ = 0,62). Questi risultati hanno sostenuto la strategia prespecificata di selezione delle variabili per il modello parsimonioso a effetti misti, secondo cui variabili strutturali altamente correlate sono state escluse per minimizzare la multicollinearità preservando l'interpretabilità del modello.
Nell'insieme, le analisi principali a effetti misti, le analisi di sensibilità, la diagnostica dei residui e la valutazione esplorativa delle correlazioni hanno dimostrato risultati internamente coerenti. In tutte le analisi, la zona funzionale ha mostrato la più forte associazione indipendente con il rumore ambientale, mentre il numero di visitatori e la capacità massima di posti letto del reparto sono risultati anch'essi indipendentemente associati all'LAeq. La maggior parte delle altre caratteristiche operative registrate non ha mostrato un'associazione indipendente con il rumore ambientale dopo l'aggiustamento reciproco.
Disponibilità dei dati:
I dati a supporto dei risultati di questo studio sono disponibili all'interno dell'articolo e dei suoi Materiali Supplementari. L'insieme completo dei dati grezzi delle misurazioni acustiche, i record di calibrazione, l'insieme di dati sulle variabili operative, i dati sorgente delle figure e il codice per l'analisi statistica in R sono stati depositati nel repository Zenodo e sono pubblicamente disponibili all'indirizzo https://doi.org/10.5281/zenodo.21404807. Tabelle Supplementari 1–3 e Figure Supplementari 1 e 2 forniscono ulteriori insiemi di dati elaborati, analisi di sensibilità, diagnosi del modello e analisi di supporto.
Figura supplementare 1: Grafici diagnostici dei residui per il modello lineare a effetti misti parsimonioso. Grafici diagnostici per la valutazione delle assunzioni del modello di regressione lineare a effetti misti parsimonioso. (A) Grafico quantile–quantile normale (Q–Q) dei residui condizionali che mostra l'approssimata normalità della distribuzione dei residui. (B) Grafico dei residui condizionali rispetto ai valori stimati, che non evidenzia deviazioni significative dall'omoschedasticità né schemi sistematici lungo l'intervallo dei valori stimati. Il reparto ospedaliero è stato incluso come intercetta casuale, mentre la zona funzionale, il numero di visitatori e la capacità letti sono stati inclusi come effetti fissi. I controlli diagnostici sui residui confermano l'adeguatezza delle assunzioni del modello per l'analisi di regressione principale. I grafici diagnostici sono stati generati utilizzando il pacchetto DHARMa in R. LAeq, livello equivalente continuo di pressione sonora pesata A.Cliccare qui per scaricare questo file.
Figura supplementare 2: Matrice di Correlazione di Spearman delle Variabili Operative Incluse nelle Analisi Esplorative a Effetti Misti. Matrice di correlazione di Spearman che mostra le associazioni a coppie tra le variabili operative dipartimentali valutate durante le analisi esplorative di regressione a effetti misti. I valori delle celle rappresentano i coefficienti di correlazione per ranghi di Spearman (ρ), e la scala cromatica indica l'intensità e la direzione della correlazione, con il blu che indica correlazioni positive e il rosso correlazioni negative. L'analisi di correlazione è stata effettuata per valutare la possibile multicollinearità tra le variabili esplicative candidate prima dello sviluppo del modello. Le variabili che mostravano una forte collinearità sono state prese in considerazione durante la selezione del modello al fine di minimizzare la multicollinearità nel modello di regressione parsimonioso. Le variabili operative incluse erano la capacità massima di letti del reparto, i letti occupati, la superficie della stanza, il volume della stanza, il numero di personale, il numero di visitatori, il numero di allarmi e il numero di conversazioni. Spearman’s ρ, coefficiente di correlazione per ranghi di Spearman.Cliccare qui per scaricare questo file.
Tabella Supplementare 1: Esiti Acustici Secondari in Base alla Zona Funzionale dell'Ospedale. Riassunto descrittivo degli esiti acustici secondari misurati nelle quattro zone funzionali ospedaliere standardizzate. I valori sono espressi come media ± deviazione standard (SD) e rappresentano osservazioni aggregate provenienti da tutti e 15 i reparti di degenza. Il livello equivalente continuo di pressione sonora ponderato A (LAeq) è stato l'esito acustico primario, mentre il livello massimo di pressione sonora ponderato A (LAFmax) e il livello di picco di pressione sonora ponderato C (LCpeak) sono stati registrati come misure acustiche secondarie per caratterizzare i livelli sonori massimi e di picco durante le normali operazioni diurne ospedaliere. LAeq, livello equivalente continuo di pressione sonora ponderato A; LAFmax, livello massimo di pressione sonora ponderato A; LCpeak, livello di picco di pressione sonora ponderato C; SD, deviazione standard.Cliccare qui per scaricare il file.
Tabella Supplementare 2: Modello di regressione lineare a effetti misti completamente aggiustato per l'analisi di sensibilità. Risultati del modello di regressione lineare a effetti misti completamente aggiustato per l'analisi di sensibilità, volto a valutare le associazioni tra il livello equivalente continuo di pressione sonora ponderata A (LAeq) e tutte le caratteristiche operative registrate. Il reparto ospedaliero è stato incluso come intercetta casuale per tenere conto del raggruppamento delle osservazioni ripetute. Zona funzionale, periodo temporale, capacità letti del reparto, letti occupati, superficie della stanza, volume della stanza, numero di operatori, numero di visitatori, numero di allarmi e numero di conversazioni sono stati inclusi come effetti fissi. Il corridoio è stato utilizzato come categoria di riferimento per la zona funzionale, e il pomeriggio come categoria di riferimento per il periodo temporale. Sono riportati i coefficienti di regressione (β), gli errori standard (SE), gli intervalli di confidenza al 95% calcolati con il test di Wald (IC) e i valori P a due code. Questo modello è stato eseguito come analisi di sensibilità per valutare la robustezza del modello semplificato presentato nella Tabella 4, dopo aggiustamento per tutte le caratteristiche operative registrate. LAeq, livello equivalente continuo di pressione sonora ponderata A; β, coefficiente di regressione; SE, errore standard; IC, intervallo di confidenza.Cliccare qui per scaricare il file.
Tabella Supplementare 3: Prestazioni del Modello e Statistiche Diagnostiche per il Modello Primario a Effetti Misti Lineari. Riassunto delle prestazioni del modello e delle statistiche diagnostiche per il modello primario a effetti misti lineari che valuta i livelli equivalenti continui di pressione sonora ponderati A (LAeq) nelle diverse zone funzionali dell'ospedale. La deviazione standard dei residui (σ) rappresenta la variabilità non spiegata all'interno del modello, mentre la deviazione standard dell'intercetta casuale quantifica la variabilità tra i reparti. Il coefficiente di correlazione intraclasse (ICC) stima la proporzione di varianza totale attribuibile al raggruppamento all'interno dei reparti ospedalieri. L'R2 marginale rappresenta la proporzione di varianza spiegata soltanto dagli effetti fissi, mentre l'R2 condizionale rappresenta la proporzione di varianza spiegata sia dagli effetti fissi che da quelli casuali. Il criterio di informazione di Akaike (AIC), il criterio di informazione bayesiano (BIC) e la log-verosimiglianza sono forniti come misure dell'adattamento del modello. Tutti i fattori di inflazione della varianza (VIF) erano <5, indicando l'assenza di multicollinearità problematica; inoltre, i valori della distanza di Cook erano <1,0 per tutte le osservazioni, indicando che nessuna osservazione singola ha esercitato un'influenza eccessiva sulle stime del modello. LAeq, livello equivalente continuo di pressione sonora ponderato A; ICC, coefficiente di correlazione intraclasse; VIF, fattore di inflazione della varianza.Cliccare qui per scaricare il file.