participanți & procedura
participanții au fost extrași din valurile 6 și 7 (2014-2015/2015-2016) ale monitorului angajamentului cu mediul Natural (MENE) sondaj (singurele valuri în care rezultatele noastre cheie au fost măsurate în mod constant). Sondajul, care face parte din statisticile naționale ale guvernului britanic, este repetat transversal (diferite persoane participă la fiecare val) și se desfășoară în întreaga Anglie și pe tot parcursul anului (aprox. 4 000 de persoane pe săptămână) pentru a reduce potențialele tendințe geografice și sezoniere49. Ca parte a statisticilor oficiale din Marea Britanie, protocoalele de eșantionare sunt extinse, pentru a asigura un eșantion cât mai reprezentativ al populației engleze adulte. Detalii complete pot fi găsite în rapoartele tehnice anuale MENE49 cu caracteristici cheie, inclusiv: (a) „un sistem computerizat de eșantionare care integrează fișierul adresei Oficiului Poștal cu datele recensământului din 2001 privind zonele mici la nivelul zonei de ieșire. Acest lucru permite valuri reproduse de eșantioane stratificate în mai multe etape”; (b) „zonele din fiecare regiune Standard sunt stratificate în benzi de densitate a populației și în bandă, în ordine descrescătoare în funcție de procentul populației din gradele socio-economice I și II”; (c) ” maximizați precizia statistică a eșantionării, undele secvențiale de lucru pe teren sunt alocate sistematic în cadrul eșantionării pentru a; (d) „pentru a asigura un eșantion echilibrat de adulți în cadrul adreselor contactate efectiv, se stabilește o cotă în funcție de sex (bărbat, femeie casnică, femeie necasnică); în cadrul cotei de femeie casnică, prezența copiilor și statutul de muncă și în cadrul cotei de sex masculin, statutul de lucru”; și (e) „datele sondajului sunt ponderate pentru a se asigura că eșantionul este reprezentativ pentru populația Regatului Unit în ceea ce privește caracteristicile demografice standard” (ref.49, p. 5). Datele sunt colectate folosind interviuri față-în-față la domiciliu cu răspunsuri înregistrate folosind software-ul de intervievare personală asistată de calculator (CAPI).
deși eșantionul total pentru acești ani a fost n = 91.190, întrebările privind sănătatea și bunăstarea au fost adresate doar în fiecare a patra săptămână de eșantionare (adică lunar, mai degrabă decât săptămânal), rezultând un eșantion redus de n = 20.264. Pentru a ține seama de orice prejudecăți reziduale ale eșantionării la acest nivel lunar, în setul de date sunt incluse ponderi speciale ale sondajului lunar. Acestea au fost aplicate în analiza actuală pentru a se asigura că rezultatele au rămas generalizabile pentru întreaga populație adultă din Anglia. Toate datele au fost anonimizate de Natural England și sunt accesibile publicului la: http://publications.naturalengland.org.uk/publication/2248731?category=47018. Aprobarea etică nu a fost necesară pentru această analiză secundară a statisticilor naționale disponibile public.
rezultate: sănătate auto-raportată & bunăstare subiectivă
sănătatea auto-raportată (de acum înainte: sănătate) a fost evaluată folosind un singur articol: „cum este sănătatea ta în general?'(uneori denumit ‘SF1’). Opțiunile de răspuns au fost: ‘foarte rău’, ‘rău’, ‘corect’, ‘bun’și’ foarte bun’. Răspunsurile sunt puternic asociate cu utilizarea serviciilor medicale50 și mortalitate51; și, în mod esențial, pentru scopurile actuale, zona verde de vecinătate13. În urma lucrărilor anterioare, am dihotomizat răspunsurile în ‘Bun’ (‘Bun/foarte bun’, ponderat = 76,5%) și’nu bun’ (‘corect/rău/foarte rău’, 23,5%)52. Bunăstarea subiectivă (de acum înainte: bunăstarea) a fost evaluată folosind măsura satisfacției de viață, una dintre măsurile naționale de bunăstare din Marea BRITANIE53: ‘în general, cât de mulțumit sunteți de viața din zilele noastre?’cu răspunsuri variind de la 0’ deloc ‘la 10’complet’. Din nou, în urma studiilor anterioare am dihotomizat răspunsurile în ‘mare’ (8-10, 60.2%) și bunăstare scăzută (0-7, 39,8%) 54. Histogramele distribuțiilor (non-normale) pentru ambele variabile de rezultat sunt prezentate în Anexa A. De notat deși punctele de dihotomizare s-au bazat pe cercetări anterioare, acestea sunt în concordanță cu datele actuale; percentila 50 pentru sănătate a fost în răspunsul ‘bun’ și pentru bunăstare în ‘8’. Analizele de sensibilitate efectuate asupra variațiilor ordinale (atât pentru sănătate, cât și pentru bunăstare) și liniare (numai pentru bunăstare) ale acestor variabile sunt prezentate în anexa E.
expunere: Contactul recreativ cu natura în ultimele 7 zile
contactul recreativ cu natura sau timpul petrecut în medii naturale în ultima săptămână a fost derivat prin înmulțirea numărului de vizite recreative raportate pe săptămână cu durata unei vizite selectate aleatoriu în ultima săptămână. Participanții au fost prezentați la sondaj după cum urmează: „vă voi întreba despre ocaziile din ultima săptămână când v-ați petrecut timpul în afara ușilor. Prin out of doors înțelegem spații deschise în și în jurul orașelor și orașelor, inclusiv parcuri, canale și zone naturale; coasta și plajele; și rural, inclusiv terenuri agricole, păduri, dealuri și râuri. Acest lucru ar putea fi orice, de la câteva minute până la toată ziua. Poate include timpul petrecut aproape de casa sau locul de muncă, mai departe sau în timpul vacanței în Anglia. Cu toate acestea, aceasta nu include: excursii de cumpărături de rutină sau; timpul petrecut în propria grădină.”Apoi au fost întrebați „de câte ori, dacă ați făcut acest tip de vizită ieri/pe <ziua> ” pentru fiecare dintre cele șapte zile anterioare. Nouăzeci și opt la sută dintre respondenți au raportat săptămâna trecută vizite la 7. Restul de 2% au fost limitate la 7 vizite pentru a evita distorsionarea dramatică a estimărilor săptămânale.
după înregistrarea detaliilor de bază ale fiecărei vizite (până la 3 pe zi), o singură vizită a fost selectată aleatoriu de către software-ul CAPI, pentru ca intervievatorul să pună întrebări suplimentare, inclusiv: „cât a durat această vizită?”(Ore & minute). Datorită selecției aleatorii, chiar dacă vizita selectată nu a fost neapărat reprezentativă pentru o anumită persoană, procedura de randomizare ar trebui să reducă potențialele prejudecăți la nivelul populației la care au fost efectuate analizele noastre. Estimările săptămânale ale duratei au fost astfel obținute prin înmulțirea duratei acestei vizite selectate aleatoriu cu numărul de vizite declarate în ultimele șapte zile (plafonat la 7). În urma abordării studiilor anterioare de expunere-răspuns în domeniu (de exemplu, Shanahan și colab., 2016), durata a fost clasificată în 7 categorii: 0 minute (n = 11.668) ; 1-59 minute (n = 355); 60-119 minute (n = 1.113); 120-179 minute (n = 1.290); 180-239 minute (n = 1.014); 240-299 minute (n = 882); 300 minute (n = 3.484). O bandă alternativă la 30 de minute a fost problematică din cauza Ns foarte scăzut pentru unele benzi (de exemplu, 1-29 minute, n = 85), reflectând faptul că estimările duratei săptămânale s-au grupat în jurul marcajelor de oră, de exemplu, 78% din observațiile neponderate din banda de 120-179 minute au fost exact 120 de minute (a se vedea anexa a, figura C pentru histograma duratei). Cea mai mare bandă a fost plafonată la 300 de minute de la sută din cauza înclinării mari pozitive a datelor.
variabilele de Control
sănătatea și bunăstarea sunt asociate cu caracteristicile socio-demografice și de mediu atât la nivel de vecinătate (de exemplu, privarea de zonă), cât și la nivel individual (de exemplu, starea relației) 55. Deoarece multe dintre aceste variabile pot fi, de asemenea, legate de expunerea la natură, acestea au fost controlate în analizele ajustate.
variabile de control la nivel de zonă
datele covariate la nivel de zonă au fost atribuite la nivelul spațial al recensământului din 2001 zonele de super-ieșire cu Strat Inferior (LSOAs) în care trăiau indivizii. În Anglia erau 32.482 LSOAs, fiecare conținând aproximativ 1.500 de persoane într-o zonă fizică medie de 4km2.
Neighbourhood greenspace
pentru a înțelege cât de mult este greenspace în vecinătatea unui individ, am obținut o valoare a densității zonei folosind baza de date generalizată de utilizare a terenurilor (GLUD)56. GLUD oferă, pentru fiecare LSOA din Anglia, zona acoperită de spațiu verde și grădini domestice. Acestea au fost însumate și împărțite la suprafața totală lsoa pentru a furniza metrica densității spațiului verde. Această valoare a fost alocată fiecărui individ din eșantion, pe baza lsoa de reședință. În urma literaturii anterioare, indivizii au fost repartizați la una dintre cele cinci chintile de spațiu verde pe baza acestei definiții (variind de la cel mai puțin verde la cel mai verde)33. În loc să derivăm chintile de spațiu verde din eșantionul curent (adică împărțiți eșantionul curent în cinci părți egale pe baza procentului de spațiu verde din LSOA), am atribuit indivizi în schimb uneia dintre cele cinci chintile de spațiu verde predeterminate pe baza distribuției spațiului verde în toate cele 32.482 LSOAs din Anglia. Deși acest lucru a însemnat că nu am obținut exact cote egale de 20% din eșantionul nostru actual în chintilele greenspace (deși datorită protocolului de eșantionare eram încă foarte aproape de acest lucru, a se vedea anexa B) această abordare a permis inferențe să se facă în întreaga țară, mai degrabă decât pur și simplu la eșantionul actual. În analizele de sensibilitate exploratorie am definit greenspace ca fiind doar categoria GLUD ‘greenspace’, cu categoria GLUD’ gardens ‘ exclusă. Acest lucru a produs rezultate foarte similare, așa că ne-am concentrat pe o definiție mai incluzivă, incluzând ambele aspecte. În analizele ulterioare de sensibilitate exploratorie, am atribuit indivizilor cinci categorii de spațiu verde definite prin intervale egale de acoperire a spațiului verde (de exemplu, 0-20%, 21-40%, 41-60% etc.) mai degrabă decât chintile pe baza procentelor populației. Acest lucru a produs, de asemenea, rezultate foarte similare, așa că din nou am decis să mergem cu cea mai comună abordare. În analizele ulterioare, cea mai mică chintilă verde a acționat ca categorie de referință.
deprivarea zonei
fiecare LSOA din Anglia este evaluată în funcție de mai mulți parametri ai deprivării, inclusiv șomajul și criminalitatea, nivelurile educaționale, veniturile, valorile de sănătate, barierele în calea locuințelor și serviciilor și mediul de viață. Un indice total al scorului privării Multiple (IMD) este derivat din aceste subdomenii57. În urma studiilor anterioare52, am repartizat indivizii în chintile de privare pe baza LSOA în care trăiau. Ca și în cazul greenspace, punctele de tăiere pentru cvintilele de privare a zonei s-au bazat, de asemenea, pe toate LSOAs-urile din Anglia, mai degrabă decât pe cele din eșantionul actual, pentru a permite inferența populației în ansamblu (cvintila cea mai lipsită = ref).
poluarea aerului
o măsură indicativă a poluării aerului a fost operaționalizată ca fundal LSOA PM10 atribuit tertilelor tuturor LSOAs din Anglia (cea mai mică concentrație de particule = ref). Concentrațiile de PM10, bazate pe simulările modelului de cartografiere a climei de poluare (PCM) 58, au fost medii în perioada 2002-2012 și agregate de la rezoluția de 1 km pătrați la LSOAs.
controale la nivel Individual
controale la nivel Individual comparabile cu studiile anterioare din acest domeniu6, 7, 12,13, 15 au inclus: sex (bărbat = ref); vârstă (clasificată ca 16-64 = ref; 65+); grad social ocupațional (AB (cea mai mare, de exemplu, managerial), C1, C2 și DE (cea mai mică, de exemplu, muncă necalificată, = ref) ca; statutul ocupării forței de muncă (full-time, part-time, în educație, pensionar, nu lucrează/șomeri = ref); statutul relației (căsătorit/coabitare; singur/separat/divorțat/văduv = ref); etnie (Alb Britanic; altele = ref); numărul de copii din gospodărie (1 vs.0 = ref); și dreptul de proprietate câine (Da; Nu = ref).
două variabile de control suplimentare au fost deosebit de importante. În primul rând, sondajul a întrebat: Aveți vreo boală de lungă durată, o problemă de sănătate sau un handicap care vă limitează activitățile zilnice sau tipul de muncă pe care o puteți face?'(‘Funcționare restricționată’: Da; Nu = ref). Inclusiv această variabilă, cel puțin parțial, controlează cauzalitatea inversă. Dacă se constată asocieri similare între expunerea la natură și sănătate și bunăstare atât pentru cei cu sau fără funcționare restricționată, acest lucru ar susține ideea că asociațiile nu se datorează doar persoanelor mai sănătoase și mai mobile care vizitează natura mai des.
am controlat, de asemenea, pentru numărul de zile pe săptămână de persoane raportate angajarea în activitate fizică > 30 de minute; în analiza actuală dihotomizat ca fie întâlnire sau nu îndeplinesc liniile directoare de 150 de minute pe săptămână (adică. 5 zile în săptămână cu activitate fizică>30 de minute). Unii oameni realizează această orientare deși activitatea fizică în medii naturale35, astfel, orice asociere între timpul petrecut în natură și sănătate se poate datora pur și simplu activității fizice angajate în aceste setări. Credem că acest lucru nu este cazul în contextul actual, deoarece corelația (ordinea rangului) dintre contactul săptămânal cu natura și numărul de zile pe săptămână ale unei persoane angajate în >30 de minute de activitate fizică a fost doar rs = 0,27. Cu toate acestea, prin controlul nivelurilor de activitate săptămânale, relațiile modelate între timpul în natură și sănătate au mai puține prejudecăți din această sursă și, prin urmare, estimări îmbunătățite ale asocierii cu expunerea naturală în sine.
controale temporale
datorită naturii combinate multianuale a datelor, anul / valul a fost, de asemenea, controlat pentru. Analiza preliminară nu a găsit niciun efect al sezonului în care au fost colectate datele, astfel încât aceasta a fost exclusă din analizele finale.
strategia de analiză
regresiile logistice binomiale ponderate în sondaj au fost utilizate pentru a prezice șansele relative ca un individ să aibă o sănătate bună sau o bunăstare ridicată în funcție de expunerea săptămânală a naturii în ceea ce privește categoriile de durată pe săptămână. Modelul fit a fost furnizat de pseudo R2; aici estimarea mai conservatoare Cox și Snell. Variabilele binare rezultate au fost regresate mai întâi în funcție de categoriile de durată a expunerii pentru a testa relațiile directe; modelele ajustate au fost apoi specificate pentru a include variabilele de control individuale și de nivel de zonă. Din cauza lipsei datelor privind nivelul zonei pentru o mică minoritate de participanți (n = 456), eșantioanele noastre de estimare pentru aceste modele ajustate au fost n = 19.808. Analiza preliminară a constatat că proporțiile descriptive ponderate în rândul acestui eșantion de estimare redus diferă doar neglijabil de cele din toate observațiile disponibile în eșantionul mai larg MENE, sugerând că abordarea noastră completă de analiză a cazurilor nu a denaturat reprezentativitatea populației eșantionului de estimare. Eșantionul complet n = 20.264 a fost menținut pentru modelul neajustat pentru a oferi cea mai precisă și ponderată reprezentare a datelor, deoarece reducerea modelelor neajustate la n = 19.808 a produs rezultate practic identice. Deși analizele noastre principale au utilizat categorii de durată ale contactului săptămânal cu natura, o analiză exploratorie a utilizat modele aditive generalizate care încorporează un spline de durată de regresie cubică penalizată ca variabilă continuă (ajustând pentru același set de covariabile). Acest lucru ne-a permis să producem un grafic ‘mai lin’ al datelor. Analizele și trasarea au fost realizate folosind versiunea R 3.4.1, folosind pachete mgcv și visreg59.
pentru a explora generalizabilitatea oricărui model în diferite grupuri socio-demografice, am stratificat a priori analizele pe mai multe zone și covariabile individuale (așa cum sunt definite mai sus) care s-au dovedit a fi importante în studiile anterioare: (a) urbanitate; (b) spațiu verde de Vecinătate; (c) privarea zonei; (d) Sex; (e) vârstă; (f) funcționare restricționată; (g) statut socio-economic Individual (SES); (f) etnie; și (g) activitate fizică. În cazul celor trei predictori multi-categorie (spațiu verde/deprivare, SES individual), au fost derivate clasificări binare pentru analizele stratificate pentru a menține dimensiuni robuste ale eșantioanelor în fiecare categorie. În cazul lsoa greenspace și privarea, divizările binare au fost făcute pe baza punctului de tăiere median pentru toate LSOAs din Anglia; SES a fost dihotomizat prin prăbușirea categoriilor de grade sociale în mod standard, A/B / C1 vs. C2/D / E.