Charakterystyka badania i ramy pomiarowe
Pomiary hałasu środowiskowego przeprowadzono w 15 oddziałach stacjonarnych dużego trzeciorzędnego szpitala klinicznego między marcem a sierpniem 2025 roku, stosując zstandaryzowany hierarchiczny protokół próbkowania (Rysunek 1). W każdym oddziale oceniono cztery zdefiniowane strefy funkcjonalne (salę chorych, obszar aktywności pacjentów, korytarz oraz dyżurkę pielęgniarską). W każdej strefie funkcjonalnej wyznaczono trzy stałe punkty próbkowania, a pomiary wykonywano podczas trzech zstandaryzowanych dziennych okresów obserwacyjnych (09:0–19:00). W każdym punkcie i okresie obserwacyjnym zapisano trzy kolejne techniczne powtórzenia pomiarów przy użyciu skalibrowanego miernika poziomu dźwięku klasy 1, generując dane akustyczne dla LAeq, LAFmax oraz LCpeak. Równolegle z każdą obserwacją akustyczną rejestrowano zmienne operacyjne, w tym liczbę łóżek, liczbę łóżek zajętych, liczbę personelu, liczbę odwiedzających, zdarzenia alarmowe, zdarzenia konwersacyjne, powierzchnię pomieszczenia oraz objętość pomieszczenia. Zaplanowany schemat próbkowania obejmował 1 620 obserwacji akustycznych (15 oddziałów × 4 strefy funkcjonalne × 3 punkty próbkowania × 3 okresy obserwacyjne × 3 powtórzenia techniczne). Jedna zaplanowana obserwacja w obszarze aktywności pacjentów Oddziału Intensywnej Terapii II nie mogła zostać zrealizowana z powodu tymczasowych ograniczeń dostępu do opieki ratunkowej i została potraktowana zgodnie z analizą przypadków kompletnych bez imputacji danych, zgodnie z wcześniej określonym planem analizy statystycznej. Dodatkowe pomiary w pomieszczeniach wewnętrznych oraz w ogólnodostępnych obszarach zewnętrznych wykonano w celu opisowego porównania kontekstowego, jednak zostały one wyłączone z analizy efektów mieszanych dla oddziałów (Rysunek 1). Charakterystyka uczestniczących oddziałów została podsumowana w Tabeli 1. Liczba łóżek wynosiła od 21 do 72, a odpowiadająca im średnia liczba łóżek zajętych wahała się od 17,8 do 60,6. Średnia liczba personelu była stosunkowo spójna we wszystkich oddziałach (5,5–5,9 pracowników na obserwację), podczas gdy liczba odwiedzających wynosiła średnio około pięciu osób na obserwację. Powierzchnie pomieszczeń na oddziałach wynosiły od 145,8 do 50,3 m2, a objętości od 466,5 do 1 601,0 m3, co odzwierciedla znaczne różnice w charakterystyce fizycznej badanych oddziałów.

