Studiekenmerken en meetkader
Metingen van omgevingsgeluid werden uitgevoerd in 15 klinische afdelingen van een groot tertiair opleidingsziekenhuis tussen maart en augustus 2025 met behulp van een gestandaardiseerd hiërarchisch bemonsteringsprotocol (Figuur 1). Binnen elke afdeling werden vier vooraf gedefinieerde functionele zones (zaal, patiëntenactiviteitsruimte, gang en verpleegpost) geëvalueerd. Binnen elke functionele zone werden drie vaste bemonsteringslocaties vastgesteld, en metingen werden verkregen gedurende drie gestandaardiseerde observatieperioden overdag (09:00–19:00). Op elke locatie en tijdens elke observatieperiode werden drie opeenvolgende technische replica-metingen geregistreerd met een gekalibreerde Klasse 1 geluidsdrukmeter, waarbij akoestische gegevens voor LAeq, LAFmax en LCpeak werden gegenereerd. Operationele variabelen, waaronder bedcapaciteit, bezette bedden, personeelsaantal, bezoekersaantal, alarmgebeurtenissen, gespreksgebeurtenissen, kameroppervlakte en kamervolume, werden gelijktijdig met elke akoestische observatie geregistreerd. Het geplande bemonsteringskader bestond uit 1.620 akoestische observaties (15 afdelingen × 4 functionele zones × 3 bemonsteringslocaties × 3 observatieperioden × 3 technische replica's). Eén geplande observatie in de patiëntenactiviteitsruimte van de afdeling Intensive Care II kon niet worden voltooid vanwege tijdelijke beperkingen in de toegang voor spoedeisende zorg en werd behandeld met behulp van complete-case-analyse zonder data-imputatie, zoals vooraf gespecificeerd in het statistische analyseplan. Aanvullende metingen in binnen- en buitenruimtes voor het publiek werden verkregen voor beschrijvende contextuele vergelijking, maar werden uitgesloten van de mixed-effects-analyses op afdelingsniveau (Figuur 1). De kenmerken van de deelnemende afdelingen zijn samengevat in Tabel 1. De bedcapaciteit varieerde van 21 tot 72 bedden, met een overeenkomstig gemiddeld aantal bezette bedden variërend van 17,8 tot 60,6. De gemiddelde personeelsaantallen waren relatief consistent tussen de afdelingen (5,5–5,9 personen per observatie), terwijl het aantal bezoekers gemiddeld ongeveer vijf personen per observatie bedroeg. De oppervlakten van de afdelingskamers varieerden van 145,8 tot 500,3 m2, en de kamervolumes varieerden van 466,5 tot 1.601,0 m3, wat wijst op aanzienlijke variatie in de fysieke kenmerken van de deelnemende afdelingen.

Figuur 1: Studieontwerp en hiërarchisch protocol voor het meten van omgevingsgeluid. Overzicht van het studieontwerp en het hiërarchische protocol voor het meten van omgevingsgeluid. Het omgevingsgeluid werd beoordeeld in 15 klinische afdelingen van een tertiair opleidingsziekenhuis met behulp van een meertraps hiërarchisch bemonsteringskader. Metingen werden uitgevoerd in vier gestandaardiseerde functionele zones (ziekenzaal, gang, verpleegpost en patiëntenactiviteitsruimte), op drie vooraf gedefinieerde bemonsteringslocaties (ingang, midden en einde), gedurende drie observatieperioden overdag (ochtend, middag en namiddag), en in drie opeenvolgende technische replica's. Het equivalent continu A-gewogen geluidsdrukniveau (LAeq) was de primaire uitkomstmaat, en het maximale A-gewogen geluidsdrukniveau (LAFmax) en het piek C-gewogen geluidsdrukniveau (LCpeak) werden geregistreerd als secundaire akoestische uitkomsten. Operationele variabelen werden gelijktijdig vastgelegd, en alle metingen ondergingen kalibratie en vooraf gedefinieerde kwaliteitscontroleprocedures voorafgaand aan de statistische analyse. LAeq, equivalent continu A-gewogen geluidsdrukniveau; LAFmax, maximaal A-gewogen geluidsdrukniveau; LCpeak, piek C-gewogen geluidsdrukniveau. Klik hier om een grotere versie van deze figuur te bekijken.
| Afdeling | Bedcapaciteit | Gemiddeld aantal bezette bedden | Kameroppervlakte (m²) | Kamervolume (m³) | Gemiddeld aantal medewerkers | Gemiddeld aantal bezoekers | Gemiddeld aantal alarmen | Gemiddeld aantal gesprekken |
| Traumatologische orthopedie | 40 | 33.8 | 280.1 | 896.3 | 5.5 | 4.8 | 2.4 | 6.5 |
| Centrum voor Borst- en Schildklierdiagnostiek en -behandeling | 60 | 50.6 | 416.8 | 1333.8 | 5.7 | 4.8 | 2.4 | 6.2 |
| Centrum voor Aangeboren Hartafwijkingen | 35 | 29.4 | 243.6 | 779.5 | 5.7 | 5 | 2.4 | 6.5 |
| Keel- Neus- en Oorheelkunde en Hoofd-Halschirurgie | 50 | 42.2 | 348.1 | 1113.9 | 5.8 | 5 | 2.5 | 6.7 |
| Afdeling Gastro-enterologische Diagnostiek en Behandeling I | 38 | 31.9 | 264.8 | 847.3 | 5.8 | 5 | 2.5 | 6.6 |
| Afdeling Gastroenterologische Diagnostiek en Behandeling II | 52 | 43.9 | 361.7 | 1157.4 | 5.7 | 4.9 | 2.4 | 6.5 |
| Hepatobiliaire en pancreaschirurgie | 48 | 40.5 | 334.1 | 1069.2 | 5.8 | 5.1 | 2.5 | 6.8 |
| Afdeling Intensive Care II | 21 | 17.8 | 145.8 | 466.5 | 5.7 | 4.8 | 2.5 | 6.5 |
| Afdeling Kniegewricht II | 45 | 37.8 | 312.7 | 1000.6 | 5.8 | 4.9 | 2.5 | 6.5 |
| Afdeling ziekte van Parkinson | 55 | 46.2 | 382.2 | 1223 | 5.9 | 5.1 | 2.5 | 6.8 |
| Afdeling Psychosomatische Geneeskunde I | 68 | 57.1 | 472.6 | 1512.3 | 5.8 | 5 | 2.5 | 6.7 |
| Afdeling Pediatrie | 40 | 33.6 | 278.6 | 891.5 | 5.8 | 5.1 | 2.5 | 6.8 |
| Afdeling Urologie | 62 | 52.1 | 431.2 | 1380 | 5.8 | 5 | 2.5 | 6.6 |
| Afdeling Vasculaire Zaken en Neurologie I | 72 | 60.6 | 500.3 | 1601 | 5.9 | 5.1 | 2.5 | 6.9 |
| Afdeling Vasculaire Zaken en Neurologie II | 60 | 50.4 | 417.5 | 1336 | 5.7 | 4.8 | 2.4 | 6.3 |
Tabel 1: Kenmerken van de ziekenhuisafdelingen die in de studie zijn opgenomen. Basiskenmerken van de bedrijfsvoering van de 15 klinische afdelingen die zijn meegenomen in de enquête naar omgevingslawaai. Variabelen op afdelingsniveau omvatten de bedcapaciteit, het gemiddeld aantal bezette bedden, het kameroppervlak, het kamervolume en het gemiddelde aantal medewerkers, bezoekers, alarmen en gesprekken dat is geregistreerd tijdens gestandaardiseerde observatieperioden overdag. Operationele variabelen vertegenwoordigen gemiddelde waarden die tijdens de vooraf gedefinieerde observatieperioden zijn verkregen en werden geëvalueerd als verkennende variabelen in de mixed-effects regressieanalyses. Kameroppervlak en kamervolume worden respectievelijk gerapporteerd in vierkante meters (m2) en kubieke meters (m3).
Beschrijvende akoestische kenmerken over functionele zones
Over alle 1.620 waarnemingen vertoonden de verpleegposten de hoogste omgevingsgeluidsniveaus, gevolgd door gangen, patiëntenactiviteitsruimten en zalen (Tabel 2; Figuur 2 en 3). De gemiddelde LAeq-waarden waren 66,90 ± 2,34 dB(A) bij de verpleegposten, 63,97 ± 2,38 dB(A) in de gangen, 61,38 ± 2,17 dB(A) in de patiëntenactiviteitsruimten en 60,38 ± 2,22 dB(A) in de zalen. Vergelijkbare ruimtelijke gradiënten werden waargenomen voor LAFmax en LCpeak, wat wijst op consistente verschillen over alle drie de akoestische meetwaarden. Gedetailleerde beschrijvende samenvattingen van de secundaire akoestische uitkomsten (LAFmax en LCpeak) worden gegeven in Aanvullende Tabel 1.
| Functionele zone | Observaties (n) | LAeq, dB(A) Gemiddelde ± SD | LAFmax, dB(A) Gemiddelde ± SD | LCpeak, dB(C) Gemiddelde ± SD | Mediaan LAeq (IQR), dB(A) |
| Afdeling | 405 | 60.38 ± 2.22 | 67.95 ± 2.37 | 79.21 ± 3.31 | 60.4 (58.9–61.9) |
| Gang | 405 | 63.97 ± 2.38 | 71.46 ± 2.59 | 82.76 ± 3.56 | 64.0 (62.5–65.5) |
| Verpleegpost | 405 | 66.90 ± 2.34 | 74.43 ± 2.56 | 85.72 ± 3.56 | 66.7 (65.4–68.3) |
| Activiteitsruimte patiënten | 405 | 61.38 ± 2.17 | 69.00 ± 2.42 | 80.20 ± 3.37 | 61.5 (60.0–62.7) |
Tabel 2: Beschrijvende akoestische kenmerken per functionele zone van het ziekenhuis. Beschrijvende akoestische kenmerken voor de vier gestandaardiseerde functionele zones van het ziekenhuis. Waarden worden weergegeven als gemiddelde ± standaarddeviatie (SD), tenzij anders aangegeven. De mediaan en het interkwartielbereik (IQR) worden gerapporteerd voor het equivalent continu A-gewogen geluidsdrukniveau (LAeq). Elke functionele zone droeg 405 waarnemingen bij aan de hiërarchische dataset (totaal n = 1.620). LAeq, equivalent continu A-gewogen geluidsdrukniveau; LAFmax, maximaal A-gewogen geluidsdrukniveau; LCpeak, piek C-gewogen geluidsdrukniveau; IQR, interkwartielbereik.

Figuur 2: Heatmap van gemiddelde equivalente continue A-gewogen geluidsdrukniveaus over ziekenhuisafdelingen en functionele zones. Heatmap die de gemiddelde equivalente continue A-gewogen geluidsdrukniveaus (LAeq) weergeeft voor elke ziekenhuisafdeling en functionele zone. Celwaarden vertegenwoordigen de gemiddelde LAeq in A-gewogen decibels [dB(A)], en de kleurenschaal geeft toenemende geluidsniveaus aan van lagere naar hogere waarden. Functionele zones omvatten zalen, gangen, verpleegposten en patiëntactiviteitsgebieden. CHD, aangeboren hartafwijking; ICU, intensive care unit; LAeq, equivalent continu A-gewogen geluidsdrukniveau. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: Verdeling van de equivalent-continue A-gewogen geluidsdrukniveaus over gestandaardiseerde functionele zones in het ziekenhuis. Boxplots die de verdeling van equivalent-continue A-gewogen geluidsdrukniveaus (LAeq) over vier gestandaardiseerde functionele zones in het ziekenhuis tonen. De boxen vertegenwoordigen het interkwartielbereik (25e–75e percentiel), de horizontale lijn in elke box geeft de mediaan aan, en het ×-symbool duidt het gemiddelde aan. De whiskers reiken tot 1,5 keer het interkwartielbereik (IQR), en waarnemingen buiten de whiskers worden weergegeven als uitschieters. Elke functionele zone bevat 405 waarnemingen (totaal n = 1.620). Verschillen tussen functionele zones werden geëvalueerd met het primaire lineaire mixed-effects model met Tukey-gecorrigeerde paarsgewijze vergelijkingen. LAeq, equivalent-continu A-gewogen geluidsdrukniveau; IQR, interkwartielbereik. Klik hier om een grotere versie van deze figuur te bekijken.
De heat map in Figuur 2 laat een consistent ruimtelijk patroon zien over de 15 deelnemende afdelingen. Verpleegstations vertoonden consistent de hoogste gemiddelde LAeq-waarden, ongeacht het klinische specialisme, terwijl ziekenzaalgebieden de laagste equivalent continue A-gewogen geluidsdrukniveaus vertoonden. Afdelingen met een grotere bedcapaciteit, waaronder Afdeling Vasculaire en Zenuwaandoeningen I, Afdeling Psychosomatische Geneeskunde I en de Afdeling Urologie, vertoonden over het algemeen hogere gemiddelde LAeq-waarden over de functionele zones dan afdelingen met een kleinere bedcapaciteit, zoals Afdeling Intensive Care II en het Centrum voor Aangeboren Hartziekten. Hoewel er een bescheiden interdepartementale variabiliteit werd waargenomen, bleef de relatieve rangschikking van de functionele zones consistent. Gemiddelde LAeq-waarden bij verpleegstations varieerden van 64,1 dB(A) op Afdeling Intensive Care II tot 68,9 dB(A) op Afdeling Vasculaire en Zenuwaandoeningen I, terwijl metingen in de ziekenzalen varieerden van 57,9 tot 61,6 dB(A). Patiëntenactiviteitsgebieden en gangen vertoonden in bijna alle afdelingen tussenliggende waarden (Figuur 2). Figuur 3 illustreert de verdeling van de LAeq-metingen binnen elke functionele zone. De mediane LAeq-waarden waren het hoogst bij de verpleegstations (66,7 dB[A]), gevolgd door de gangen (64,0 dB[A]), patiëntenactiviteitsgebieden (61,5 dB[A]) en ziekenzalen (60,4 dB[A]). Hoewel er in alle functionele zones incidenteel waarnemingen met hoge waarden werden geïdentificeerd, waren de interkwartielafstanden relatief smal, wat duidt op een beperkte variabiliteit binnen elke functionele omgeving. De boxplots tonen een duidelijke scheiding in de verdeling van de LAeq-waarden tussen de vier functionele zones (Figuur 3).
Mixed-effects vergelijking van functionele zones
Om rekening te houden met de hiërarchische structuur van herhaalde akoestische metingen binnen ziekenhuisafdelingen, maakte de primaire analyse gebruik van een lineair mixed-effects model met de ziekenhuisafdeling als random intercept en de functionele zone als fixed effect (Tabel 3). Deze modelleringsaanpak hield rekening met clustering binnen afdelingen, terwijl de gecorrigeerde verschillen in LAeq tussen de vier vooraf gedefinieerde functionele zones werden geschat.
| Vaste effecten (Uitkomst: LAeq) |
| Vast effect | β | SE | 95% BI | p-waarde |
| Intercept (Corridor) | 63.972 | 0.266 | 63,451 tot 64,493 | <0.001 |
| Verpleegpost | 2.931 | 0.146 | 2,644 tot 3,218 | <0.001 |
| Activiteitszone voor patiënten | −2.594 | 0.146 | −2,881 tot −2,307 | <0.001 |
| Afdeling | −3.589 | 0.146 | −3,876 tot −3,302 | <0.001 |
| Geschatte Marginale Gemiddelden |
| Functionele zone | Geschatte gemiddelde LAeq, dB(A) |
| Verpleegpost | 66.9 |
| Gang | 63.97 |
| Activiteitszone voor patiënten | 61.38 |
| Afdeling | 60.38 |
| Tukey-gecorrigeerde paarsgewijze vergelijkingen |
| Vergelijking | Gemiddeld verschil, dB(A) | Gecorrigeerde P-waarde |
| Verpleegpost vs. Gang | 2.93 | <0.001 |
| Verpleegpost versus patiëntenactiviteitengebied | 5.52 | <0.001 |
| Zusterpost versus afdeling | 6.52 | <0.001 |
| Gang versus patiëntenactiviteitsruimte | 2.59 | <0.001 |
| Gang versus Afdeling | 3.59 | <0.001 |
| Activiteitsruimte voor patiënten versus afdeling | 0.99 | <0.001 |
| Modelprestaties |
| Statistiek | Waarde |
| Intraclasscorrelatiecoëfficiënt (ICC) | 0.172 |
| Marginale R² | 0.53 |
| Conditionele R² | 0.61 |
Tabel 3: Primair lineair mixed-effects model ter vergelijking van omgevingsgeluid over verschillende functionele ziekenhuiszones. Resultaten van het primaire lineaire mixed-effects model dat verschillen in equivalente continue A-gewogen geluidsdrukniveaus (LAeq) evalueert tussen de vier gestandaardiseerde functionele ziekenhuiszones. De ziekenhuisafdeling werd gespecificeerd als een random intercept om rekening te houden met de clustering van herhaalde waarnemingen, en de functionele zone werd opgenomen als een fixed effect, waarbij de gang als referentiecategorie diende. Regressiecoëfficiënten (β), standaardfouten (SE), Wald 95% betrouwbaarheidsintervallen (CI) en tweezijdige P-waarden worden gepresenteerd voor de fixed effects. Geschatte marginale gemiddelden werden vergeleken met behulp van Tukey-gecorrigeerde paarsgewijze vergelijkingen. De modelprestaties zijn samengevat met de intraclass correlatiecoëfficiënt (ICC) en de marginale en conditionele determinatiecoëfficiënten (R2), berekend volgens de methode van Nakagawa en Schielzeth. LAeq, equivalent continu A-gewogen geluidsdrukniveau; β, regressiecoëfficiënt; SE, standaardfout; CI, betrouwbaarheidsinterval; ICC, intraclass correlatiecoëfficiënt.
De mixed-effects analyse toonde een zeer significant algemeen effect van de functionele zone op de omgevingsgeluidsniveaus aan (algemene likelihood-ratio test, P < 0,001). Met de gang als referentiecategorie vertoonden posten van verpleegkundigen significant hogere LAeq-waarden (β = 2,931 dB[A], SE = 0,146, 95% BI = 2,644–3,218, P < 0,001), terwijl patiëntenactiviteitsruimtes (β = −2,594 dB[A], SE = 0,146, 95% BI = −2,881 tot −2,307, P < 0,001) en ziekenzalen (β = −3,589 dB[A], SE = 0,146, 95% BI = −3,876 tot −3,302, P < 0,001) significant lagere LAeq-waarden vertoonden (Tabel 3). De geschatte marginale gemiddelden afgeleid uit het mixed-effects model toonden een duidelijk hiërarchisch patroon van omgevingsgeluid over de functionele zones. De gecorrigeerde gemiddelde LAeq was het hoogst bij posten van verpleegkundigen (66,90 dB[A]), gevolgd door gangen (63,97 dB[A]), patiëntenactiviteitsruimtes (61,38 dB[A]) en ziekenzalen (60,38 dB[A]). Alle paarsgewijze vergelijkingen bleven statistisch significant na Tukey-correctie (alle gecorrigeerde P < 0,001), wat wijst op statistisch significante verschillen in LAeq tussen alle vier de functionele zones. De gemiddelde verschillen varieerden van 0,99 dB(A) tussen patiëntenactiviteitsruimtes en ziekenzalen tot 6,52 dB(A) tussen posten van verpleegkundigen en ziekenzalen. Figuur 4 vat de model-gecorrigeerde schattingen samen. Geschatte marginale gemiddelden en hun 95% betrouwbaarheidsintervallen vertoonden minimale overlap tussen de functionele zones, en de Tukey-groepering wees elke zone toe aan een unieke significantiegroep (A–D). Posten van verpleegkundigen werden toegewezen aan groep A met de hoogste gecorrigeerde LAeq, gevolgd door gangen (groep B), patiëntenactiviteitsruimtes (groep C) en ziekenzalen (groep D). Variantiecomponentanalyse toonde aan dat ongeveer 17,2% van de totale variabiliteit in LAeq toe te schrijven was aan verschillen tussen ziekenhuisafdelingen (intraclass correlatiecoëfficiënt = 0,172), wat de inclusie van een random intercept op afdelingsniveau ondersteunt. De fixed effects verklaarden 53% van de variabiliteit in omgevingsgeluid (marginale R2 = 0,53), terwijl het gecombineerde fixed- en random-effects model 61% van de totale variabiliteit verklaarde (conditionele R2 = 0,61), wat duidt op een goede verklarende prestatie van het hiërarchische model (Tabel 3).

Figuur 4: Geschatte Marginale Gemiddelden van Equivalent Continue A-Gewogen Geluidsdrukken in Functionele Zones van het Ziekenhuis. Geschatte marginale gemiddelden (EMMs) van equivalent continue A-gewogen geluidsdrukken (LAeq) in de vier functionele zones van het ziekenhuis, afgeleid uit het primaire lineaire mixed-effects model. De punten representeren de door het model geschatte marginale gemiddelden, en de horizontale foutbalken geven 95% betrouwbaarheidsintervallen (CIs) aan. De ziekenhuisafdeling werd opgenomen als random effect en de functionele zone als fixed effect. Paarsgewijze vergelijkingen werden uitgevoerd met behulp van Tukey-gecorrigeerde geschatte marginale gemiddelden; functionele zones met verschillende letters representeren statistisch significante verschillen (P < 0,05). Elke functionele zone droeg 405 waarnemingen bij aan de analyse. LAeq, equivalent continue A-gewogen geluidsdruk; EMM, geschat marginaal gemiddelde; CI, betrouwbaarheidsinterval. Klik hier om een grotere versie van deze figuur te bekijken.
Parsimoniose mixed-effects regressie van operationele kenmerken geassocieerd met LAeq
Om de operationele kenmerken die geassocieerd zijn met omgevingslawaai te onderzoeken en tegelijkertijd multicollineariteit te minimaliseren, werd een vooraf gespecificeerd parsimonious mixed-effects model toegepast, waarbij functionele zone, bezoekersaantal en bedcapaciteit van de afdeling waren opgenomen als vaste effecten en de ziekenhuisafdeling werd behouden als random intercept (Tabel 4). De selectie van deze variabelen was gebaseerd op de verkennende correlatiematrix (Aanvullende Figuur 1), waaruit sterke correlaties bleken tussen verschillende structurele variabelen, in het bijzonder kameroppervlak en kamervolume (Spearman’s ρ = 0,94). Bijgevolg werden sterk gecorreleerde structurele variabelen uitgesloten van het parsimonious model om de interpreteerbaarheid en modelstabiliteit te verbeteren. Na correctie voor de operationele covariabelen bleef de functionele zone de sterkste onafhankelijke factor die geassocieerd was met omgevingslawaai. Ten opzichte van gangen waren verpleegposten geassocieerd met een gecorrigeerde toename van 2,933 dB(A) (SE = 0,146, 95% BI = 2,647–3,220, P < 0,001), terwijl patiëntenactiviteitsgebieden (β = −2,602 dB[A], SE = 0,146, 95% BI = −2,888 tot −2,315, P < 0,001) en zalen (β = −3,606 dB[A], SE = 0,146, 95% BI = −3,892 tot −3,319, P < 0,001) significant stiller bleven dan gangen.
| Predictor | β-coëfficiënt | SE | 95% BI | P-waarde |
| Intercept (Gang) | 63.97 | 0.27 | 63.45 tot 64.49 | <0.001 |
| Verpleegpost (vs. Gang) | 2.933 | 0.146 | 2.647 tot 3.220 | <0.001 |
| Activiteitsruimte patiënt (vs. Gang) | −2.602 | 0.146 | −2.888 tot −2.315 | <0.001 |
| Zaal (vs. Gang) | −3.606 | 0.146 | −3.892 tot −3.319 | <0.001 |
| Aantal bezoekers (per extra bezoeker) | 0.06 | 0.026 | 0.008 tot 0.112 | 0.024 |
| Bedcapaciteit (per extra bed) | 0.064 | 0.005 | 0.055 tot 0.073 | <0.001 |
| Model: LAeq ~ Functionele Zone + Aantal bezoekers + Bedcapaciteit + (1 | Afdeling) |
Tabel 4: Parsimoniouse lineaire mixed-effects regressie van operationele kenmerken geassocieerd met omgevingslawaai.Resultaten van het vooraf gespecificeerde parsimoniouse lineaire mixed-effects regressiemodel dat de associaties evalueert tussen het equivalent continu A-gewogen geluidsdrukniveau (LAeq) en geselecteerde operationele kenmerken. Het model bevatte de functionele zone, het aantal bezoekers en de bedcapaciteit van de afdeling als vaste effecten en de ziekenhuisafdeling als een random intercept om rekening te houden met de clustering van herhaalde waarnemingen (model: LAeq ~ Functional Zone + Visitor Count + Bed Capacity + (1 | Department)). De gang werd gespecificeerd als de referentiecategorie voor de functionele zone, en het aantal bezoekers en de bedcapaciteit werden gemodelleerd als continue variabelen. Regressiecoëfficiënten (β), standaardfouten (SE), Wald 95% betrouwbaarheidsintervallen (CI) en tweezijdige P-waarden worden gerapporteerd. Regressiecoëfficiënten voor het aantal bezoekers en de bedcapaciteit vertegenwoordigen de geschatte verandering in LAeq geassocieerd met respectievelijk elke extra bezoeker en elk extra ziekenhuisbed. Het intercept vertegenwoordigt de geschatte gemiddelde LAeq voor de referentiecategorie (gang) wanneer de continue predictoren gelijk zijn aan nul. Deze tabel presenteert het vooraf gespecificeerde parsimoniouse model dat is ontwikkeld om overfitting van het model te minimaliseren; het volledig aangepaste sensitiviteitsmodel inclusief alle geregistreerde operationele kenmerken wordt gepresenteerd in Aanvullende Tabel S2. LAeq, equivalent continu A-gewogen geluidsdrukniveau; β, regressiecoëfficiënt; SE, standaardfout; CI, betrouwbaarheidsinterval.
Onder de continue operationele predictoren vertoonde het aantal bezoekers een onafhankelijk positief verband met LAeq. Elke extra bezoeker die aanwezig was tijdens de meting was geassocieerd met een geschatte toename van 0,060 dB(A) (SE = 0,026, 95% CI = 0,008–0,112, P = 0,024). Op dezelfde manier bleef de bedcapaciteit van de afdeling onafhankelijk geassocieerd met omgevingsgeluid, waarbij elk extra ziekenhuisbed geassocieerd was met een geschatte toename van 0,064 dB(A) in LAeq (SE = 0,005, 95% CI = 0,055–0,073, P < 0,001).
Figuur 5 illustreert de gefittede voorspellingen gegenereerd vanuit het parsimonieuze mixed-effects model. De voorspelde LAeq nam ongeveer lineair toe over het waargenomen bereik van de bedcapaciteit per afdeling, en hogere bezoekersaantallen verschoven de voorspelde geluidsniveaus consequent naar boven over het gehele capaciteitsbereik. De parallelle gefittede regressielijnen duiden op additieve effecten van het bezoekersaantal en de bedcapaciteit binnen het gefitte model na correctie voor functionele zone en clustering op afdelingsniveau. Modelgebaseerde voorspellingen toonden aan dat voor een afdeling met ongeveer 20 bedden de voorspelde LAeq varieerde van 61,6 dB(A) bij een laag bezoekersaantal (drie bezoekers) tot 62,4 dB(A) bij een hoger bezoekersaantal (zeven bezoekers). Aan het bovenste uiteinde van het waargenomen bereik van de bedcapaciteit (72 bedden) steeg de overeenkomstige voorspelde LAeq respectievelijk naar 64,5 dB(A) en 65,3 dB(A), wat de gecombineerde bijdrage van de bedcapaciteit per afdeling en het bezoekersaantal binnen het gefitte model illustreert.

Figuur 5: Modelgebaseerde voorspellingen van equivalente continue A-gewogen geluidsdrukniveaus op basis van de bedcapaciteit van de afdeling. Voorspelde equivalente continue A-gewogen geluidsdrukniveaus (LAeq) over het waargenomen bereik van de bedcapaciteit van de afdeling, afgeleid van het parsimone lineaire mixed-effects regressiemodel. De gekleurde lijnen vertegenwoordigen de door het model voorspelde LAeq-waarden bij een laag (3 bezoekers), gemiddeld (5 bezoekers) en hoog (7 bezoekers) aantal bezoekers, terwijl alle andere modeltermen constant worden gehouden. De gearceerde gebieden geven de 95% betrouwbaarheidsintervallen aan voor de voorspelde waarden, en de punten vertegenwoordigen de waargenomen metingen. Voorspellingen worden alleen getoond binnen het waargenomen bereik van de bedcapaciteit van de afdeling die in de studie is opgenomen. LAeq, equivalent continu A-gewogen geluidsdrukniveau; CI, betrouwbaarheidsinterval. Klik hier om een grotere versie van deze figuur te bekijken.
Sensitiviteitsanalyses met behulp van het volledig gecorrigeerde mixed-effectsmodel (Aanvullende Tabel 2), waarin daarnaast rekening werd gehouden met bezette bedden, kameroppervlakte, kamervolume, aantal personeelsleden, alarmgebeurtenissen, gespreksgebeurtenissen en meetperiode, leverden zeer consistente schattingen op voor de functionele zone en het aantal bezoekers. Verpleegstations bleven de luidruchtigste functionele zone, en het aantal bezoekers bleef onafhankelijk geassocieerd met een hogere LAeq (β = 0,060, P = 0,023), terwijl bezette bedden, kameroppervlakte, kamervolume, aantal personeelsleden, alarmgebeurtenissen, gespreksgebeurtenissen en de meetperiode na wederzijdse correctie niet onafhankelijk geassocieerd waren met omgevingsruis. Deze bevindingen ondersteunen de robuustheid van het parsimonieuze model en geven aan dat de functionele zone en het aantal bezoekers de belangrijkste variabelen waren die in deze studie geassocieerd waren met omgevingsruis.
Beschrijvende vergelijking met referentieakoestische criteria
Om de waargenomen omgevingsgeluidsniveaus in context te plaatsen, werden de gemeten LAeq-waarden beschrijvend vergeleken met drie veel geciteerde akoestische referentiecriteria: de WHO-richtlijnwaarde overdag voor patiëntzorgomgevingen in ziekenhuizen (35 dB[A]), de EN ISO 16032-referentiewaarde voor algemene binnenomgevingen (40 dB[A]) en de Chinese milieukwaliteitsnorm GB3096-2022 Klasse 1-referentiewaarde (45 dB[A]). Omdat deze referentiewaarden verschillen in reikwijdte, beoogde toepassing, meetmethodologie en middelingstermijn, werden de vergelijkingen uitsluitend voor beschrijvende doeleinden uitgevoerd en vormen zij geen formele beoordelingen van de naleving van regelgeving (Figuur 6). Over alle 1.620 akoestische waarnemingen heen overschreden de omgevingsgeluidsniveaus elk van de geselecteerde referentiewaarden. In totaal overschreden 1.615 waarnemingen (99,7%) de WHO-richtlijnwaarde overdag van 35 dB(A), overschreden 1.586 waarnemingen (97,9%) de EN ISO 16032-referentiewaarde van 40 dB(A) en overschreden 1.513 waarnemingen (93,2%) de Chinese GB3096-2022 Klasse 1-referentiewaarde van 45 dB(A).

Figuur 6: Percentage van akoestische waarnemingen dat geselecteerde internationale en nationale akoestische referentiecriteria overschrijdt in verschillende functionele ziekenhuiszones. Percentage van de waarnemingen van het equivalent continu A-gewogen geluidsdrukniveau (LAeq) dat geselecteerde internationale en nationale akoestische referentiecriteria overschrijdt in de vier gestandaardiseerde functionele ziekenhuiszones. De staven representeren het aandeel waarnemingen dat de referentiedrempels van de richtlijn van de Wereldgezondheidsorganisatie (WHO) (>35 dB[A]), EN ISO 16032 (>40 dB[A]) en de Chinese nationale standaard GB3096-2022 (>45 dB[A]) overschrijdt. Percentages werden berekend als het aantal waarnemingen dat elke referentiedrempel overschreed, gedeeld door het totale aantal waarnemingen binnen elke functionele zone (n = 405). Omdat deze referentiedocumenten verschillen in hun beoogde toepassing, meetmethodologie, omgevingsinstelling en middelingstermijn, worden de vergelijkingen uitsluitend voor beschrijvende doeleinden gepresenteerd en vormen zij geen formele nalevingsbeoordelingen. LAeq, equivalent continu A-gewogen geluidsdrukniveau; WHO, Wereldgezondheidsorganisatie. Klik hier om een grotere versie van deze figuur te bekijken.
Er werden duidelijke verschillen waargenomen tussen de functionele zones. Verpleegposten vertoonden het hoogste aandeel waarnemingen dat elke referentiewaarde overschreed, waarbij 100,0% zowel de WHO- als de EN ISO-referentiewaarden overschreed en 99,0% de Chinese referentiewaarde overschreed. Gangen vertoonden eveneens hoge overschrijdingspercentages ( respectievelijk 100,0%, 99,3% en 96,0%). In vergelijking hiermee vertoonden ziekenzalen de laagste — maar nog steeds aanzienlijke — aandelen waarnemingen die de drie referentiewaarden overschreden (respectievelijk 99,3%, 95,3% en 87,9%), terwijl patiëntenactiviteitsgebieden intermediaire aandelen vertoonden (respectievelijk 99,5%, 97,0% en 90,6%). Deze beschrijvende bevindingen waren consistent met de mixed-effects analyses, waarbij verpleegposten de hoogste omgevingsgeluidsniveaus vertoonden en ziekenzalen de laagste van de vier functionele zones. De hoge aandelen waarnemingen die alle drie de referentiewaarden overschreden, geven aan dat verhoogd omgevingsgeluid overdag in het gehele ziekenhuis werd waargenomen en niet beperkt bleef tot individuele afdelingen of functionele zones.
Gevoeligheidsanalyses en modeldiagnostiek
De robuustheid van de belangrijkste bevindingen werd geëvalueerd met behulp van een volledig aangepast sensitiviteitsmodel met gemengde effecten, waarin daarnaast rekening werd gehouden met bezette bedden, kameroppervlakte, kamervolume, personeelsaantal, alarmgebeurtenissen, gesprekgebeurtenissen en de meetperiode (Aanvullende Tabel 2). De schattingen voor de functionele zones bleven zeer consistent met die waargenomen in het vooraf gespecificeerde parsimone model. Verpleegposten vertoonden consistent een significant hogere LAeq dan gangen, terwijl patiëntactiviteitsgebieden en zalen significant stiller bleven (allemaal P < 0,001). Het aantal bezoekers bleef ook onafhankelijk geassocieerd met een hogere LAeq (β = 0,060, P = 0,023), terwijl bezette bedden, kameroppervlakte, kamervolume, personeelsaantal, alarmgebeurtenissen, gesprekgebeurtenissen en de meetperiode na simultane correctie niet onafhankelijk geassocieerd waren met omgevingsgeluid (Aanvullende Tabel 2). Deze resultaten geven aan dat de hoofdbevindingen robuust waren voor aanvullende correctie voor gecorreleerde operationele kenmerken. Modelveronderstellingen werden geëvalueerd met behulp van residuele diagnostische procedures (Aanvullende Figuur 1). De normal quantile-quantile (Q–Q) plot toonde een nauwe overeenstemming tussen de waargenomen en theoretische residuele kwantielen, wat duidt op geen substantiële afwijking van normaliteit. De residu-versus-fitted plot toonde een willekeurige spreiding rond nul zonder waarneembare trends of bewijs van heteroscedasticiteit, wat de aannames van lineariteit en constante residuele variantie ondersteunt. In overeenstemming met deze visuele bevindingen toonde de Shapiro–Wilk test geen statistisch significante afwijking van normaliteit aan (W = 0,992, P = 0,084). Aanvullende statistieken voor modelprestaties zijn samengevat in Aanvullende Tabel 3. De residuele standaarddeviatie was 1,92 dB(A) en de standaarddeviatie van het random-intercept op afdelingsniveau was 0,90 dB(A). De intraclass correlatiecoëfficiënt was 0,172, wat duidt op een matige clustering van waarnemingen binnen ziekenhuisafdelingen. Het model met gemengde effecten verklaarde 53% van de variabiliteit door de vaste effecten alleen (marginale R2 = 0,53) en 61% wanneer de random effects op afdelingsniveau werden meegenomen (conditionele R2 = 0,61). De variantie-inflatiefactoren waren uniform laag (maximale VIF = 2,14), wat aangeeft dat er geen bewijs is voor problematische multicollineariteit tussen de behouden predictoren. Bovendien waren alle Cook’s distance-waarden <1,0 (maximum = 0,18), wat suggereert dat geen enkele individuele waarneming een ongebruikelijke invloed uitoefende op het fitted model. De verkennende correlatiematrix (Aanvullende Figuur 2) toonde de verwachte positieve associaties tussen verschillende operationele kenmerken. Sterke correlaties werden waargenomen tussen kameroppervlakte en kamervolume (Spearman’s ρ = 0,94), bedcapaciteit en bezette bedden (ρ = 0,86), en bedcapaciteit en kamervolume (ρ = 0,81). Matige correlaties werden waargenomen tussen het aantal bezoekers en bezette bedden (ρ = 0,71) en tussen het personeelsaantal en gesprekgebeurtenissen (ρ = 0,62). Deze bevindingen ondersteunden de vooraf gespecificeerde variabeleselectiestrategie voor het parsimone model met gemengde effecten, waarbij sterk gecorreleerde structurele variabelen werden uitgesloten om multicollineariteit te minimaliseren en tegelijkertijd de interpreteerbaarheid van het model te behouden.
Over het algemeen vertoonden de primaire mixed-effects analyses, sensitiviteitsanalyses, residualen-diagnostiek en de exploratieve correlatiebeoordeling intern consistente resultaten. In alle analyses vertoonde de functionele zone de sterkste onafhankelijke associatie met omgevingsgeluid, terwijl het aantal bezoekers en de bedcapaciteit van de afdeling ook onafhankelijk geassocieerd bleven met LAeq. De meeste andere geregistreerde operationele kenmerken waren na wederzijdse correctie niet onafhankelijk geassocieerd met omgevingsgeluid.
Beschikbaarheid van gegevens:
De gegevens die de bevindingen van deze studie ondersteunen, zijn beschikbaar in het artikel en de bijbehorende aanvullende materialen. De volledige ruwe dataset van akoestische metingen, kalibratieregisters, de dataset van operationele variabelen, de brongegevens van de figuren en de R-code voor de statistische analyse zijn gedeponeerd in de Zenodo-repository en zijn publiekelijk toegankelijk via https://doi.org/10.5281/zenodo.21404807. Aanvullende tabellen 1–3 en Aanvullende figuren 1 en 2 bieden aanvullende verwerkte datasets, sensitiviteitsanalyses, modeldiagnostiek en ondersteunende analyses.
Aanvullende figuur 1: Residuele diagnostische plots voor het parsimonieuze lineaire mixed-effects model. Diagnostische plots ter evaluatie van de aannames van het parsimonieuze lineaire mixed-effects regressiemodel. (A) Normaal quantile-quantile (Q–Q) plot van de conditionele residuen, die een benaderde normaliteit van de residuverdeling aantoont. (B) Plot van de conditionele residuen versus de gefitte waarden, die geen substantiële afwijking van homoscedasticiteit of systematisch patroon over het gefitte bereik laat zien. Ziekenhuisafdeling werd opgenomen als een random intercept, en functionele zone, bezoekersaantal en bedcapaciteit werden opgenomen als fixed effects. Residuele diagnostiek ondersteunt de adequaatheid van de modelaannames voor de primaire regressieanalyse. Diagnostische plots werden gegenereerd met behulp van het DHARMa-pakket in R. LAeq, equivalent continu A-gewogen geluidsdrukniveau.Klik hier om dit bestand te downloaden.
Aanvullende Figuur 2: Spearman-correlatiematrix van operationele variabelen opgenomen in de exploratieve mixed-effects analyses. Spearman-correlatiematrix die de paarsgewijze associaties toont tussen departementale operationele variabelen die zijn geëvalueerd tijdens exploratieve mixed-effects regressieanalyses. Celwaarden vertegenwoordigen de Spearman-rangcorrelatiecoëfficiënten (ρ), en de kleurenschaal geeft de sterkte en richting van de correlatie aan, waarbij blauw positieve correlaties vertegenwoordigt en rood negatieve correlaties. De correlatieanalyse werd uitgevoerd om potentiële multicollineariteit tussen kandidaat-verklarende variabelen te evalueren vóór de modelontwikkeling. Variabelen die sterke collineariteit vertoonden, werden meegewogen tijdens de modelselectie om multicollineariteit in het parsimonieuze regressiemodel te minimaliseren. Operationele variabelen omvatten de bedcapaciteit van de afdeling, bezette bedden, kameroppervlakte, kamervolume, personeelsbezetting, aantal bezoekers, aantal alarmen en aantal gesprekken. Spearman’s ρ, Spearman-rangcorrelatiecoëfficiënt.Klik hier om dit bestand te downloaden.
Aanvullende tabel 1: Secundaire akoestische uitkomsten per functionele ziekenhuiszone. Beschrijvende samenvatting van de secundaire akoestische uitkomsten gemeten over de vier gestandaardiseerde functionele ziekenhuiszones. Waarden worden gepresenteerd als gemiddelde ± standaarddeviatie (SD) en vertegenwoordigen samengevoegde waarnemingen van alle 15 klinische afdelingen. Het equivalent continu A-gewogen geluidsdrukniveau (LAeq) was de primaire akoestische uitkomst, terwijl het maximaal A-gewogen geluidsdrukniveau (LAFmax) en het piek C-gewogen geluidsdrukniveau (LCpeak) werden geregistreerd als secundaire akoestische maten om de maximale en piek geluidsniveaus tijdens de routinematige dagelijkse ziekenhuisactiviteiten te karakteriseren. LAeq, equivalent continu A-gewogen geluidsdrukniveau; LAFmax, maximaal A-gewogen geluidsdrukniveau; LCpeak, piek C-gewogen geluidsdrukniveau; SD, standaarddeviatie.Klik hier om dit bestand te downloaden.
Aanvullende tabel 2: Volledig aangepast sensitiviteitsmodel voor lineaire mixed-effects regressie. Resultaten van het volledig aangepaste sensitiviteitsmodel voor lineaire mixed-effects regressie ter evaluatie van de associaties tussen het equivalent continu A-gewogen geluidsdrukniveau (LAeq) en alle geregistreerde operationele kenmerken. De ziekenhuisafdeling werd opgenomen als een random intercept om rekening te houden met de clustering van herhaalde waarnemingen. Functionele zone, tijdsperiode, bedcapaciteit van de afdeling, bezette bedden, kameroppervlakte, kamervolume, personeelsaantal, bezoekersaantal, aantal alarmen en aantal gesprekken werden opgenomen als fixed effects. De gang diende als referentiecategorie voor de functionele zone, en de middag diende als referentiecategorie voor de tijdsperiode. Regressiecoëfficiënten (β), standaardfouten (SE), Wald 95% betrouwbaarheidsintervallen (CI) en tweezijdige P-waarden worden gerapporteerd. Dit model werd uitgevoerd als sensitiviteitsanalyse om de robuustheid van het parsimone model in Tabel 4 te evalueren na correctie voor alle geregistreerde operationele kenmerken. LAeq, equivalent continu A-gewogen geluidsdrukniveau; β, regressiecoëfficiënt; SE, standaardfout; CI, betrouwbaarheidsinterval.Klik hier om dit bestand te downloaden.
Aanvullende tabel 3: Modelprestaties en diagnostische statistieken voor het primaire lineaire mixed-effects model. Overzicht van de modelprestaties en diagnostische statistieken voor het primaire lineaire mixed-effects model ter evaluatie van equivalente continue A-gewogen geluidsdrukken (LAeq) over verschillende functionele zones in het ziekenhuis. De residuele standaarddeviatie (σ) vertegenwoordigt de onverklaarde variabiliteit binnen het model, terwijl de standaarddeviatie van het random-intercept de variabiliteit tussen afdelingen kwantificeert. De intraclass correlatiecoëfficiënt (ICC) schat het aandeel van de totale variantie toe te schrijven aan clustering binnen ziekenhuisafdelingen. Marginale R2 vertegenwoordigt het aandeel van de variantie dat enkel door de vaste effecten wordt verklaard, terwijl conditionele R2 het aandeel van de variantie vertegenwoordigt dat door zowel de vaste als de random effecten wordt verklaard. Het Akaike-informatiecriterium (AIC), het Bayesiaanse informatiecriterium (BIC) en de log-likelihood zijn opgenomen als maten voor de model fit. De variantie-inflatiefactoren (VIF's) waren allemaal <5, wat duidt op de afwezigheid van problematische multicollineariteit, en de Cook's distance-waarden waren <1,0 voor alle waarnemingen, wat aangeeft dat geen enkele individuele waarneming een onevenredige invloed uitoefende op de schattingen van het model. LAeq, equivalent continu A-gewogen geluidsdrukniveau; ICC, intraclass correlatiecoëfficiënt; VIF, variantie-inflatiefactor.Klik hier om dit bestand te downloaden.