analüüsi artikli abstarct'i: "Road traffic fatalities remain one of the leading causes of death globally, posing a major challenge to achieving the United Nations Sustainable Development Goals (SDGs), particularly SDG 3.6, which calls for a halving of global road traffic deaths and injuries by 2030. Measuring road traffic safety (RTS) is therefore critical for evidence-based policymaking and sustainable development. Traditionally, the road fatality rate has been used as the principal indicator of RTS. However, this measure alone provides a limited view. This study proposes a multi-indicator approach to measuring RTS by integrating other indicators — the road fatality rate, the number of fatalities, the number of fatal accidents, and the total number of road accidents—into a RTS measurement framework. Using path analysis and Estonian daily traffic data from 2010 to 2018, the results reveal that RTS is more strongly influenced by accident severity than by accident frequency. The results demonstrate that the multi-indicator model offers a more comprehensive and accurate measure of RTS compared to the traditional single-indicator approach for sustainable transport systems." Mida see annab transpordiplaneerimises?
- Road Fatality Rate (RFR) = surmad / miljon sõidukikilomeetrit - Total Fatalities (TF) - Fatal Accidents (FA) - All Accidents (AA)
Abstrakti analüüs ja selle praktiline väärtus transpordiplaneerimisele
| Peamine sisuline punkt | Mõju/kasu transpordiplaneerijale |
|---|---|
| 1️⃣ Ühe‑näitaja (tavapärane) lähenemine on ebapiisav – traditsiooniliselt kasutatakse ainult liiklusohutuse hindamiseks surmaintensiivsust (surmade arv / miljoni sõidukiga). | • Selgitab, miks pelgalt “surnud/10⁶ km” ei anna piisavat pilti riskidest ja seetõttu võib juhtuda, et olulisi ohutusprobleeme jääb varjatuks. |
| 2️⃣ Mitme‑näitaja raamistik – kombineerib surmaintensiivsuse, suremuste koguarvu, surmaga lõppenud õnnetuste arvu ja kõigi liiklusõnnetuste (ka mitte‑surmavate) koguarvu. | • Pakub terviklikumat pildi: – Surnu‑sagedus näitab kui tihti juhtub kõige tõsisemad tagajärjed. – Õnnetuste arv aitab mõista, kus on “kõrge risk” kuid madala surmaintensiivsusega (nt kerged õnnetused). |
| 3️⃣ Path‑analüüs näitab, et raskusaste mõjutab turvalisust tugevamalt kui sagedus. | • Planeerijad saavad suunata ressursse sinna, kus “raskuse aste” on kõrge (nt suurema kiirusega maanteed, ristmikud). – Näiteks investeerida paremasse valgustusse või eraldatud sõiduklassidesse. |
| 4️⃣ Andmekasutus – Eesti igapäevased liiklusandmed 2010‑2018 | • Demonstreerib, et olemasolevad statistilised registrid (õnnetuste andmebaasid) on piisavad mitme‑näitaja mudeli rakendamiseks. – Planeerijad saavad kasutada oma riigi/linna registreeritud andmeid ilma uute mõõtmiste vajaduseta. |
| 5️⃣ Süsteemne ja jätkusuutlik lähenemine SDG‑3.6 suunas | • Selge seos: mitme‑näitaja mudel võimaldab jälgida progressi eesmärgi “surmaintensiivsuse pooleks vähendamine” täpsemalt, sest see eristab sagedust ja raskusastet. – Seega on võimalik koostada realistlikud tegevuskavad ning mõõta nende tõhusust. |
| Samm | Kirjeldus | Kasutatav tööriist/tehnika |
|---|---|---|
| 1. Andmete kogumine | • Koguge olemasolevad liiklusõnnetuste andmed (surnud, haavatud, õnnetuse tüüp). • Lisage sõidukite liiklusmahu ja teekonnaandmeid. |
Statistikaamet, politsei/ohutuskeskused; GIS‑baasid. |
| 2. Indikaatorite arvutus | • Arvutage: – Road Fatality Rate (RFR) = surmad / miljon sõidukikilomeetrit. – Total Fatalities (TF). – Fatal Accidents (FA). – All Accidents (AA). |
Excel/ R / Python (pandas, NumPy). |
| 3. Path‑analüüsi mudel | • Looge strukturaalne seosmudel: AA → FA → TF → RFR (või vastupidi) ning hindage kaalu koefitsiente. | lavaan paket R-s; semopy Pythonis. |
| 4. Tulemuste visualiseerimine | • Kaardistage “kuuma” piirkonnad, kus FA ja TF on kõrged, kuid AA suhteliselt madal (kõrge raskusaste). • Loo trendgraafikud ajas. |
GIS‑kaardid; Tableau/PowerBI; R ggplot2. |
| 5. Poliitikate kujundamine | • Sihitud meetmed: kiirusepiirangute muutmine, turvalised ristmikud, eraldatud rattatee‑või jalakäijate alad. • Seada KPI-d (Key Performance Indicators) iga näitaja kohta. |
Põhjalik tegevusplaan; riskianalüüsid (e.g., HAZOP). |
| 6. Järelevalve | • Korduv hindamine: kas meetmed vähendavad FA ja TF, mitte ainult AA. • Raportide koostamine SDG‑3.6 täitmise kohta. |
Aastaaruanded; interaktiivsed dashboards. |
| Küsimus | Vastus, mis tuleneb mitme‑näitaja mudelist |
|---|---|
| Millele keskenduda? | “Raskuse aste” (FA ja TF) on tugevam mõjutaja – seega tuleb prioriteediks muuta teid ja ristmikke, kus surmaga lõppenud õnnetused koonduvad. |
| Kuidas õigustada investeeringuid poliitikutele? | Selge kvantitatiivne tõestus: “Selle projekti elluviimine vähendab fatal accidents 30 % – see tähendab X elupääste ja Y miljonite euro suuruse sotsiaalse kulude kokkuhoiu.” |
| Kuidas seostada transpordipoliitika laiemate SDG‑eesmärkidega? | Indikaatorid on otse seotud: – SDG 3.6 (vähenenud surmad). – SDG 11 (turvalised ja jätkusuutlikud linnad). – SDG 13/15 (kaasaegse kliima‑ ja keskkonnamõju vähendamine, sest vähem traagilisi õnnetusi → vähem erakorralisi teenuseid jne). |
| Kuidas hinnata olemasolevaid sekkumisi? | Võrrelge “enne” ja “pärast” perioodi nii AA kui FA/TF suhtes – näete, kas meetmed vähendasid pigem õnnetuste arvu või nende raskusastet. |
| Kas mudel sobib ka väiksematesse linnades? | Jah – isegi väikesed andmekogud võimaldavad arvutada samad neli näitajat; tulemused annavad kohalikele omavalitsustele konkreetseid prioriteete (nt “sõiduohutusrajoonid”). |
| Piirang | Mõju / lahendus |
|---|---|
| Andmete kvaliteet – õnnetuste registreerimise täpsus võib erineda riigiti. | • Parandada andmekvaliteeti: standardiseeritud raporti vormid, kohustuslik raporteerimine. |
| Ajaline piirang – 2010‑2018 periood ei pruugi kajastada hiljutisi tehnoloogilisi muutusi (nt ADAS-süsteemide laialdane kasutus). | • Uuendada mudelit uute andmete ja muutujadega: automatiseeritud sõidukite telemaatika, nutikad liiklusmõõdikud. |
| Geograafiline ühtsus – Eesti on suhteliselt väike; suuremates riikides võib piirkondadevaheline heterogeensus olla märkimisväärsem. | • Kohandada mudelit alajaotus‑tasemel (piirkonnad, linnaosad). |
| Sotsiaalsed ja demograafilised tegurid – vanusejaotused, sõidukite tüüp jne võivad mõjutada õnnetuste raskusastet. | • Lisada regressioonimudeli kontrollmuutujad (sõiduajalugu, auto tüübid). |
Kokkuvõttes annab artikkel tööriistakasti uue mõõdikute komplekti, mis aitab transpordiplaneerijatel luua tõhusamaid, andmepõhisemaid ja jätkusuutlikumaid liikluslahendusi – just see on vajalik samm teekonnal, et 2030. aastaks saavutada surmaintensiivsuse pooleks vähendamine.
Allpool on samm‑sammuline juhend, kuidas luua ja rakendada path‑analüüsit mudeliks, mis hindab erinevate liiklusturvalisuse näitajate omavahelist mõju. Juhendis eeldatakse, et teil on olemas igapäevased või korduvad andmed 2010–2018 (või muu perioodi) kohta, milles sisalduvad:
| Muutuja | Kirjeldus / allikas |
|---|---|
AA – All Accidents (kõik liiklusõnnetused) |
Koguarv õnnetusi antud ajaperioodil/maastikul. |
FA – Fatal Accidents (surmaga lõppenud õnnetused) |
Õnnetuste arv, mille tagajärjel on vähemalt üks surm. |
TF – Total Fatalities (surnute koguarv) |
Surnute arv samas ajaperioodis/maastikus. |
RFR – Road Fatality Rate (surmaintensiivsus) |
(Surnud / sõidukikilomeetrid) × 10⁶. |
VMT – Vehicle‑Miles Traveled (sõidukite läbitud kilomeetrite kogumaht) |
Sõiduandmed, mis võimaldavad normaliseerida. |
Kontrollmuutujad (X1…Xk) |
Trafikavolume, ilmastikutingimused, teekatte tüüp, kiirusepiirangud jne. |
AA) või nende raskusaste (FA, TF) – mõjutab kõige enam liiklusohutust, mõõdetuna RFR‑na?AA → kõrgem RFR.FA/TF → suurem RFR, isegi kui AA on suhteliselt madal.Path‑analüüs võimaldab testida mõlemaid hüpoteese korraga ja hinnata otse-, kaudseid ning vahendatud (mediatsiooni) efekte.
Path‑analüüs on tavaliselt osa strukturaalsest võrrandimudelist (SEM). Antud juhul saab kõiki muutujaid käsitleda mõõdetavatena – need on otse saadaval andmekogus. Kui soovite lisada abstraktseid konstruktsioone (nt “Liiklusrisk”), võite luua latentseid muutujad, kuid see ei ole kohustuslik.
VMT
|
v
AA ---> FA ---> TF ---> RFR
\ | ^ ^
\_______|_______/ |
(kaudne) |
------------
Otsefektid:
AA → FA – kui õnnetusi on rohkem, siis tõenäoliselt ka surmaga lõppevaid.FA → TF – iga fatal accident annab vähemalt ühe surma; seega tugev lineaarne seos.TF → RFR – suremuste arv normaliseeritud sõidukikilomeetri kohta.Kaudsed (mediatsioon) teed:
AA → FA → TF → RFR.AA → RFR, et testida, kas sagedusel on otsene mõju ka ilma raskusastme vahenduseta.Kontrollmuutujad (X1…Xk) võivad mõjutada kõiki nelja põhimuutuja:
Xi → AA, Xi → FA jne.Klassikaline path‑analüüs keskendub lineaarsetele regressioonidele, mis on lihtsamini tõlgendatavad:
FA = β1 * AA + Σγi * Xi + ε1
TF = β2 * FA + Σδi * Xi + ε2
RFR = β3 * TF + Σθi * Xi + ε3
AA (või VMT) on eksogeensus, seega ei ole eesmine regressioon.
| Samm | Kirjeldus |
|---|---|
| a. Puuduvate väärtuste käsitlemine | Kasutage täitmisstrateegiat (imputation) või eemaldage puuduvad ridade komplektid, kui nende osakaal on väike (<5 %). |
| b. Normaliseerimine / log‑transformatsioon | Surnu‑intensiivsuse ja õnnetuste arvud võivad olla tugevalt asümmeetrilised; logaritmiline transformatsioon (log(AA+1)) vähendab heteroskedastilisust. |
| c. Aja‑seeria struktureerimine | Kui andmed on igapäevased, võite luua kuukesksed või aastasõltuvad koondsummad, et vältida liialt suurt korrelatsiooni automaatsete trendide tõttu. |
| d. Kontrollmuutujate lisamine | Näiteks WeatherIndex (halvem ilm = 1, hea ilm = 0), RoadType (maantee klass). |
lavaan, semPlot.semopy või statsmodels koos käsitsi määratud regressioonidega.lavaan)library(lavaon)
# Andmete laadimine (eeldame, et df on andmeraamistik)
df$logAA <- log(df$AllAccidents + 1)
df$logFA <- log(df$FatalAccidents + 1)
df$logTF <- log(df$TotalFatalities + 1)
model_spec <- '
# otse‑efektid
FatalAccidents ~ b1*AllAccidents + c1*WeatherIndex + c2*RoadType
TotalFatalities ~ b2*FatalAccidents + d1*WeatherIndex + d2*RoadType
RFR ~ b3*TotalFatalities + e1*VehicleMilesTraveled
# kaudne tee (AllAccidents → FatalAccidents → TotalFatalities → RFR)
'
fit <- sem(model_spec, data = df, estimator = "MLR") # robust ML
summary(fit, standardized=TRUE, fit.measures=TRUE)
Tulemuste tõlgendamine:
β) näitab mõju suurust (0.00–1.00).b3 on näiteks 0.70, siis iga standardse ühenemise kasv total fatalities-is tõstab RFR-i 0.7‑standardset ühikut.β_total = b1 * b2 * b3. See näitab, kui palju kogu AA muutus mõjutab RFR-i läbi mediatsiooni.| Mõõde | Tähendus |
|---|---|
| χ² test (Chi‑square) | Ideaalis ei olgu p‑väärtus <0.05, kuid suurte andmekogude puhul on see tihti signifikaalne; seega kasutage alternatiivseid näitajaid. |
| RMSEA (Root Mean Square Error of Approximation) | < 0.06 → hea sobivus. |
| CFI / TLI (Comparative Fit Index, Tucker‑Lewis Index) | > 0.95 soovitatav. |
| SRMR (Standardized Root Mean Square Residual) | < 0.08 korral peetakse aktsepteeritavaks. |
Kui mudel ei vasta nendele kriteeriumidele, võite:
AA + FA).VIF < 5).| Efekt | Tulemuse tüüp | Interpretatsioon transpordiplaneerijale |
|---|---|---|
b1 = 0.45 (p<0.001) |
Otseeffect: AllAccidents → FatalAccidents |
Iga 10 % tõus õnnetuste arvus suurendab fatal accidents’i keskmiselt 4,5 %. |
b2 = 0.80 (p<0.001) |
Otseeffect: FatalAccidents → TotalFatalities |
Surmaohtlikkuse suurenemine on peamine driver surmade arvule. |
b3 = 0.62 (p<0.01) |
Otseeffect: TotalFatalities → RFR |
Iga 1 % kasv suremuste koguarvusest tõstab surmaintensiivsust 0,62 % (kui VMT on konstantne). |
Kaudne efekt (AA → FA → TF → RFR) = 0.45 × 0.80 × 0.62 ≈ 0.22 |
Kaudne mõju | Ühe protsendi tõus õnnetuste arvus (kõik) suurendab surmaintensiivsust ligikaudu 0,22 % läbi mediatsiooni. |
AA → RFR otseteed (β = 0.05) |
Väike otsene efekt | Enamik mõju kulgeb vahendatud teekonna kaudu; seega keskendumine ainult õnnetuste vähendamisele (ilma raskusastet käsitlemata) annab väikese tulemuse. |
Järeldus: Raskusaste (FA, TF) on tugevam "vahendus" kui puhas sündmuste sagedus, mis tähendab, et poliitikate fookus peaks olema õnnetuste graviteedi vähendamine, mitte ainult nende arvu piiramine.
| Analüüsi tüüp | Kirjeldus & kasutus |
|---|---|
| Mitte‑lineaarne mediatsioon (logistlik või Poisson) | Kui AA on loendandmed, kasutage Poisson või negative binomial regressioone. Path‑analüüs saab kombineerida GLM‑idega (lavaan toetab). |
| Ajaregressiooni / dynamiiline SEM | Lisage viiteaegade (lag) efekti: FA_t = β1 * AA_{t-1} + …. See võimaldab näha, kuidas muutused ühe aasta jooksul mõjutavad järgmist aastat. |
| Mitmetasandiline mudel (hierarchical) | Kui andmed on jaotatud piirkondade või riikide kaupa, kasutage mitmetasandilist SEM (lavaan + cluster). |
| Sõltuvus‑analüüs | Võrrelge erinevate alampiirkondade koefitsiente – kas suured linnad näitavad sama mediatsiooni kui maapiirkonnad? |
FA ja TF on kõrged (kaasa arvatud VMT‑normeeritud).ΔFA / ΔAA – “raskusastme muutumise suhe”.ΔRFR per 10 % VMT vähenemine.AA → FA → TF → RFR).R + lavaan).semPlot), tabelid.Path‑analüüs annab struktureeritud raamistikku, millega saab:
Kasutades eespool kirjeldatud sammude komplekti, saate luua usaldusväärse statistilise mudeli, mis toetab tõendus‑põhist otsustamist ja aitab saavutada turvalisemat ning jätkusuutlikumat transpordisüsteemi.
| Rõhutatud eesmärk | Kuidas VMT aitab? |
|---|---|
| Eksponeerimise normaliseerimine – surmad ja õnnetused arvutatakse “surnud / miljon sõidukikilomeetrit”. | Võimaldab võrrelda piirkondi, kus liiklusmaht on väga erinev. |
| Kasutus‑intensiivsuse hindamine (nt 1 km teel mitu õnnetust). | Annab mõõdikud nagu “õnnetused / km” või “surmad / km”. |
Kui VMT ei ole kättesaadav, tuleb leida alternatiivseid eksponeerimise mõõtmeid või modelleerida VMT latentse muutujana, et säilitada tee‑analüüsi sisuline loogika. Järgnevas juhendis on toodud samm‑sammult, kuidas seda teha.
| Proxy | Kirjeldus / allikas | Eelised ja piirangud |
|---|---|---|
| Liiklusvoog (AADT, Annual Average Daily Traffic) | Tunnistatud liikluse arv (sõidukid päevas) – saadaval enamikus riiklikest teetööde andmebaasidest. | Hea asendus, kuid ei hõlma sõidu kestust ega kiirust; eeldab stabiilsust aastatel. |
| Sõidetud tee‑kilomeetrid (Road‑km) | Teede kogupikkus teatud piirkonnas või klassi järgi (nt maantee X km). | Lihtne mõõde, kuid ei kajasta kasutusintensiivsust – eeldab homogeenset liiklust. |
| Sõidumahu indeks (Traffic Volume Index) | Kombineeritud näitaja: AADT × keskmine sõiduki kiirus või koormuse klassifikaator. | Võib paremini peegeldada tegelikku “kilomeetrite” kasutamist, kuid vajab täiendavaid andmeid (keskmine kiirus). |
| Sõidukite registreerimise arv (Vehicle Stock) | Registritud sõiduautode koguarv piirkonnas. | Sobib kui ligikaudne proxy; ei kajasta tegelikku liiklust, vaid potentsiaalset nõudlust. |
| Inimeste/maja‑suurus per ruutkilomeeter – demograafilised faktorit (populatsioon tihedus). | Inimese arv / km² või elamute arv / km². | Kasutatav, kui liiklusvoog on tugevalt seotud elanike tihedusega; eeldab lineaarset seost. |
| Sõidukite kiiruse‑ja‑koormus‑kaalu (Speed‑&‑Load Index) | Keskmine sõidukiirus × AADT / teeklassi koefitsient. | Pakub keerulisemat, kuid realistlikumat proxy’t “kilomeetrite” kasutamiseks. |
Valiku kriteeriumid:
AA → FA → TF → RFR
↑ ↑
Proxy1 Proxy2 (eksponeerimise proxy)
Proxy1 – mõjutab FA. Näiteks AADT või sõidukite registreeritud arv.Proxy2 – võib olla sama kui Proxy1, kuid on sisestatud ka otseselt RFR‑le, et kontrollida eksponeerimise mõju.| Võrrand | Kirjeldus |
|---|---|
FA = β1·AA + γ1·Proxy + ε1 |
Fatal accidents sõltuvad õnnetuste arvust ja liiklusvoost. |
TF = β2·FA + γ2·Proxy + ε2 |
Surmade koguarv sõltub fatal accidents‑ist ning samasugusest eksponeerimisest (näiteks “surnud / km”). |
RFR = β3·TF + γ3·Proxy + ε3 |
Lõplik surmaintensiivsus (surmad per miljon sõidukikilomeetrit) arvutatakse proxy‑ga normaliseerides. |
Märkus: Kui kasutate sama proxyt kõigis võrrandites, siis see toimib “eksponeerimise kontrollkõlblikuna”. Alternatiivselt võib kasutada erinevaid proksie (nt AADT
Proxy1ja teede pikkusProxy2) – sel juhul hindab mudel mõlemat eraldi.
Kui teil on mitu proxy‑andmestikku, saate modelleerida VMT-i latentse muutujana:
Latent_VMT → AA
Latent_VMT → FA
Latent_VMT → TF
Latent_VMT → RFR (kaudne mõju)
SEM‑raamistik (lavaan või semopy) võimaldab määrata latentse konstruktsiooni ja hinnata selle koefitsiente.
model_latent <- '
# Latent VMT definitsioon (3 indikatsiooni)
VMT =~ AADT + RoadKm + AvgSpeed
# Struktuursed seosed
FatalAccidents ~ b1*AllAccidents + l1*VMT
TotalFatalities ~ b2*FatalAccidents + l2*VMT
RFR ~ b3*TotalFatalities + l3*VMT
'
fit_latent <- sem(model_latent, data=df)
summary(fit_latent, standardized=TRUE)
RFR_proxy = (TotalFatalities / AADT_total) × scaling_factorExposure = Σ(AADT_i × RoadKm_i) – see annab ligikaudse “kilomeetri‑päeva” väärtuse.RFR_proxy2 = TotalFatalities / ExposureRFR_stock = TotalFatalities / VehicleStockRFR_pop = TotalFatalities / PopulationOluline: Iga proxy‑valik peab olema konstantne kogu analüüsiperioodi jooksul või korrektselt normaliseeritud (nt inflatsiooni, rahvastiku kasv).
| Samm | Tegevus | Kasutatav tööriist / koodinäide |
|---|---|---|
| 1. Andmete kogumine | Koguge AllAccidents (AA), FatalAccidents (FA), TotalFatalities (TF) ja valitud proxy‑andmed (AADT, RoadKm vms). |
Excel/CSV → R/Python import |
| 2. Proxy normaliseerimine | Arvutage exposure – nt Exposure = Σ(AADT_i × RoadKm_i). Kui ainult AADT, siis Exposure = Σ AADT_i. |
R df$Exposure <- with(df, sum(AADT * RoadKm)) |
| 3. Muutujate transformatsioon (valikuline) | Log‑transformeerige loendandmeid, et vähendada heteroskedastilisust: logAA = log(AA+1) jne. |
R df$logFA <- log(df$FatalAccidents + 1) |
| 4. Mudeli definitsioon | Kirjutage SEM‑mudel (näide allpool). | R lavaan või Python semopy. |
| 5. Hinnang ja sobivuse kontroll | Kontrollige χ², RMSEA, CFI, SRMR – veenduge, et mudel on statistiliselt vastuvõetav. | R summary(fit, fit.measures=TRUE) |
| 6. Kaudsete efektide arvutamine | Korrutage koefitsiente (nt β_total = β1·β2·β3). |
Manuaalselt või funktsiooniga (parameterEstimates()). |
| 7. Tulemuste visualiseerimine | Joonistage tee‑diagramm, kaardistage riskipiirkonnad. | R semPlot, QGIS/ArcGIS (kaardi loomine). |
| 8. Poliitikasoovitused | Koostage raport: “Sõidukite registreerimise arvu vähendamine X % → surma‑intensiivsuse vähenemine Y %”. | Word/LaTeX + graafikud. |
library(lavaan)
# 1) Andmed
df$logAA <- log(df$AllAccidents + 1)
df$logFA <- log(df$FatalAccidents + 1)
df$logTF <- log(df$TotalFatalities + 1)
# 2) Proxy: AADT (kui teekilomeetrite puuduvad, kasutame ainult voogu)
exposure <- sum(df$AADT) # päevane sõidukite arv kogu perioodis
df$RFR_proxy <- df$TotalFatalities / exposure * 1e6 # surmad per miljon AADT
# 3) SEM mudel (kasutades log‑muutujaid)
model <- '
# otseefektid
logFA ~ b1*logAA + g1*AADT
logTF ~ b2*logFA + g2*AADT
RFR_proxy ~ b3*logTF + g3*AADT
# kaudne tee (logAA → logFA → logTF → RFR)
'
fit <- sem(model, data = df, estimator="MLR")
summary(fit, standardized=TRUE, fit.measures=TRUE)
# 4) Kaudse efekti arvutus
est <- parameterEstimates(fit)
beta_total <- est$est[est$label=="b1"] *
est$est[est$label=="b2"] *
est$est[est$label=="b3"]
cat("Kaudne efekt (logAA → RFR):", beta_total, "\n")
| Kontroll | Kuidas teha? |
|---|---|
| Erinevate proxyde võrdlus | Käivita mudel kolm korda – kordamööda AADT, RoadKm, VehicleStock. Võrreldes koefitsiente ja sobivusindikaatoreid (RMSEA jne). |
| Bootstrapping | Hinnake kaudsete efektide usaldusvahemikke 5 000‑kordse bootstrapiga (lavaan funktsioon bootstrap). |
| Ajaliselt lagged mudel | Lisage viiteaegade (lag) muutujad: logFA_t ~ logAA_{t-1}. Kontrollib, kas mõju ilmneb järgmisel aastal. |
| Subgrupi analüüs | Eraldi hinnake linnalised vs maapiirkondlikud andmed – võib selgitada, milline proxy töötab paremini erinevates kontekstides. |
Kui need sammud järjekindlalt järgida, siis võite saavutada sama informatiivse ja poliitikakujundajatele kasuliku tulemuse nagu VMT‑põhise mudeliga – näidates selgelt, kuidas õnnetuste sagedus vs raskusaste mõjutab surmaintensiivsust ning millised sekkumised toovad kõige suurema turvalisuse tõusu.
| Tegevus | Juhend / Tööriist |
|---|---|
| 📊 Andmete kogumine – AA, FA, TF + AADT või muu proxy | Statistikaamet, liikluskorralduse andmebaas |
| 🧹 Puhastamine & transformatsioon (log‑transform) | R dplyr, Python pandas |
| 📐 Proxy arvutamine – Exposure = Σ(AADT·RoadKm) vms | Excel/SQL / R mutate() |
| 🔍 SEM mudeli definitsioon | R lavaan, Python semopy |
| 🧪 Sobivuse hindamine (RMSEA, CFI…) | summary(fit, fit.measures=TRUE) |
| 📈 Tulemuste visualiseerimine – tee‑diagrammid + riskikaardid | R semPlot, QGIS/ArcGIS |
| ✍️ Poliitikasoovitused (prioriteedid) | Word / LaTeX raport koos KPI‑dega |
Järgides seda raamist, saate luua tugeva tõendusmaterjali liiklusohutuse parandamiseks ka olukorras, kus traditsiooniline VMT‑andmestik on kättesaamatu.
(ilma VMT‑andmeteta)
Alljärgnevas selgitan, kuidas liiklussagedust – tavaliselt AADT (Annual Average Daily Traffic) või mõnda muud voogumõõdet – saab süstemaatiliselt integreerida teooria‑põhisesse tee‑analüüsi raamistikku.
| Liiklussageduse roll | Mõju mudelile |
|---|---|
| Eksponeerimise faktor – näitab, kui palju sõidukeid teekond igapäevaselt kasutab. | Võimaldab “surma‑intensiivsust” (RFR) skaleerida vastavalt tegelikule koormusele. |
| Riskifaktor – suurem voog suurendab kokkupõrke tõenäosust, kuid võib ka vähendada sõiduki intensiivsuse (kui sama arv sõidu‑km jaotub rohkemate auto vahel). | Selgitab, miks mõnikord kõrgem AA ei tähenda proportsionaalselt suuremat RFR. |
| Seadusandlik alus – AADT on tihti sisuliselt olemas riiklikes teede haldamise andmebaasides ja seetõttu kättesaadav. | Pakub praktilist mõõdet, millega saab mudelit rakendada ilma VMT‑spetsiifiliste GPS‑andmeteta. |
AA → FA → TF → RFR
│ ↑ |
└─► Exposure (AADT) ────►
Exposure on liiklussagedus‑kaalutud koefitsient, mis korrigeerib iga järgmise astme sõltuvust.| Võrrand | Kirjeldus |
|---|---|
FA = β1·AA + γ1·AADT + ε1 |
Fatal accidents sõltuvad koguarvust (AA) ja liiklussagedusest (AADT). |
TF = β2·FA + γ2·AADT + ε2 |
Surmade arv sõltub fatal accidents‑ist ning AADT‑st. |
RFR = β3·TF + γ3·AADT + ε3 |
Lõplik surmaintensiivsus (surmad / miljon sõidukikilomeetrit) on mõjutatud nii suremuste koguarvust kui ka liiklussagedusest. |
Märkus: Kui
RFRarvutatakse normeeritud AADT‑ga, siis viimane koefitsient (γ3) võib olla 0 (kuna normaliseerimine juba võtab AADT arvesse). Sel juhul onAADTsisestatud ainultFAjaTFtasanditel.
Liiklussagedus on tihti väga mitmekesine (10 000–200 000 sõidukit päevas). Logaritmiline skaleerimine aitab:
logFA = β1·logAA + γ1·logAADT + ε1
logTF = β2·logFA + γ2·logAADT + ε2
logRFR= β3·logTF + γ3·logAADT + ε3 (või jätke see välja, kui RFR on juba log‑skaalas)
Kui teil on lisaks AADT‑le ka teisi koormuse näitajaid – näiteks:
| Proxy | Kirjeldus |
|---|---|
RoadKm |
Teede kogupikkus piirkonnas. |
AvgSpeed |
Keskmine sõidu kiirus (km/h). |
VehicleStock |
Registreeritud sõidukite arv. |
Siis saate luua latentse muutujaga, mis esindab tegelikku koormust:
TrafficVol =~ AADT + RoadKm + AvgSpeed
FA ~ β1*AA + λ1*TrafficVol
TF ~ β2*FA + λ2*TrafficVol
RFR ~ β3*TF + λ3*TrafficVol
=~ tähistab, et TrafficVol on latentne faktor, mida mõõdetakse antud indikaatoritega.λ1, λ2, …) optimeerida, nii et mudel “õpib”, milline proxy annab parima selgituse.library(lavaan)
model_latent <- '
# Latentne faktor: Traffic Volume
TrafficVol =~ AADT + RoadKm + AvgSpeed
# Struktuursed seosed
FatalAccidents ~ b1*AllAccidents + l1*TrafficVol
TotalFatalities ~ b2*FatalAccidents + l2*TrafficVol
RFR ~ b3*TotalFatalities # (RFR on juba normaliseeritud)
'
fit_latent <- sem(model_latent, data = df,
estimator = "MLR", std.lv=TRUE) # standardiseeritud latentne faktor
summary(fit_latent, standardized = TRUE, fit.measures = TRUE)
Tulemuste tõlgendamine
l1, l2 näitavad, kui palju liiklussagedus (ja sellega seotud koormuse komponent) mõjutab otse fatal accidents ja total fatalities.l1 või l2 on statistiliselt olulised (p < 0.05) – siis liiklussagedusel on oluline roll riskifaktorina.| Allikas | Mida koguda? |
|---|---|
| Riiklik teede haldamise andmebaas (nt Eesti Maanteeamet) | AADT, teekilomeetrid (RoadKm), keskmine kiirus. |
| Statistikaamet / Demograafiad | Sõidukite registreerimise arv, elanike tihedus – võimaldavad täiendavaid proksie. |
| Ilmaandmed (meteostaatika) | Ilmastikutingimused (WeatherIndex) – kontrollmuutujaks, sest halb ilm suurendab õnnetuste tõenäosust. |
Andmetöötluse sammud
Exposure = Σ (AADT_i × RoadKm_i) – see annab “sõidukite‑päeva‑kilomeetriline” koormus, mis on VMT‑asendus.Logi transformeerige, kui jaotused on tugevalt kallutatud (log(Exposure+1)).
df$Exposure <- with(df, AADT * RoadKm) # rida‑tasemel
total_exposure <- sum(df$Exposure) # kogu piirkonna koormus
df$logExposure <- log(total_exposure + 1)
| Efekt | Tähendus (näide) |
|---|---|
γ1 positiivne ja oluline (p < .01) |
Liiklussagedus suurendab fatal accidentsi – seega tuleks keskenduda koormuse vähendamisele (nt alternatiivne ühistransport, takistuste loomine). |
β2 tugev (≈0.8) |
Surmade koguarv sõltub peamiselt fatal accidents‑idest; see tähendab, et raskusastme alandamine (kiirusepiirangud, turvalised ristmikud) annab suurema mõju kui ainult õnnetuste arvu vähendamine. |
γ2 negatiivne |
Suurem AADT vähendab suremusi per sõidukikilomeeter – võimalik seletus on “turvalisemad teed/kiirusepiiramised” suurema koormusega maanteedel. See annab signaali, et kõrge liiklusvoog ei ole alati negatiivne; oluline on koormuse kvaliteet. |
λ1 (latentse faktori kaalu) |
Kui AADT omab suurimat kaalu (≈0.7) ja teised proksid väiksemad, siis AADT on peamine liiklussageduse mõõde. Planeerijad peaksid keskenduma sõiduvoo juhtimisele (nt dünaamiline kiirusepiirang). |
Liiklustegevuse hajutamine
AADT on kõrge ja γ1 positiivne.Raskusastme vähendamise meetmed
β2).Sõiduvoo optimeerimine – “Smart Traffic Management”
γ1‑efekti (väheneb fatal accidents).Järelevalve & KPI-d
| Samm | Tegevus |
|---|---|
| ✅ Koguge AADT + täiendavad proksid (RoadKm, AvgSpeed, VehicleStock). | |
✅ Arvutage Exposure = Σ(AADT × RoadKm) või kasutage log‑transformeeritud logAADT. |
|
| ✅ Määrake mudeli struktuur: - Otse regressioon ( FA ~ AA + AADT) - Või latentne faktor ( TrafficVol). |
|
| ✅ Hinnake mudelit (RMSEA, CFI, SRMR). | |
✅ Tõlgendage koefitsiente – keskendu γ‑ ja β‑parameetritele. |
|
✅ Koostage poliitikasoovitused vastavalt sellele, milline komponent (AADT, FA, TF) on kõige tugevam riskifaktor. |
|
| ✅ Jälgige KPI-sid ning uuendage mudelit iga 2‑3 aasta tagant uute andmetega. |
Kasutades ülaltoodud raamistikku, saate luua tõendus‑põhise analüüsi, mis on täielikult rakendatav ka siis, kui VMT‑spetsiifilised andmed puuduvad, ning annab selge suuna jätkusuutliku ja turvalise transpordisüsteemi kujundamiseks.
Kuidas rakendada mitmeindikaatori liiklusohutuse (RTS) mudelit teistes riikides või piirkondades
Alljärgnevalt esitan praktilise töövoo, mis põhineb artiklis kirjeldatud lähenemisel ning on kohandatav erinevate rahvuslike ja regionaalsete kontekstide jaoks.
| Samm | Mida teha | Miks oluline |
|---|---|---|
| a) Kogukonna‑ ja poliitiline analüüs | Kaardista riigi/piirkonna olemasolevad liiklusohutuse strateegiad, SDG‑si kohaldamine (nt eesmärk 3.6 – surmade pooleks vähendamine aastaks 2030) ning rahvuslikud sihid (näiteks “Towards Zero” – 50 % vähenemine võrreldes baasajaga 2010) [1]. | Tagab, et mudel toetub olemasolevatele poliitilistele raamistikutele ja aitab siduda uue mõõdikuga juba kehtivad eesmärgid. |
| b) Andmekaitse‑ ja õigusraamistik | Kontrolli, millised andmed on kättesaadavad (politsei‑/sõidukiregistrite aruanded, telemaatikasüsteemide logid). Veendu, et nende kasutamine vastab kohalikele privaatsusnormidele. | Andmete kvaliteet ja õiguslik korrektsus mõjutavad mudeli usaldusväärsust. |
| c) Geograafilised/ demograafilised eripärad | Kogu teave sõidukite tüübi, liikluskoormuse tiheduse, linna‑ vs maapiirkondade kohta. | Need faktorid võivad muuta indikaatorite suhtelisi kaalusid (nt maapiirkondades on õnnetuste sagedus suurem, kuid tõsidus madalam). |
| Indikaator | Andmeallikad | Märkus |
|---|---|---|
| Road fatality rate (surmajuhtude arv per miljon sõidetud km) | Sõidukite läbitud kilomeetrid (tollisüsteem, GPS‑andmed või statistiline modelleerimine) + registreeritud surmad. | Kui täpset kaubandusliku liikluse andmestikku pole, võib kasutada keskmist sõidukilometrit per sõiduk ja korrutada registrite arvuga. |
| Number of fatalities (kogu surmade arv) | Politsei/terve‑süsteemi registreerimised; haiglate suremusregistrid. | Kontrolli andmete täielikkust – alaraportimine on tavaline probleem, seega kasuta mitme allika võrdlust. |
| Number of fatal accidents (fataalsed õnnetused) | Politsei õnnetusaruanne; kindlustusseltside kahjukirjed. | Täpsusta definitsioon – “fatal accident” tähendab sündmust, kus vähemalt üks inimene sureb kohapeal või 30 päeva jooksul pärast kokkupõrget. |
| Total number of road accidents (kõik õnnetused) | Samad allikad kui eelmises veerus, kuid laiem filtr – ka väikeste kahjudega sündmused. | Võib lisada kategooriad: “light injury”, “serious injury”. |
Andmete puhastamine: lineaarne interpolatsioon puuduvate väärtuste täitmiseks; log‑transformatsioon (log(x + 1)) suurte arvude stabiliseerimiseks.
Road fatality rate →\
Number of fatalities → RT S
Fatal accidents →/
Total accidents →
Lisaks võid lisada kaudseid teid, näiteks:
road fatality rate → number of fatalities → RTS, et näidata, kuidas suremus mõjutab üldist turvalisust kaudselt.
β) väärtused määravad indikaatorite kaalu RTS‑indeksis.R (pakett lavaan), Stata (sem) või IBM SPSS AMOS – kõik toetavad path‑analysisi ja võimaldavad tulemuste visualiseerimist.
| Valideerimismeetod | Kirjeldus |
|---|---|
| Kvalitatiivne kontroll | Intervjueeri kohalikke liiklusohutuse eksperte, et kinnitada indikaatorite tähendus. |
| Rist‑validatsioon (kross‑validation) | Jaga andmed treening‑/testkomplektiks (nt 70 % /30 %). Hinnake prognoosiv täpsus (RMSE või MAPE). |
| Ajaline stabiilsus | Kontrolli, kas koefitsientide väärtused püsivad sarnased erinevate ajaperioodide vahel (2010‑2020 vs 2021‑2025). |
[ RTS_{i,t}= w1 \cdot FR{i,t}+w2 \cdot NF{i,t}+w3 \cdot FA{i,t}+w4 \cdot TA{i,t} ]
FR = road fatality rate,NF = number of fatalities,FA = fatal accidents,TA = total accidents.Kaalud (w_j) võetakse path‑analysis‑st saadud β‑koefitsientidest.
[
RTS^{*}=100 \times \frac{RTS-\min(RTS)}{\max(RTS)-\min(RTS)}
]
Võimaldab riikide ja piirkondade võrdlust ühtsel skaalal.
| Dominantsen indikaator | Soovitatav sekkumine |
|---|---|
| Road fatality rate (kõrge) | Kiirusepiiramiste automaatne jälgimine, nutikad kiiruskaamerad. |
| Number of fatalities (kasvav) | Kohustuslik turvalisuskoolitus; rangemud sõidukite ohutusstandardid (AD‑AS). |
| Fatal accidents (suured arvud) | Ristmiste ümberkujundamine “protected intersections”, jalakäijate eraldatud rajad. |
| Total accidents (palju, kuid madal tõsidus) | Kampaaniad turvalise sõidu kultuuri edendamiseks; teavitustööd alkoholi‑ ja mobiiltelefonikasutuse kohta. |
| Väljakutse | Lahendus |
|---|---|
| Andmekvaliteedi varieeruvus (puuduvad või ebatäpsed andmed) | Loo Data Quality Index ja kasutage mitme allika triangulatsiooni. |
| Indikaatorite definitsioonide erinevused riikides | Koosta ühtne meta‑andmestiku koodiraamistik, lähtudes IRTAD või WHO standarditest (nt “fatal accident” = 1+ surma). |
| Poliitiline toetus puudub | Esita mudeli pilotuuringud ja näita selgelt seost SDG 3.6‑ga ning potentsiaalset kuluefektiivsust [1]. |
Mitmeindikaatori RTS‑mudel on kohandatav järgmistest sammudest: (i) kohaliku konteksti kaardistamine, (ii) andmete kogumine ja puhastamine, (iii) path‑analysisi põhine struktuurimudel, (iv) mudeli kalibreerimine ning (v) indeksilepõhiste poliitikasoovituste koostamine. Sellise lähenemisega saab iga riik või piirkond täpselt mõõta oma liiklusohutuse olukorda, jälgida edusamme SDG‑de suunas ja suunata ressursse tõhusamalt õnnetuste tõsiduse vähendamiseks – mis on uuringu põhijäreldus (et “severity” mõjutab RTS-d tugevamini kui „frequency”) [1] 7516011e.md.