Rycina 1: Projekt badania i hierarchiczny protokół pomiaru hałasu środowiskowego. Przegląd projektu badania i hierarchicznego protokołu pomiaru hałasu środowiskowego. Hałas środowiskowy oceniano w 15 oddziałach stacjonarnych trzeciostopniowego szpitala klinicznego z wykorzystaniem wielostopniowego hierarchicznego schematu próbkowania. Pomiary przeprowadzono w czterech ustandaryzowanych strefach funkcjonalnych (sala chorych, korytarz, punkt pielęgniarski oraz strefa aktywności pacjentów), w trzech zdefiniowanych lokalizacjach próbkowania (wejście, środek i koniec), w trzech dziennych okresach obserwacji (rano, w południe i po południu) oraz w trzech kolejnych powtórzeniach technicznych. Głównym punktem końcowym był równoważny ciągły poziom ciśnienia akustycznego z filtrem A (LAeq), natomiast jako drugorzędne parametry akustyczne rejestrowano maksymalny poziom ciśnienia akustycznego z filtrem A (LAFmax) oraz szczytowy poziom ciśnienia akustycznego z filtrem C (LCpeak). Równolegle rejestrowano zmienne operacyjne, a wszystkie pomiary poddano kalibracji i zdefiniowanym procedurom kontroli jakości przed analizą statystyczną. LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A; LAFmax, maksymalny poziom ciśnienia akustycznego z filtrem A; LCpeak, szczytowy poziom ciśnienia akustycznego z filtrem C. Prosimy o kliknięcie tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| Katedra | liczba łóżek | Średnia liczba zajętych łóżek | Powierzchnia pomieszczenia (m²) | Objętość pomieszczenia (m³) | Średnia liczba personelu | Średnia liczba odwiedzających | Średnia liczba alarmów | Średnia liczba rozmów |
| Ortopedia urazowa | 40 | 33.8 | 280.1 | 896.3 | 5.5 | 4.8 | 2.4 | 6.5 |
| Centrum Diagnostyki i Leczenia Piersi i Tarczycy | 60 | 50.6 | 416.8 | 1333.8 | 5.7 | 4.8 | 2.4 | 6.2 |
| Centrum Chorób Wrodzonych Serca | 35 | 29.4 | 243.6 | 779.5 | 5.7 | 5 | 2.4 | 6.5 |
| Chirurgia otolaryngologiczna i szyjna | 50 | 42.2 | 348.1 | 1113.9 | 5.8 | 5 | 2.5 | 6.7 |
| Oddział Diagnostyki i Leczenia Chorób Układu Pokarmowego I | 38 | 31.9 | 264.8 | 847.3 | 5.8 | 5 | 2.5 | 6.6 |
| II Oddział Diagnostyki i Leczenia Chorób Układu Pokarmowego | 52 | 43.9 | 361.7 | 1157.4 | 5.7 | 4.9 | 2.4 | 6.5 |
| Chirurgia wątroby, dróg żółciowych i trzustki | 48 | 40.5 | 334.1 | 1069.2 | 5.8 | 5.1 | 2.5 | 6.8 |
| Oddział Intensywnej Terapii II | 21 | 17.8 | 145.8 | 466.5 | 5.7 | 4.8 | 2.5 | 6.5 |
| Oddział Stawu Kolanowego II | 45 | 37.8 | 312.7 | 1000.6 | 5.8 | 4.9 | 2.5 | 6.5 |
| Oddział Choroby Parkinsona | 55 | 46.2 | 382.2 | 1223 | 5.9 | 5.1 | 2.5 | 6.8 |
| Katedra Medycyny Psychosomatycznej I | 68 | 57.1 | 472.6 | 1512.3 | 5.8 | 5 | 2.5 | 6.7 |
| Oddział Pediatryczny | 40 | 33.6 | 278.6 | 891.5 | 5.8 | 5.1 | 2.5 | 6.8 |
| Katedra Urologii | 62 | 52.1 | 431.2 | 1380 | 5.8 | 5 | 2.5 | 6.6 |
| Zakład Naczyń i Nerwów I | 72 | 60.6 | 500.3 | 1601 | 5.9 | 5.1 | 2.5 | 6.9 |
| Oddział Naczyniowy i Neurologiczny II | 60 | 50.4 | 417.5 | 1336 | 5.7 | 4.8 | 2.4 | 6.3 |
Tabela 1: Charakterystyka oddziałów szpitalnych objętych badaniem.Podstawowe cechy operacyjne 15 oddziałów stacjonarnych objętych badaniem hałasu środowiskowego. Zmienne na poziomie oddziału obejmują liczbę łóżek, średnią liczbę zajętych łóżek, powierzchnię pomieszczeń, objętość pomieszczeń oraz średnią liczbę personelu, odwiedzających, alarmów i rozmów zarejestrowanych podczas zestandaryzowanych okresów obserwacji w ciągu dnia. Zmienne operacyjne reprezentują wartości średnie uzyskane w predefinednych okresach obserwacji i zostały ocenione jako zmienne eksploracyjne w analizach regresji z efektami mieszanymi. Powierzchnia i objętość pomieszczeń są podane odpowiednio w metrach kwadratowych (m2) i metrach sześciennych (m3).
Opisowe charakterystyki akustyczne w strefach funkcjonalnych
Wśród wszystkich 1620 obserwacji najwyższe poziomy hałasu środowiskowego odnotowano na stanowiskach pielęgniarskich, a następnie w korytarzach, obszarach aktywności pacjentów i salach chorych (Tabela 2; Rysunki 2 i 3). Średnie wartości LAeq wynosiły 6,90 ± 2,34 dB(A) na stanowiskach pielęgniarskich, 63,97 ± 2,38 dB(A) w korytarzach, 61,38 ± 2,17 dB(A) w obszarach aktywności pacjentów oraz 60,38 ± 2,22 dB(A) na salach chorych. Podobne gradienty przestrzenne zaobserwowano dla LAFmax i LCpeak, co wskazuje na spójne różnice we wszystkich trzech metrykach akustycznych. Szczegółowe podsumowania opisowe wtórnych wyników akustycznych (LAFmax i LCpeak) przedstawiono w Tabeli uzupełniającej 1.
| Strefa funkcjonalna | Obserwacje (n) | LAeq, dB(A) Średnia ± SD | LAFmax, dB(A) średnia ± SD | LCpeak, dB(C) Średnia ± SD | Mediana LAeq (IQR), dB(A) |
| Oddział | 405 | 60.38 ± 2.22 | 67.95 ± 2.37 | 79.21 ± 3.31 | 60.4 (58.9–61.9) |
| Korytarz | 405 | 63.97 ± 2.38 | 71.46 ± 2.59 | 82.76 ± 3.56 | 64.0 (62.5–65.5) |
| Stanowisko pielęgniarskie | 405 | 66.90 ± 2.34 | 74.43 ± 2.56 | 85.72 ± 3.56 | 66.7 (65.4–68.3) |
| Strefa aktywności pacjenta | 405 | 61.38 ± 2.17 | 69.00 ± 2.42 | 80.20 ± 3.37 | 61.5 (60.0–62.7) |
Tabela 2: Opisowe charakterystyki akustyczne w zależności od funkcjonalnej strefy szpitala.Opisowe charakterystyki akustyczne dla czterech ustandaryzowanych funkcjonalnych stref szpitala. Wartości przedstawiono jako średnia ± odchylenie standardowe (SD), chyba że zaznaczono inaczej. Dla równoważnego ciągłego poziomu ciśnienia akustycznego z filtrem A (LAeq) podano medianę oraz rozstęp międzykwartylny (IQR). Każda strefa funkcjonalna dostarczyła 405 obserwacji do hierarchicznego zbioru danych (łącznie n = 1 620). LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A; LAFmax, maksymalny poziom ciśnienia akustycznego z filtrem A; LCpeak, szczytowy poziom ciśnienia akustycznego z filtrem C; IQR, rozstęp międzykwartylny.

Rycina 2: Mapa cieplna średnich równoważnych ciągłych poziomów ciśnienia akustycznego z korekcją A w różnych oddziałach i strefach funkcjonalnych szpitala. Mapa cieplna przedstawiająca średnie równoważne ciągłe poziomy ciśnienia akustycznego z filtrem A (LAeq) dla każdego oddziału szpitalnego i strefy funkcjonalnej. Wartości w komórkach reprezentują średnie LAeq w decybelach z filtrem A [dB(A)], a skala kolorów wskazuje narastający poziom dźwięku od wartości niższych do wyższych. Strefy funkcjonalne obejmują sale chorych, korytarze, stacje pielęgniarek oraz obszary aktywności pacjentów. CHD, wrodzone wady serca; ICU, oddział intensywnej terapii; LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Rozkład równoważnych ciągłych poziomów ciśnienia akustycznego ważonego krzywą A w ustandaryzowanych strefach funkcjonalnych szpitala. Wykresy pudełkowe przedstawiające rozkład równoważnych ciągłych poziomów ciśnienia akustycznego z filtrem A (LAeq) w czterech ustandaryzowanych strefach funkcjonalnych szpitala. Pudełka reprezentują rozstęp międzykwartylowy (percentyle 25.–75.), pozioma linia wewnątrz każdego pudełka oznacza medianę, a symbol × oznacza średnią. Wąsy rozciągają się do 1,5 krotności rozstępu międzykwartylowego (IQR), a obserwacje poza wąsami są przedstawione jako wartości odstające. Każda strefa funkcjonalna obejmuje 405 obserwacji (łącznie n = 1 620). Różnice między strefami funkcjonalnymi oceniano za pomocą podstawowego liniowego modelu efektów mieszanych z parami porównań skorygowanych testem Tukeya. LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A; IQR, rozstęp międzykwartylowy. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Mapa ciepła przedstawiona na Rysunku 2 wykazuje spójny wzorzec przestrzenny we wszystkich 15 uczestniczących oddziałach. Stanowiska pielęgniarskie konsekwentnie wykazywały najwyższe średnie wartości LAeq niezależnie od specjalizacji klinicznej, natomiast obszary sal chorych wykazywały najniższe równoważne ciągłe poziomy ciśnienia akustycznego ważone krzywą A. Oddziały o większej liczbie łóżek, w tym I Oddział Naczyniowy i Neurologiczny, I Oddział Medycyny Psychosomatycznej oraz Oddział Urologiczny, wykazywały ogólnie wyższe średnie wartości LAeq w poszczególnych strefach funkcjonalnych niż oddziały o mniejszej liczbie łóżek, takie jak II Oddział Intensywnej Terapii oraz Centrum Wad Wrodzonych Serca. Mimo zaobserwowania niewielkiej zmienności międzyoddziałowej, względny ranking stref funkcjonalnych pozostał spójny. Średnie wartości LAeq na stanowiskach pielęgniarskich wahały się od 64,1 dB(A) na II Oddziale Intensywnej Terapii do 68,9 dB(A) na I Oddziale Naczyniowym i Neurologicznym, podczas gdy pomiary w salach chorych mieściły się w zakresie od 57,9 do 61,6 dB(A). Obszary aktywności pacjentów i korytarze wykazywały wartości pośrednie w prawie wszystkich oddziałach (Rysunek 2). Rysunek 3 ilustruje rozkład pomiarów LAeq w obrębie każdej strefy funkcjonalnej. Mediany wartości LAeq były najwyższe na stanowiskach pielęgniarskich (6,7 dB[A]), następnie na korytarzach (64,0 dB[A]), w obszarach aktywności pacjentów (61,5 dB[A]) oraz w salach chorych (60,4 dB[A]). Chociaż we wszystkich strefach funkcjonalnych zidentyfikowano sporadyczne obserwacje wysokich wartości, rozstępy międzykwartylowe były stosunkowo wąskie, co wskazuje na ograniczoną zmienność w każdym środowisku funkcjonalnym. Wykresy pudełkowe wykazują wyraźną separację rozkładu wartości LAeq pomiędzy czterema strefami funkcjonalnymi (Rysunek 3).
Analiza porównawcza stref funkcjonalnych z wykorzystaniem modeli mieszanych
Aby uwzględnić hierarchiczną strukturę powtarzalnych pomiarów akustycznych w obrębie oddziałów szpitalnych, w analizie głównej zastosowano liniowy model mieszany, w którym oddział szpitalny określono jako losową część stałej (random intercept), a strefę funkcjonalną jako efekt stały (Tabela 3). Takie podejście do modelowania pozwoliło uwzględnić grupowanie wewnątrz oddziałów przy jednoczesnym oszacowaniu skorygowanych różnic w LAeq pomiędzy czterema zdefiniowanymi strefami funkcjonalnymi.
| Efekty stałe (Wynik: LAeq) |
| Efekt stały | β | SE | 95% CI | Wartość P |
| Punkt przecięcia (Korytarz) | 63.972 | 0.26 | 63.451 do 64.493 | <0.01 |
| Stanowisko pielęgniarskie | 2.931 | 0.146 | 2.64 do 3.218 | <0.01 |
| Strefa aktywności pacjentów | −2.594 | 0.146 | −2.81 do −2.307 | <0.01 |
| Sala chorych | −3.589 | 0.146 | −3.876 do −3.302 | <0.01 |
| Szacowane średnie marginalne |
| Strefa funkcjonalna | Szacowana średnia LAeq, dB(A) |
| Stanowisko pielęgniarskie | 6.9 |
| Korytarz | 63.97 |
| Strefa aktywności pacjentów | 61.38 |
| Sala chorych | 60.38 |
| Porównania parowe z korektą Tukeya |
| Porównanie | Różnica średnich, dB(A) | Skorygowana wartość P |
| Stanowisko pielęgniarskie vs Korytarz | 2.93 | <0.01 |
| Stanowisko pielęgniarskie vs Strefa aktywności pacjentów | 5.52 | <0.001 |
| Stanowisko pielęgniarskie vs Sala chorych | 6.52 | <0.01 |
| Korytarz vs Strefa aktywności pacjentów | 2.59 | <0.001 |
| Korytarz vs Sala chorych | 3.59 | <0.01 |
| Strefa aktywności pacjentów vs Sala chorych | 0.9 | <0.001 |
| Wydajność modelu |
| Statystyka | Wartość |
| Współczynnik korelacji wewnątrzklasowej (ICC) | 0.172 |
| Marginalny R² | 0.53 |
| Warunkowy R² | 0.61 |
Tabela 3: Główny liniowy model mieszany efektów porównujący hałas środowiskowy w funkcjonalnych strefach szpitala.Wyniki głównego liniowego modelu mieszanego efektów oceniającego różnice w równoważnych ciągłych poziomach ciśnienia akustycznego z filtrem A (LAeq) pomiędzy czterema ustandaryzowanymi funkcjonalnymi strefami szpitala. Oddział szpitalny określono jako losowy punkt przecięcia (random intercept), aby uwzględnić grupowanie powtórzonych obserwacji, natomiast strefę funkcjonalną włączono jako efekt stały, przyjmując korytarz jako kategorię referencyjną. Dla efektów stałych przedstawiono współczynniki regresji (β), błędy standardowe (SE), 95-procentowe przedziały ufności (CI) Walda oraz dwustronne wartości P. Szacowane średnie marginalne porównano za pomocą parowych porównań z korektą Tukeya. Wydajność modelu podsumowano za pomocą współczynnika korelacji wewnątrzklasowej (ICC) oraz marginalnego i warunkowego współczynnika determinacji (R2), obliczonych zgodnie z metodą Nakagawy i Schielzetha. LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A; β, współczynnik regresji; SE, błąd standardowy; CI, przedział ufności; ICC, współczynnik korelacji wewnątrzklasowej.
Analiza efektów mieszanych wykazała wysoce istotny ogólny wpływ strefy funkcjonalnej na poziom hałasu środowiskowego (ogólny test ilorazu wiarygodności, P < 0,01). Przy przyjęciu korytarza jako kategorii referencyjnej, na stanowiskach pielęgniarskich odnotowano istotnie wyższe wartości LAeq (β = 2,931 dB[A], SE = 0,146, 95% CI = 2,64–3,218, P < 0,01), natomiast obszary aktywności pacjentów (β = −2,594 dB[A], SE = 0,146, 95% CI = −2,81 do −2,307, P < 0,01) oraz sale chorych (β = −3,589 dB[A], SE = 0,146, 95% CI = −3,876 do −3,302, P < 0,01) wykazały istotnie niższe wartości LAeq (Tabela 3). Oszacowane średnie marginalne wynikające z modelu efektów mieszanych wykazały wyraźny hierarchiczny wzorzec hałasu środowiskowego w poszczególnych strefach funkcjonalnych. Skorygowana średnia LAeq była najwyższa na stanowiskach pielęgniarskich (6,90 dB[A]), następnie na korytarzach (63,97 dB[A]), obszarach aktywności pacjentów (61,38 dB[A]) oraz w salach chorych (60,38 dB[A]). Wszystkie porównania parowe pozostały istotne statystycznie po korekcie Tukeya (wszystkie skorygowane P < 0,01), co wskazuje na statystycznie istotne różnice w LAeq pomiędzy wszystkimi czterema strefami funkcjonalnymi. Różnice średnich wynosiły od 0,9 dB(A) pomiędzy obszarami aktywności pacjentów a salami chorych do 6,52 dB(A) pomiędzy stanowiskami pielęgniarskimi a salami chorymi. Rysunek 4 podsumowuje szacunki skorygowane zgodnie z modelem. Oszacowane średnie marginalne oraz ich 95% przedziały ufności wykazały minimalne nakładanie się pomiędzy strefami funkcjonalnymi, a grupowanie Tukeya przypisało każdą strefę do odrębnej grupy istotności (A–D). Stanowiska pielęgniarskie przypisano do grupy A z najwyższym skorygowanym LAeq, następnie korytarze (grupa B), obszary aktywności pacjentów (grupa C) i sale chorych (grupa D). Analiza składników wariancji wykazała, że około 17,2% całkowitej zmienności LAeq wynikało z różnic pomiędzy oddziałami szpitalnymi (współczynnik korelacji wewnątrzklasowej = 0,172), co uzasadniało wprowadzenie losowego punktu przecięcia na poziomie oddziału. Efekty stałe wyjaśniły 53% zmienności hałasu środowiskowego (marginalny R2 = 0,53), podczas gdy połączony model efektów stałych i losowych wyjaśnił 61% całkowitej zmienności (warunkowy R2 = 0,61), co wskazuje na dobrą zdolność wyjaśniającą modelu hierarchicznego (Tabela 3).

Rysunek 4: Szacowane średnie marginalne równoważnych ciągłych poziomów ciśnienia akustycznego z korekcją A w różnych strefach funkcjonalnych szpitala. Szacowane średnie marginalne (EMM) równoważnych ciągłych poziomów ciśnienia akustycznego z poprawką A (LAeq) w czterech strefach funkcjonalnych szpitala, wyznaczone na podstawie głównego liniowego modelu efektów mieszanych. Punkty reprezentują szacowane przez model średnie marginalne, a poziome słupki błędów oznaczają 95-procentowe przedziały ufności (CI). Oddział szpitalny został uwzględniony jako efekt losowy, a strefa funkcjonalna jako efekt stały. Porównania parowe przeprowadzono przy użyciu szacowanych średnich marginalnych z korektą Tukeya; strefy funkcjonalne oznaczone różnymi literami reprezentują różnice istotne statystycznie (P < 0,05). Każda strefa funkcjonalna dostarczyła 405 obserwacji do analizy. LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A; EMM, szacowana średnia marginalna; CI, przedział ufności. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Oszczędna regresja z mieszanymi efektami dla charakterystyk operacyjnych powiązanych z LAeq
Aby zbadać charakterystyki operacyjne związane z hałasem środowiskowym przy jednoczesnym zminimalizowaniu współliniowości, zastosowano uprzednio określony oszczędny model mieszany, w którym strefę funkcjonalną, liczbę odwiedzających oraz liczbę łóżek na oddziale uwzględniono jako efekty stałe, a oddział szpitalny pozostawiono jako losowy punkt przecięcia (Tabela 4). Wybór tych zmiennych oparto na eksploracyjnej macierzy korelacji (Rycina uzupełniająca 1), która wykazała silne korelacje pomiędzy kilkoma zmiennymi strukturalnymi, w szczególności powierzchnią i objętością pomieszczenia (ρ Spearmana = 0,94). W związku z tym wysoko skorelowane zmienne strukturalne wykluczono z modelu oszczędnego, aby poprawić jego interpretowalność i stabilność. Po skorygowaniu o operacyjne współzmienne strefa funkcjonalna pozostała najsilniejszym niezależnym czynnikiem związanym z hałasem środowiskowym. W porównaniu z korytarzami, stanowiska pielęgniarek wiązały się ze skorygowanym wzrostem o 2,93 dB(A) (SE = 0,146, 95% CI = 2,647–3,20, P < 0,01), natomiast obszary aktywności pacjentów (β = −2,602 dB[A], SE = 0,146, 95% CI = −2,8 do −2,315, P < 0,01) i sale chorych (β = −3,606 dB[A], SE = 0,146, 95% CI = −3,892 do −3,319, P < 0,01) pozostawały znacznie cichsze niż korytarze.
| Predyktor | Współczynnik β | SE | 95% CI | Wartość p |
| Przesiekanie (Korytarz) | 63.97 | 0.27 | 63,45 do 64,49 | <0.001 |
| Stanowisko pielęgniarskie (w porównaniu z korytarzem) | 2.933 | 0.146 | 2,647 do 3,20 | <0.001 |
| Strefa aktywności pacjenta (w porównaniu z korytarzem) | −2.602 | 0.146 | −2,88 do −2,315 | <0.001 |
| Sala chorych (w porównaniu do korytarza) | −3.606 | 0.146 | −3,892 do −3,319 | <0.001 |
| Liczba odwiedzających (na każdego dodatkowego odwiedzającego) | 0.06 | 0.026 | 0,08 do 0,12 | 0.024 |
| Wydajność łóżek (na każde dodatkowe łóżko) | 0.064 | 0.005 | 0,05 do 0,073 | <0.001 |
| Model: LAeq ~ Strefa Funkcjonalna + Liczba Odwiedzających + Liczba Łóżek + (1 | Oddział) |
Tabela 4: Oszczędna regresja liniowa z efektami mieszanymi dla charakterystyk operacyjnych związanych z hałasem środowiskowym.Wyniki uprzednio zdefiniowanego oszczędnego modelu regresji liniowej z efektami mieszanymi, oceniającego związki między równoważnym ciągłym poziomem ciśnienia akustycznego z poprawką A (LAeq) a wybranymi charakterystykami operacyjnymi. Model obejmował strefę funkcjonalną, liczbę odwiedzających oraz liczbę łóżek oddziałowych jako efekty stałe, a oddział szpitalny jako losowy punkt przecięcia w celu uwzględnienia grupowania powtarzających się obserwacji (model: LAeq ~ Functional Zone + Visitor Count + Bed Capacity + (1 | Department)). Korytarz określono jako kategorię referencyjną dla strefy funkcjonalnej, a liczbę odwiedzających oraz liczbę łóżek modelowano jako zmienne ciągłe. Przedstawiono współczynniki regresji (β), błędy standardowe (SE), 95% przedziały ufności (CI) Walda oraz dwustronne wartości P. Współczynniki regresji dla liczby odwiedzających i liczby łóżek reprezentują szacowaną zmianę LAeq związaną z każdym dodatkowym odwiedzającym oraz każdym dodatkowym łóżkiem szpitalnym, odpowiednio. Punkt przecięcia reprezentuje szacowaną średnią wartość LAeq dla kategorii referencyjnej (korytarz), gdy predyktory ciągłe są równe zero. Tabela ta przedstawia uprzednio zdefiniowany oszczędny model opracowany w celu zminimalizowania nadmiernego dopasowania modelu; w pełni skorygowany model wrażliwości obejmujący wszystkie zarejestrowane charakterystyki operacyjne przedstawiono w tabeli uzupełniającej S2. LAeq, równoważny ciągły poziom ciśnienia akustycznego z poprawką A; β, współczynnik regresji; SE, błąd standardowy; CI, przedział ufności.
Wśród ciągłych predyktorów operacyjnych liczba odwiedzających wykazała niezależny pozytywny związek z LAeq. Każdy dodatkowy odwiedzający obecny podczas pomiaru wiązał się z szacowanym wzrostem o 0,060 dB(A) (SE = 0,026, 95% CI = 0,08–0,12, P = 0,024). Podobnie liczba łóżek na oddziale pozostała niezależnie związana z hałasem środowiskowym, przy czym każde dodatkowe łóżko szpitalne wiązało się z szacowanym wzrostem LAeq o 0,064 dB(A) (SE = 0,05, 95% CI = 0,05–0,073, P < 0,001).
Rysunek 5 przedstawia dopasowane predykcje wygenerowane z oszczędnego modelu mieszanego. Przewidywane LAeq rosło w przybliżeniu liniowo w obserwowanym zakresie liczby łóżek na oddziale, a większa liczba odwiedzających konsekwentnie przesuwała przewidywane poziomy dźwięku w górę w całym zakresie liczby łóżek. Równoległe dopasowane linie regresji wskazują na addytywne efekty liczby odwiedzających oraz liczby łóżek w dopasowanym modelu po skorygowaniu o strefę funkcjonalną i grupowanie na poziomie oddziału. Predykcje oparte na modelu wykazały, że dla oddziału z około 20 łóżkami przewidywane LAeq wynosiło od 61,6 dB(A) przy niskiej liczbie odwiedzających (trzech odwiedzających) do 62,4 dB(A) przy wyższej liczbie odwiedzających (siedmiu odwiedzających). W górnej granicy obserwowanego zakresu liczby łóżek (72 łóżka) odpowiadające im przewidywane wartości LAeq wzrosły odpowiednio do 64,5 dB(A) i 65,3 dB(A), co obrazuje połączony wkład liczby łóżek na oddziale oraz liczby odwiedzających w dopasowanym modelu.

Rysunek 5: Prognozy oparte na modelu w zakresie równoważnych ciągłych poziomów ciśnienia akustycznego z korekcją A w zależności od liczby łóżek na oddziale. Przewidywane równoważne ciągłe poziomy ciśnienia akustycznego z korekcją A (LAeq) w obserwowanym zakresie liczby łóżek oddziałowych, wyznaczone na podstawie oszczędnego liniowego modelu regresji z efektami mieszanymi. Kolorowe linie przedstawiają wartości LAeq przewidywane przez model dla niskiej (3 odwiedzających), średniej (5 odwiedzających) i wysokiej (7 odwiedzających) liczby odwiedzających, przy zachowaniu stałych wartości pozostałych parametrów modelu. Zacieniowane obszary wskazują 95-procentowe przedziały ufności dla wartości przewidywanych, a punkty reprezentują pomiary obserwacyjne. Przewidywania przedstawiono wyłącznie w obserwowanym zakresie liczby łóżek oddziałowych uwzględnionych w badaniu. LAeq, równoważny ciągły poziom ciśnienia akustycznego z korekcją A; CI, przedział ufności. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.
Analizy wrażliwości z wykorzystaniem w pełni skorygowanego modelu efektów mieszanych (Supplementary Table 2), który dodatkowo uwzględniał liczbę zajętych łóżek, powierzchnię pokoju, objętość pokoju, liczbę personelu, zdarzenia alarmowe, zdarzenia konwersacyjne oraz okres pomiaru, dostarczyły wysoce spójnych szacunków dla strefy funkcjonalnej i liczby odwiedzających. Stanowiska pielęgniarskie pozostały najgłośniejszą strefą funkcjonalną, a liczba odwiedzających nadal była niezależnie skorelowana z wyższym LAeq (β = 0.060, P = 0.023), podczas gdy liczba zajętych łóżek, powierzchnia pokoju, objętość pokoju, liczba personelu, zdarzenia alarmowe, zdarzenia konwersacyjne oraz okres pomiaru nie wykazały niezależnego związku z hałasem środowiskowym po wzajemnej korekcie. Wyniki te potwierdzają rzetelność modelu oszczędnego i wskazują, że strefa funkcjonalna oraz liczba odwiedzających były głównymi zmiennymi powiązanymi z hałasem środowiskowym w niniejszym badaniu.
Porównanie opisowe z referencyjnymi kryteriami akustycznymi
Aby osadzić zaobserwowane poziomy hałasu środowiskowego w odpowiednim kontekście, zmierzone wartości LAeq porównano opisowo z trzema powszechnie cytowanymi referencyjnymi kryteriami akustycznymi: dzienną wartością wytyczną WHO dla środowisk opieki nad pacjentami w szpitalach (35 dB[A]), wartością referencyjną EN ISO 16032 dla ogólnych środowisk wewnętrznych (40 dB[A]) oraz chińską standardową wartością referencyjną jakości środowiska GB3096-202 Klasa 1 (45 dB[A]). Ponieważ wartości referencyjne te różnią się zakresem, przeznaczeniem, metodologią pomiaru oraz okresem uśredniania, porównania przeprowadzono wyłącznie w celach opisowych i nie stanowią one formalnych ocen zgodności z przepisami (Rysunek 6). We wszystkich 1620 obserwacjach akustycznych poziomy hałasu środowiskowego przekroczyły każdą z wybranych wartości referencyjnych. Ogółem 1615 obserwacji (9,7%) przekroczyło dzienną wartość wytyczną WHO wynoszącą 35 dB(A), 1586 obserwacji (97,9%) przekroczyło wartość referencyjną EN ISO 16032 wynoszącą 40 dB(A), a 1513 obserwacji (93,2%) przekroczyło chińską wartość referencyjną GB3096-202 Klasa 1 wynoszącą 45 dB(A).

Rycina 6: Procent obserwacji akustycznych przekraczających wybrane międzynarodowe i krajowe referencyjne kryteria akustyczne w poszczególnych strefach funkcjonalnych szpitala. Procent obserwacji równoważnego ciągłego poziomu ciśnienia akustycznego z filtrem A (LAeq) przekraczających wybrane międzynarodowe i krajowe referencyjne kryteria akustyczne w czterech zestandaryzowanych strefach funkcjonalnych szpitala. Słupki reprezentują proporcję obserwacji przekraczających wartości referencyjne wytycznych Światowej Organizacji Zdrowia (WHO) (>35 dB[A]), EN ISO 16032 (>40 dB[A]), oraz chińskiej normy krajowej GB3096-202 (>45 dB[A]). Procenty obliczono jako liczbę obserwacji przekraczających dany próg referencyjny podzieloną przez całkowitą liczbę obserwacji w każdej strefie funkcjonalnej (n = 405). Ponieważ dokumenty referencyjne te różnią się przeznaczeniem, metodologią pomiaru, warunkami środowiskowymi oraz okresem uśredniania, porównania przedstawiono wyłącznie w celach opisowych i nie stanowią one formalnych ocen zgodności. LAeq, równoważny ciągły poziom ciśnienia akustycznego ważony krzywą A; WHO, Światowa Organizacja Zdrowia. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Zaobserwowano wyraźne różnice pomiędzy strefami funkcjonalnymi. W punktach pielęgniarskich odnotowano najwyższy odsetek pomiarów przekraczających każdą z wartości referencyjnych, przy czym 10,0% przekroczyło wartości referencyjne WHO oraz EN ISO, a 9,0% przekroczyło chińską wartość referencyjną. Podobnie wysokie odsetki przekroczeń zaobserwowano w korytarzach (odpowiednio 10,0%, 9,3% i 96,0%). W porównaniu z nimi, w salach chorych odnotowano najniższe — choć wciąż znaczące — odsetki pomiarów przekraczających trzy wartości referencyjne (odpowiednio 9,3%, 95,3% i 87,9%), natomiast w strefach aktywności pacjentów odsetki te były średnie (odpowiednio 9,5%, 97,0% i 90,6%). Te wyniki opisowe były zgodne z analizami efektów mieszanych, w których punkty pielęgniarskie wykazały najwyższe poziomy hałasu środowiskowego, a sale chorych najniższe spośród czterech stref funkcjonalnych. Wysokie odsetki pomiarów przekraczających wszystkie trzy wartości referencyjne wskazują, że podwyższony dzienny hałas środowiskowy występował w całym szpitalu, a nie tylko w poszczególnych oddziałach czy strefach funkcjonalnych.
Analizy wrażliwości i diagnostyka modelu
Krzewność głównych wyników oceniono przy użyciu w pełni skorygowanego modelu mieszanego wrażliwości, który dodatkowo uwzględniał liczbę zajętych łóżek, powierzchnię pomieszczeń, objętość pomieszczeń, liczbę personelu, zdarzenia alarmowe, zdarzenia konwersacyjne oraz okres pomiarowy (Supplementary Table 2). Szacunki dla stref funkcjonalnych pozostały wysoce zgodne z tymi zaobserwowanymi w uprzednio określonym modelu oszczędnym. Punkty pielęgniarskie nadal wykazywały istotnie wyższe LAeq niż korytarze, podczas gdy strefy aktywności pacjentów i sale chorych pozostały istotnie cichsze (wszystkie P < 0,01). Liczba odwiedzających również pozostała niezależnie powiązana z wyższym LAeq (β = 0,060, P = 0,023), natomiast liczba zajętych łóżek, powierzchnia pomieszczeń, objętość pomieszczeń, liczba personelu, zdarzenia alarmowe, zdarzenia konwersacyjne oraz okres pomiarowy nie były niezależnie powiązane z hałasem środowiskowym po jednoczesnej korekcie (Supplementary Table 2). Wyniki te wskazują, że główne wnioski były odporne na dodatkową korektę o skorelowane cechy operacyjne. Założenia modelu oceniono za pomocą procedur diagnostyki reszt (Supplementary Figure 1). Wykres kwantyli normalnych (Q–Q) wykazał bliską zgodność między zaobserwowanymi a teoretycznymi kwantylami reszt, co wskazuje na brak istotnego odejścia od rozkładu normalnego. Wykres reszt względem wartości dopasowanych wykazał losowy rozrzut wokół zera bez dostrzegalnych trendów ani dowodów na heteroscedastyczność, co potwierdza założenia liniowości i stałej wariancji reszt. Zgodnie z tymi obserwacjami wizualnymi, test Shapiro-Wilka nie wykazał statystycznie istotnego odchylenia od normalności (W = 0,992, P = 0,084). Dodatkowe statystyki wydajności modelu podsumowano w Supplementary Table 3. Standardowe odchylenie reszt wynosiło 1,92 dB(A), a standardowe odchylenie losowego punktu przecięcia na poziomie oddziału wynosiło 0,90 dB(A). Współczynnik korelacji wewnątrzklasowej wyniósł 0,172, co wskazuje na umiarkowaną grupowanie obserwacji w obrębie oddziałów szpitalnych. Model mieszany wyjaśnił 53% zmienności za pomocą samych efektów stałych (marginalny R2 = 0,53) oraz 61%, gdy włączono losowe efekty na poziomie oddziału (warunkowy R2 = 0,61). Czynniki inflacji wariancji były jednolicie niskie (maksymalne VIF = 2,14), co wskazuje na brak problematycznej współliniowości wśród zachowanych predyktorów. Co więcej, wszystkie wartości odległości Cooka były <1,0 (maksymalnie = 0,18), co sugeruje, że żadna pojedyncza obserwacja nie wywierała nadmiernego wpływu na dopasowany model. Eksploracyjna macierz korelacji (Supplementary Figure 2) wykazała oczekiwane pozytywne powiązania między kilkoma cechami operacyjnymi. Silne korelacje zaobserwowano między powierzchnią a objętością pomieszczeń (ρ Spearmana = 0,94), liczbą łóżek a liczbą zajętych łóżek (ρ = 0,86) oraz liczbą łóżek a objętością pomieszczeń (ρ = 0,81). Umiarkowane korelacje zaobserwowano między liczbą odwiedzających a liczbą zajętych łóżek (ρ = 0,71) oraz między liczbą personelu a zdarzeniami konwersacyjnymi (ρ = 0,62). Wyniki te potwierdziły przyjętą strategię selekcji zmiennych dla oszczędnego modelu mieszanego, w której wysoce skorelowane zmienne strukturalne zostały wykluczone w celu zminimalizowania współliniowości przy zachowaniu interpretowalności modelu.
Ogólnie rzecz biorąc, główne analizy z efektami mieszanymi, analizy wrażliwości, diagnostyka reszt oraz eksploracyjna ocena korelacji wykazały wewnętrznie spójne wyniki. we wszystkich analizach strefa funkcjonalna wykazała najsilniejszy niezależny związek z hałasem środowiskowym, podczas gdy liczba odwiedzających oraz liczba łóżek na oddziale również pozostały niezależnie powiązane z LAeq. Większość pozostałych zarejestrowanych parametrów operacyjnych nie była niezależnie powiązana z hałasem środowiskowym po wzajemnej korekcie.
Dostępność danych:
Dane wspierające wyniki niniejszego badania są dostępne w treści artykułu oraz w jego materiałach uzupełniających. Kompletny zestaw surowych pomiarów akustycznych, zapisy kalibracji, zestaw danych zmiennych operacyjnych, dane źródłowe do rycin oraz kod statystycznej analizy w programie R zostały zdeponowane w repozytorium Zenodo i są publicznie dostępne pod adresem https://doi.org/10.5281/zenodo.21404807. Tabele uzupełniające 1–3 oraz Ryciny uzupełniające 1 i 2 zawierają dodatkowe przetworzone zestawy danych, analizy wrażliwości, diagnostykę modeli oraz analizy pomocnicze.
Rycina uzupełniająca 1: Wykresy diagnostyczne reszt dla oszczędnego liniowego modelu efektów mieszanych. Wykresy diagnostyczne oceniające założenia oszczędnego liniowego modelu regresji z efektami mieszanymi. (A) Wykres kwantylowy (Q–Q) reszt warunkowych wykazujący przybliżoną normalność rozkładu reszt. (B) Wykres reszt warunkowych w funkcji wartości dopasowanych, nie wykazujący istotnych odchyleń od homoscedastyczności ani systematycznych wzorców w całym zakresie dopasowania. Oddział szpitalny został uwzględniony jako efekt losowy (random intercept), natomiast strefa funkcjonalna, liczba odwiedzających oraz liczba łóżek zostały uwzględnione jako efekty stałe. Diagnostyka reszt potwierdza zasadność założeń modelu dla głównej analizy regresji. Wykresy diagnostyczne zostały wygenerowane przy użyciu pakietu DHARMa w programie R. LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A.Prosimy o kliknięcie tutaj, aby pobrać ten plik.
Rycina uzupełniająca 2: Macierz korelacji Spearmana dla zmiennych operacyjnych uwzględnionych w eksploracyjnych analizach z efektami mieszanymi. Macierz korelacji Spearmana przedstawiająca powiązania parzyste między operacyjnymi zmiennymi oddziałowymi ocenianymi podczas eksploracyjnych analiz regresji z efektami mieszanymi. Wartości w komórkach reprezentują współczynniki korelacji rang Spearmana (ρ), a skala kolorystyczna wskazuje siłę i kierunek korelacji, gdzie kolor niebieski oznacza korelacje dodatnie, a czerwony korelacje ujemne. Analizę korelacji przeprowadzono w celu oceny potencjalnej współliniowości między kandydującymi zmiennymi objaśniającymi przed opracowaniem modelu. Zmienne wykazujące silną współliniowość zostały uwzględnione podczas wyboru modelu, aby zminimalizować współliniowość w oszczędnym modelu regresji. Zmienne operacyjne obejmowały liczbę łóżek na oddziale, liczbę zajętych łóżek, powierzchnię pokoju, objętość pokoju, liczbę personelu, liczbę odwiedzających, liczbę alarmów oraz liczbę rozmów. Spearman’s ρ, współczynnik korelacji rang Spearmana.Proszę kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 1: Wtórne wyniki akustyczne w zależności od funkcjonalnej strefy szpitala.Podsumowanie opisowe wtórnych wyników akustycznych zmierzonych w czterech zestandaryzowanych funkcjonalnych strefach szpitala. Wartości przedstawiono jako średnia ± odchylenie standardowe (SD), reprezentując zbiorcze obserwacje ze wszystkich 15 oddziałów stacjonarnych. Równoważny ciągły poziom ciśnienia akustycznego filtrowany krzywą A (LAeq) był głównym wynikiem akustycznym, natomiast maksymalny poziom ciśnienia akustycznego filtrowany krzywą A (LAFmax) oraz szczytowy poziom ciśnienia akustycznego filtrowany krzywą C (LCpeak) zostały zapisane jako wtórne parametry akustyczne w celu charakterystyki maksymalnych i szczytowych poziomów dźwięku podczas rutynowego funkcjonowania szpitala w ciągu dnia. LAeq, równoważny ciągły poziom ciśnienia akustycznego filtrowany krzywą A; LAFmax, maksymalny poziom ciśnienia akustycznego filtrowany krzywą A; LCpeak, szczytowy poziom ciśnienia akustycznego filtrowany krzywą C; SD, odchylenie standardowe.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 2: W pełni skorygowany liniowy model regresji z efektami mieszanymi dla analizy wrażliwości. Wyniki w pełni skorygowanego modelu regresji liniowej z efektami mieszanymi w analizie wrażliwości, oceniającego powiązania między równoważnym ciągłym poziomem ciśnienia akustycznego z filtrem A (LAeq) a wszystkimi odnotowanymi charakterystykami operacyjnymi. Oddział szpitalny został uwzględniony jako losowy punkt przecięcia w celu uwzględnienia grupowania powtarzanych obserwacji. Jako efekty stałe uwzględniono strefę funkcjonalną, przedział czasowy, liczbę łóżek na oddziale, liczbę zajętych łóżek, powierzchnię pokoju, objętość pokoju, liczbę personelu, liczbę odwiedzających, liczbę alarmów oraz liczbę rozmów. Korytarz posłużył jako kategoria referencyjna dla strefy funkcjonalnej, a popołudnie jako kategoria referencyjna dla przedziału czasowego. Współczynniki regresji (β), błędy standardowe (SE), 95% przedziały ufności (CI) Walda oraz dwustronne P podano wartości. Model ten został opracowany w ramach analizy wrażliwości w celu oceny odporności modelu oszczędnego przedstawionego w tabeli 4 po skorygowaniu o wszystkie zarejestrowane charakterystyki operacyjne. LAeq, równoważny ciągły poziom ciśnienia akustycznego z ważeniem A; β, współczynnik regresji; SE, błąd standardowy; CI, przedział ufności.Prosimy kliknąć tutaj, aby pobrać ten plik.
Tabela uzupełniająca 3: Wydajność modelu i statystyki diagnostyczne dla głównego liniowego modelu efektów mieszanych. Podsumowanie wydajności modelu i statystyk diagnostycznych dla głównego liniowego modelu mieszanego, oceniającego równoważne ciągłe poziomy ciśnienia akustycznego ważone krzywą A (LAeq) w różnych strefach funkcjonalnych szpitala. Odchylenie standardowe reszt (σ) reprezentuje niewyjaśnioną zmienność wewnątrzmodelową, natomiast odchylenie standardowe losowego punktu przecięcia określa zmienność międzyoddziałową. Współczynnik korelacji wewnątrzklasowej (ICC) szacuje proporcję całkowitej wariancji przypisaną do grupowania w obrębie oddziałów szpitalnych. Marginalny R2 reprezentuje proporcję wariancji wyjaśnionej wyłącznie przez efekty stałe, podczas gdy warunkowe R2 reprezentuje proporcję wariancji wyjaśnionej zarówno przez efekty stałe, jak i losowe. Kryterium informacyjne Akaike (AIC), kryterium informacyjne bayesowskie (BIC) oraz logarytm wiarygodności zostały podane jako miary dopasowania modelu. Czynniki inflacji wariancji (VIF) były wszystkie <5, co wskazuje na brak dowodów na problematyczną współliniowość, a wartości odległości Cooka wynosiły <1,0 dla wszystkich obserwacji, co wskazuje, że żadna pojedyncza obserwacja nie wywarła nadmiernego wpływu na szacunki modelu. LAeq, równoważny ciągły poziom ciśnienia akustycznego z filtrem A; ICC, współczynnik korelacji wewnątrzklasowej; VIF, czynnik amplifikacji wariancji.Kliknij tutaj, aby pobrać ten plik.