Tervis ja tööjõud – mõned joonised

(Lühem versioon ilmus Arenguseire Keskuse uudiskirjas.)

McKinsey rahvusvahelise instituut avaldas suvel mahuka analüüsi “Prioritizing health: A prescription for prosperity” parema tervise positiivsest mõjust majandusele. Eesti kohta tehtud näitearvutuste põhjal pakkus instituut, et parema tervisega inimeste korral võiks Eesti sisemajanduse koguprodukt olla 20 aasta pärast 11% suurem (link kaardile, kust näeb Eesti andmeid).

Seda lugedes tekkis soov tõsta kokku mõned Eesti andmeid, mis iseloomustavad tervise mõju majandusele.

Parema tervise mõju majandusele avaldub läbi erinevate kanalite:

1. Parem tervis suurendab tööjõu hulka: inimesed ei sure tööealistena, nad töötavad kõrgema vanuseni ja nad töötavad suurema töökoormusega.
2. Parema tervisega inimesed on töökohal tootlikumad.
3. Parema tervisega inimesed saavad omandada parema hariduse, mis avaldab omakorda mõju tootlikkusele.
4. Parema tervisega inimesed säästavad enam, mis loob vahendid investeerimiseks.

Terved inimesed vajavad ka vähem tervishoiu- ja hoolekandekulusid, olgu tehtud kas leibkonna või riigieelarve vahenditest, mis loob võimaluse kasutada seda raha muudeks tegevusteks.

Viimane laiem Eesti kohta tehtud uuring jääb 2006. aastasse, kui valmis mahukas Praxise, sotsiaalministeeriumi ja Maailma Terviseorganisatsiooni analüüs tervise mõjust Eesti majandusele (link). Uuring näitas, et Eestis on haigestumine oluline tegur töötajate enneaegsel hõivest lahkumisel ning kehval tervisel on ka märkimisväärselt negatiivne mõju töötundidele ja palgale. Teiste riikide kogemuse põhjal pakkus see uuring, et täiskasvanute suremuse vähendamisel 1,5% võrra aastas võiks 25 aasta pärast tulemuseks olla 14% kõrgem SKP inimese kohta.

Hilisemaid üksikuid uuringuid Eesti kohta on veel mitmeid, nt halva tervise mõju kohta tööturult eemalejäämisele (Sõrmus (2014), Karilaid-Vidder (2015)) või haiguspäevadele (Aaviksoo 2017).

Alljärgnevalt ongi toodud mõned joonised, mis iseloomustavad, kui suur osa Eesti tööjõu potetnsiaalist jääb kasutamata kehva tervise tõttu.

1. Enneaegne suremus

Esimene kanal, kuidas tervis mõjutab töötajate arvu on inimeste enneaegne suremus. Näiteks on Eesti meeste suremustõenäosused vanusevahemikus 30-70 enam kui kaks korda suuremad kui Rootsi meestel. Erinevused naistel suremustõenäosustes ei ole nii suured, jäädes enamasti alla kahe korra.

See tähendab, et kui säiliks 2018. aasta suremuskordajad, siis 2018. aastal sündinud 7345 poisslapsest elaks 65 aasta vanuseks vaid kolmveerand (5579). Kui aga Eesti noored mehed elaksid tervislikumalt ja meil oleks Rootsi meeste suremuskordajad, siis elaks pensionieani lausa tuhat poissi rohkem (6606). Kokku elaks üks põlvkond mehi oma parimas tööeas (vanuses 20-65) enam ligi 11 tuhat inimaastat rohkem (3,6%).

2. Haiguse tõttu tööturult eemal olek

Statistikaameti andmetel on Eestis aastas keskmiselt 40 tuhat tööealist inimest, kes ei tööta kehva tervise tõttu. Neile lisandub veel 20 tuhat juba pensionieas inimest, kes parema tervise puhul samuti töötaks. Seega ideaalse tervise puhul võiks lisanduda ligi 60 tuhat praegu tööturult väljas olevat inimest tööturule. See moodustab ligi 9% praegu tööturul hõivatud inimestest ning oleks oluline panus majanduskasvu.

Kui võtta võrdluseks taas näiteks Rootsi või Soome, siis näeme, et vanusrühmas 50-64 on meil 3-4% enam inimesi, kes on haiguse või vigastuse tõttu mitteaktiivsed. Kui saavutaksime näiteks selles vanusrühmas kasvõi Soome või Rootsi taseme, siis oleks tööturul rohkem umbes 1000 inimest.

3. Väiksem töölt puudumine ja töötunnid

Haigekassa statistika järgi maksti 2018. aastal haigushüvitist 3,8 miljoni päeva eest. Arvestades juurde hüvitamata päevad, siis puudusid inimesed töölt haiguse tõttu suurusjärgus 2-3% tööpäevadest.

Praxise 2006. aasta uuring leidis, et halva terviseseisundiga inimestel olid nädala töötunnid meestel keskmiselt 12 tunni võrra väiksemad ja naistel 8 tunni võrra väiksemad võrreldes hea terviseseisundiga inimestega.

4. Produktiivsus

Haigestumine mõjutab mitte ainult püsivalt tööjõust lahkumist või ka suremust, vaid ka produktiivsust töökohal. Ühest küljest võib produktiivsust otseselt suurendada paranenud füüsilise ja vaimse aktiivsuse kaudu. Teisalt võivad füüsiliselt ja vaimselt aktiivsemad isikud paremini ja efektiivsemalt kasutada tehnikat, masinaid ja seadmeid. Produktiivsust mõõdetakse tavaliselt inimese palgaga. Ülalmainitud Eesti 2006. aasta uuring, kasutades küll juba 2003. aasta andmeid leidis, et halb tervis (võrreldes hea tervisega) vähendas inimeste kuusissetulekuid oluliselt – meeste puhul peaaegu 30% ja naiste puhul umbes 20% keskmisest palgast.

Ri fail joonistega

Eesti Panga finantsstatistika allasikutamine R-iga

Eesti Pank avaldab enda kodulehel finantsstatistikat, kuid kahjuks mitte väga mugavas formaadis, et seda automaatselt Ri laadida ja siis analüüsida.

Siiski, kui sama tabelit kasutada sageli, saab ühekordse käsitööga ülejäänud protsessi automatiseerida.

Näiteks soovime vaadata majapidamiste hoiuste muutust, mis asuvad tabelis 1.1 Krediidiasutuste ja liisinguettevõtete analüütilised näitajad (miljon eurot). Tabeli esmakordsel käsitsi allalaadimisel (olgu Exceli või CSV failina) saame alla laetud faili lingist teada raporti id (Reports?id), alguskuupäeva (date1), lõpukuupäeva (date2), keele (lng=et), allalaetavad väärtused, nt väärtus, muutus või aastane muutus. (fa=,VAARTUS,MUUT,MUUT_AA). Neid saab siis juba ise hiljem automaatses päringus muuta Ris.

R-is saame kasutada käsku GET paketist httr, et andmefail alla laadida koos kõigi atribuutidega. Andmete kättesaamiseks peame eraldama sisu (käsuga content) ning seejärel saame teksti muuta andmetabeliks käsuga read.delim.

library(httr)
temp <- GET("https://statistika.eestipank.ee/Reports?id=868&format=csv&date1=2004-01-1&date2=2020-07-31&lng=et&fa=,VAARTUS")
sisu <- content(temp, "text")
dfalg <- read.delim(textConnection(sisu), 
                    header=TRUE, sep=";", 
                    strip.white=TRUE, 
                    blank.lines.skip = TRUE, skip = 4)

Et andmetabel sisaldab ka päist, siis sisse lugemata jäetud rida arv võib olla eri tabelitel erinev. Antud tabelis on esimesed 4 rida pealkirjade ja muu info päralt (skip = 4).

Sellega on andmefail olemas. Et tabel on laial kujul, siis on perioodid veergudes ning tunnuste nimed esimeses veerus. Iga konkreetne tabel võib olla kahjuks natuke erineva struktuuriga.

Meie näitetabelist saame võtta välja meid huvitavad tunnused (ridades) ja pöörata tabel ümber pikale kujule, teha kuupäeva ning teha muutuja väärtused numbriliseks. Näiteks võtan 4. ja 12. reast majapidamiste laenud ja hoiused ning teen joonise.

library(dplyr)
library(tidyr)
library(ggplot2)

df <- dfalg[c(4,12),] #hoian 4 ja 12 rea, laenud ja hoiused 
colnames(df)[1] <- "tunnus" #esimese veeru nimeks tunnus
df$tunnus[1] <- "Laenud" #nimetan ümber
df$tunnus[2] <- "Hoiused" #nimetan ümber
df <-df %  gather(key = "kpv", value = "vaartus", -tunnus) #viin pikale kujule
df <-df % 
  mutate(kpv = as.Date(gsub("X", "", kpv),"%d.%m.%Y"), #kuupäeva tegemine
         vaartus = as.numeric(gsub(",", ".", gsub(" ", "", vaartus))), #väärtusest tühik välja ja koma punktiks
         vaartus12 = (vaartus-dplyr::lag(vaartus, 12))/dplyr::lag(vaartus, 12)*100)  %>%  #aastane juurdekasv protsentides
  filter(!is.na(kpv)) #kui faili lõppu jäi mõni rida, siis viskame välja

#Joonise tegemine
#Pisut puhtamad teemad veebi jaoks
source("http://kodu.ut.ee/~avork/files/Rkoode/Rggplotwebiteemad.R")

ggplot(df %>%  dplyr::filter(kpv>=as.Date("2012-01-01")), aes(x = kpv, y = vaartus12, color = tunnus)) +
  geom_line() +  
  scale_y_continuous(labels = function(x) paste0(x, "%")) + 
  labs(x="", y = "Aastane muutus", color = "", title = "Majapidamiste hoiuste ja laenude aastane juurdekasv")  +
  theme_web_classic()

ggsave(filename = "joonis1naideEPfin.jpg", width = 8, height = 6)

Ja tulemusena näeme, et praeguse kriisi ajal kiirenes majapidamiste hoiuste kasv ning aeglustus majapidamistele antud laenude kasv.

Palju käsitööd, kuid kui uuel kuul tulevad uued andmed, siis piisab päringus muuta kuupäeva ja ongi uuendatud joonis olemas. (Eeldusel, et andmete muu struktuur ja vorm püsib muutumatuna.)

Tänusõnad Märten Veskimäele abi eest.

Ri fail http://kodu.ut.ee/~avork/files/Rkoode/saastmine_blogisse.R

Vanemahüvitise paindlikkus ja isad

Vanemahüvitis on Eesti perepoliitika üks kulukamaid meetmeid. 2019. aastal kulus sellele 262 miljonit eurot. Võrdluseks, et kõikidele teistele peretoetustele (lapsetoetus esimesele, teisele ja kolmandale lapsele, üksikvanema toetus, paljulapseliste perede toetus jm) kulus kokku 303 miljonit eurot.

Vanemahüvitist võib saada ka siis, kui lapsevanem töötab. Alates vanemahüvitise kehtestamisest on neid reegleid muudetud korduvalt ja seda vanemaid soosivamaks. Lapsevanem võib samal ajal teenida ka töötasu. Alljärgnev joonis iseloomustab arvutusreegleid, mis kehtisid erinevatel perioodidel. Pikemalt olen selgitatud seadustes sätestatud valemeid eraldi dokumendis.

Mida hilisem reform, seda rohkem jäi vanemahüvitist alles samal ajal töötamise korral. Näiteks kui 2005. aastal määrati 700-eurone vanemahüvitis ja siis saadi lisaks 2000 eurot brutopalka, jäi lõplikuks vanemahüvitise brutosummaks vaid 250 eurot. 2019. aastal aga jäi sama esialgse vanemahüvitise ja brutopalga korral vanemahüvitise lõplikuks summaks juba enam kui kaks korda rohkem ehk 530 eurot. Ning alates teatud palgast vanemahüvitist enam edasi ei vähendatud alla poole esialgse.

Taoline töötamise ja vanemahüvitise koos saamise võimaldamine on toonud kaasa ka selle, et aina rohkem on hakanud mehed, kes on ka kõrgema eelneva palga ja seega vanemahüvitisega, samaaegselt saama nii vanemahüvitist kui ka töötasu. See võimaldab perel saada suuremat vanemahüvitist, kuid jätkata isal ka töötamist. (Tõsi, ema peab mingil moel saama endale ravikindlustuse, sest lapsehoolduspuhkusel saab olla vaid üks vanem).

On näha, et alates 2018. aasta märtsist, kui kõrgema palga saamist enam ei karistatud nii palju, on järsult kasvanud isade arv, kes saavad vanemahüvitist, ulatudes 2020. aasta juuniks 14%ni vanemahüvitise saajatest. Enamus isadest samal ajal ka töötavad. Ka kevadine majanduskriis ei ole töötamisele mõju avaldanud.

Töötavate isade osakaal vanemahüvitise saajate seas on kasvanud 48%lt 2015. aasta alguses, 71%ni 2020. aasta juunis. Kui varasemalt oli töötavate isade seas pooled need, kes said palka alla alampalga ehk tegid pigem mingeid juhutöid, mis kodus lapsega olemise kõrvalt võisidki võimalikud olla, siis praeguseks saavad peaaaegu kõik töötavad ja vanemahüvitist samaaegselt saavad isad palka üle alammäära.

See kõik võib olla kooskõlas seaduse oodatava mõjuga – emad on päeval tööl, isad tegelevad samal ajal lastega ning kui lapsed lõpuks õhtul magavad, siis kirjutavad isad arvutis Javat ja teenivad kõrget palka juurde. Aga võib olla ka vastupidi, et isad on päeval tööl, lastega ei tegele, kuid pere saab nii isa vanemahüvitist kui töötasu seetõttu, et see on neile kasulikum võrreldes olukorraga, kui ema saaks vanemahüvitist.

Selleks, et neid kahte võimalikku selgitust eristada, oleks vaja, kui sotsiaalkindlustusamet avaldaks ka regulaarselt koondstatistikat, mida hakkavad tegema need ema pärast seda, kui isad hakkavad saama vanemahüvitist – kas emad lähevad (osa-ajaga) tööle või jäävad nad mitteaktiivseks.

Lõpuks ei ole ju midagi selles halba, kui vanemad oskavad hästi matemaatikat ja kasutades seaduse valemeid maksimeerivad leibkonna kogutulu. Iseasi, kas nappide sotsiaalvaldkonna vahendite puhul on siiski tegemist kõige tõhusama meetmega, kuidas suurendada isade rolli laste kasvatamises.

Viited:

• Võrk, A. (2020) Vanemahüvitise suurus töötamise korral läbi aastate, link
• R file jooniste kohta (algfaile ei ole lisatud)
• Sotsiaalkindlustusameti avalik statistika

Töötus, koroona ja Google Trends

Viimast kriisi iseloomustab andmete esitamissageduse kiirus. Uusi koroonajuhte loetakse igapäevaselt. Töötukassa hakkas aprillist avaldama töötuse statistikat igapäevaselt. Statistikaamet avaldab töötamise registri põhjal töösuhete muutuse nädala kaupa ning ettevõtluse kuist kiirstatistikat tuginedes maksudeklaratsioonidele.

Inimeste käitumise analüüsimiseks ja majanduse prognoosimiseks on üks täiendav võimalus vaadata juurde Eesti internetiotsingute märksõnade sagedust Google Trendsi abil.

Järgmisel joonisel on esitatud märksõna “töötukassa” nädalane otsingusagedus (suhtelisel skaalal) ja võrreldud seda töötute arvu muutusega alates 2018. aasta algusest.

Eriolukorra kehtestamise järgselt, 12. märtsil, kasvas märksõna “töötukassa” otsing praktiliselt momentaalselt. Tegelik töötute arvu kasv järgnes mõne nädala pärast. Praeguseks on langenud nii töötukassa otsingusagedus kui ka tegelik töötuse juurdekasv.

Sarnaselt saab vaadata märksõnade “covid” ja “koroona” arengut. Otsing nende kohta kasvas järsult nädal nädal enne eriolukorra kehtestamist, püsis mõnda aega ja seejärel langes. Viimasel paaril nädalal on taas näha tõusutrendi.

Google Trends märksõnade sagedusi on kasutatud paljudes majanduse prognoosimudelites, olgu siis töötuse, inflatsiooni või aktsiahindade prognoosimisel. Google Scholar üle tuhande teadusartikli sellel teemal. Samuti on erinevaid märksõnu kasutada ka Covid viiruse leviku prognoosimisel.

Jääb oodata, et järgmisel õppeaastal on mitmeid tudengeid, kes katsetavad, kas kriisi algust, sügavust ja lõppu aitaksid täpsemalt prognoosida inimeste internetiotsingud.

R-i fail andmete ja jooniste kohta.

Ametirühmade keskmised palgad maksuameti andmete põhjal

Eesti Maksu- ja Tolliamet publitseeris ametikohtade keskmised töötasud, et ettevõtted ja inimesed saaksid oma palku nendega võrrelda. Hetkel on kättesaadavad andmed alates 2019. aasta juulist.

Näiteks olen toonud TOP 20 ametialad (ametikohtade klassifikaatori neljanda taseme järgi). Sinna kuuluvad vaid juhtide ja tippspetsialistide pearühma kuuluvad ametialad.

Nagu oodata on, siis nende järgi tasub noortel õppida meditsiini, infotehnoloogiat või õigusteadust. Või võtta eesmärk saada tippjuhiks või kõrgeks valitsusametnikuks.

Andmete põhjal saab kergesti analüüsida ka ametirühmade keskmiste palkade vahelist hajuvust või konkreetsete ametirühmade vahelist palkade erinevust. Üle pikema perioodi vaadates võiks iseloomustada võimalikke struktuurseid muutusi tööturul.

Alljärgneval joonisel on esitatud 291 ametirühma keskmiste palkade jaotuse üle kuude kõik koos ja seejärel pearühmade kaupa. Periood on veel muidugi liiga lühike, et mingeid trende näha.

Pearühmade kaupa samu jaotusi üle kuude vaadates võib täheldada hinnatud jaotusfunktsioonide muutusi, kuid ametirühmasid enamustes pearühmades väga vähe (joonisel toodud number sulgudes iga pearühma taga). Seetõttu tasub vaadata vaid suuremaid pearühmasid, nt tippspetsialistid.

Viimaks võib püüda hinnata nende andmete põhjal ka käesoleva koroonakriisi mõju. Järgnevalt on esitatud ametirühmade keskmiste palkade muutused aprill-juuni 2020 vs juuni 2019-märts 2020, ehk siis erinevus pärast kriisi võrreldes kriisieelse ajaga.

Ametialade keskmiste palkade langemist on kõigis pearühmades. Kõige enam on palgalangusega ametialasid oskus- ja käsitööliste seas (pooltes ametirühmades palk keskmine langes). Kõige vähem on palgalangusega ametialasid tippspetsialistide seas (vaid iga kuuendas ametirühmas langes palk).

Kindlasti tuleb andmete tõlgendamisel hoolikalt lugeda EMTA metoodikakirjeldust. Andmed sisaldavad vaid töölepingulise suhtega inimeste palkasid. Tööajaerinevusi arvesse ei võeta. Seega ei ole need andmed otseselt võrreldavad Statistikaameti arvutatud keskmiste palkadega. Andmete juures ei ole antud inimeste arvu neis rühmades, mistõttu gruppid ei ole ka võrreldavad suuruse lõikes.

Allikad:

• EMTA. Ettevõtluse statistika ametinimetuste lõikes, https://www.emta.ee/et/ametikohtade-tootasu, vaadatud 30.07.2020
• Ri kood jooniste jaoks. http://kodu.ut.ee/~avork/files/Rkoode/palgadametialadeloikes31072020.R

Töötasu hüvitise olulisusest koroonakriisi ajal

Töötukassa poolt makstud töötasu hüvitise kogukulu moodustas 26. juuli seisuga 256 miljonit eurot, mis on ligi 1% selle aasta prognoositavast sisemajanduse koguproduktist. Võrdluseks võib tuua, et eelmise kriisi ajal, 2009. aastal, moodustasid kõik töötuskindlustuse kulud hüvitistele 1,14% SKPst. Seega on töötasu hüvitise kulud märkimisväärsed, sest sel aastal lisanduvad veel ka kulud tavapärastele töötuskindlustushüvitistele.

Hindamaks töötasu hüvitise olulisust tegevusvaldkondade lõikes võrdleme töötasu hüvitise brutokulu kokku märtsi, aprilli ja mai kuu eest samadel kuudel TSD deklaratsioonidel deklareeritud töölepinguga töötavate brutosummaga.

Kui vaadata töötasu hüvitiste kogukulu tegevusvaldkondade lõikes, siis on oodatult esimeste seas toitlustamise, majutuse, meelelahutusega ja reisimisega seotud tegevusalad, kus töötasu hüvitis moodustas kolmandikust kuni pooleni nende tegevusalade töötasust.

Tööstuse valdkonnast on suurimad töötasu hüvitise kasutajad tekstiilitootmine, trükindus, mööblitootmine ja nahatöötlemine.

Järgmisel joonisel on esitatud seos töötasu hüvitise osakaalu ja aastase käibe languse (märts-mai 2020 võrreldes aasta taguse perioodiga) vahel.

Oodatult on töötasu hüvitise osakaal palgas seda suurem, mida suurem on tegevusala käibe langus.

Kõige suurem käibelangus tabas EMTAKi järgi reisibüroode tegevust (EMTAK 79), majutust (EMTAK 55), muuseumid ja raamatukogud (EMTAK 91) ning loome- ja kunsti ning meelelahutusvaldkonnas (EMTAK 90).

Kui eeldada, et käibe langus peaks olema peamine töötasu hüvitise põhjus, siis võinuks oodata isegi suuremat töötasu hüvitise osakaalu näiteks tegevusalades “91 Raamatukogude, arhiivide, muuseumide ja muude kultuuriasutuste tegevus” ja “59 Kinofilmide, videote ja telesaadete tootmine; helisalvestiste ja muusika kirjastamine”. (EMTAK 19, kus oli ka suur käibe langus, kuid väga vähe töötasu hüvitise saajaid on “Koksi ja puhastatud naftatoodete tootmine”.)

Töötasu hüvitise mõte on hoida majanduskriisi ajal inimesi ettevõtetega seotud, makstes kinni suure osa nende töötasust. Selle eelduseks on see, et kriis on lühiajaline ning ei ole ette näha suuri muutusi tööturu struktuuris peale kriisi. Kriis lõppedes, oleks ettevõtetel lihtne taastada endine tootmismaht.

Kui kriis on aga pikaajaline ja sellega kaasneksid ka muutused tööturul (nt osade ametikohtade kadumine ja teiste tekkimine), siis võib olla otstarbekam maksta pigem kõrgemat töötuskindlustushüvitist ja motiveerida inimesi uusi ameteid vastu võtma.

Teise poolaasta andmete põhjal võiks juba vaadata, kuidas on muutunud hõive nendes tegevusalades, kus oli suur töötasu hüvitise osakaal. Kui vaatamata sellele, et maksti kõrget töötasu hüvitist, neis sektorites ikkagi hiljem hõive kukub oluliselt, siis võib tagantjärgi tarkusena analüüsida, kas ei oleks olnud otstarbekam kohe soodustada inimeste ümberõpet.

Tehnilise märkusena peab lisama, et andmed on kokku pandud mitte ettevõtete vaid tegevusvaldkondade tasandil eri allikatest (Töötukassa ja EMTA). Seetõttu ei pruugi olla tagatud, et samad ettevõtted on klassifitseeritud mõlema asutuse poolt sama EMTAK 2 alla.

Allikad:

• Töötukassa. Töötasu hüvitis seisuga 26.07.2020 https://www.tootukassa.ee/sites/tootukassa.ee/files/asutuste_nimekiri_26.07.2020.xlsx
• EMTA. TSD koondandmed EMTAK2 majandusharude lõikes, väljavõte 23.07.2020
• EMTA. KMD koondandmed EMTAK2 majandusharude lõikes, väljavõte 29.06.2020
• Kogu jooniste aluseks olev tabel http://kodu.ut.ee/~avork/files/Rkoode/tootasuhyvitisjakaive.csv
• Arvutuskood: http://kodu.ut.ee/~avork/files/Rkoode/palgatoetuseolulisus.R

Kui sünnitataks alati keskmiselt

Tuues põhjenduseks tervishoiuteenuste pakkumise kulu-efektiivsemaks ja kvaliteetsemaks muutmise, plaanitakse sulgeda 2018. aastal Põlva ja Valga sünnitusosakonnad. Sünnitajad on teretulnud Tartu ülikooli kliinikumi. Meedias kõlanud TÜ kliinikumi juhtide arvamuste järgi on kliinikumis ruumi piisava varuga, et pakkuda kvaliteetset teenust ja vastata sünnitajate ootustele.

Et meedias on käsitletud keskmisi näitajaid, siis sooviksin alljärgnevas lühianalüüsis rõhutada teist aspekti – sünnituste varieeruvust. Nimelt peab sünnituste arvu igapäevasel prognoosimisel ja nende jaoks ressursside planeerimisel ühe haigla piires arvestama, et tegemist on üsna volatiilse näitajaga. Alljärgneval joonisel on toodud sünnituste päevane aegrida kliinikumis jaanuar 2017-aprill 2018 koos libiseva keskmisega (must joon) ja mulle teadaolevalt avatud sünnitustubade arvuga (punane joon).

Ülaltoodud jooniselt on selgelt näha, et üksnes keskmiste sünnituste arvu võrdlemine kliinikumi sünnitustubade või perepalatite arvuga ei iseloomusta kogu pilti. Nii näiteks on vaadeldud perioodi keskmine sünnituste arv päevas 6,7, mis tõesti on väiksem kui avatud sünnitustubade arv kaheksa. Samas on veerandil päevadest sünnituste arv 8 või enam ning igal kümnendal päeval on sünnituste arv 10 või enam. (Lisaks esineb veel täiendav varieeruvus tulenevalt sellest, kui kaua üks sünnitaja sünnitustoas on, kuid nii detailseid andmeid kahjuks ei ole avalikult olemas.)

Sünnituste arvu ühes haiglas iseloomustatakse tavaliselt Poissoni jaotusega. (Vt varasemate taoliste analüüside kohta nt Gam et al. A model for the distribution of daily number of births in obstetric clinics based on a descriptive retrospective study. BMJ Open 2013;3: e002920. doi:10.1136/bmjopen-2013-002920.)

Poissoni jaotus tähendab seda, et sünnitajate keskmine arv on seotud ka hajuvusega. Kui sünnitusi tuleb juurde, siis kasvab ka nende varieeruvus. Ka TÜ kliinikumi andmete puhul on sündide arvu dispersioon 7.3 (standardhälve 2.7), mis on ligikaudu võrdne keskmise sünnitajate arvuga, nii nagu nagu Poissoni jaotusele kohane.

Kui oleksin kliinikumi sünnitusosakonna juhataja, siis minule tekitaks Põlva ja Valga haigla sünnitajate lisandumine kahetisi tundeid. Töökoormuse keskmine kasv tähendab keskmiselt suuremat personali vajadust ja seega suuremat eelarvet. Samas kasvab ka sündide varieeruvus, mis tekitab suurema korratuse igapäevatöös. Kasvab päevade arv, kui kõik sünnitustoad on täis, kui perepalateid ei jagu ning ämmaemandad jooksevad mitme sünnitustoa ja operatsiooniblokis tehtavate keisrite vahet.

Näiteks, kui liita 2017. aasta TÜ kliinikumi, Valga ja Põlva haiglate sünnituste arvud, siis kasvaks keskmine sünnituste arv ise vaid 0,8 sünnituse võrra ehk 7,5 sünnituseni päevas. Samas kasvab dispersioon 9,2ni. See toob kaasa selle, et kui kõik Valga ja Põlva haigla sünnitused oleksid Tartus, siis oleks igal neljandal päeval sünnituste arv 10 või enam ning igal kümnendal päeval 12 või enam. Seega ehkki keskmiselt on kõik korras ja sünnituste arv on keskmiselt endiselt väiksem kui sünnitustubade arv, siis kasvab risk, et sünnitajaid on enam kui normaalne töökorraldus võimaldaks.

Sünnituste keskmine arv erineb ka päevade lõikes, kuid iga päeva sees on jaotus taas Poissoni jaotuse lähedane. Nädalavahetuseti on sünnitajate keskmine arv väiksem. (Seda seletab ilmselt plaaniliste keisrite toimumine argitpäeviti.)

Alljärgnevalt on simuleeritud veel sünnitanud naiste sünnitusjärgne viibimine TÜ kliinikumis kui kõik Valga ja Põlva sünnitajad oleksid tulnud TÜ kliinikumi 2017. aastal.
Tervisearenguinstituudi statistika kohaselt Tabel RV301 oli keskmine ravikestus sünnitusabis 2016. aastal 2,8 päeva. Eeldame, et tavasünnitaja on haiglas koos sünnituspäevaga kokku 2, 3 või 4 päeva vastavalt tõenäosustega 50%, 40% ja 10%. Samad tõenäosused keisriga sünnitaja puhul oleksid vastavalt 10%, 20% ja 70%. Kui keisersünnituste osakaal on viiendik, siis saadud osakaalud annavadki keskmiseks voodipäevade arvuks täpselt 2,8.

Simulatsioonitulemused näitavad, et naiste simuleeritud arv, kes on haiglas, ületab 2,5 protsendil juhtudest voodite arvu 30 (13 perepalatit + 16 kohta kahekohalistes palatites + 1 koht ühekohalistes palatites). Eeldatavasti on see puhver, mida saab päevasisese planeerimisega kergesti korraldada, või on see lihtsalt simulatsiooniviga. Kui aga liita siia juurde kõik Valga ja Põlva haigla sünnitajad, siis ulatub naiste simuleeritud arv, kes vooditesse ära ei mahu 8 protsendini. See viitab sellele, et kliinikul on vaja planeerida endale voodi- ja inimreservi võrreldes praeguse olukorraga kui soovida, et kõik sünnitajad kliinikumis end mugavalt tunneks ka tulevikus.

Simuleeritud naiste arv haiglas päevade kaupa (mediaan, ülemine ja alumine kvartiil) on esitatud järgmisel joonisel. Kui sünnitamine oleks planeeritav, siis annaks soovituse noortele, et ajalooliste andmete põhjal on kõige rahulikum sünnitada pühapäeval ja lahkuda kolmapäeva varahommikul.

Kokkuvõtteks sooviks rohkem näha, et taoliste otsuse puhul ei piirdutaks analüüsis vaid paari keskmise arvuga, vaid näidataks (eeldatavasti on kliinikul analüüs tehtud), kuidas mõjutab Põlva ja Valga sulgemine ka kliinikumi ressursside planeerimist tulenevalt kasvavast variatsioonist.

Kui Põlva ja Valga sulgemise puhul tuuakse põhiargumendina madalast sünnituse koguarvust tulenevat võimaliku ohtu teenuse kvaliteedile, siis ei tohiks jätta näitamata, kuidas mõjutab inimeste heaolu ka see, kui sünnitaja satub Tartusse päeval, kui seal on juba ees 15 sünnitajat.

Lõpumärkus. Autoril puudub seos Põlva, Valga või TÜ kliinikumi sünnitusosakonnaga. Kõik kolm last on sündinud TÜ kliinikumis. Antud lühianalüüsi ajend oli loeng loenduvatest andmetest (count data)  ja Poissoni jaotusest majandustudengitele. Kui keegi soovib analüüsi reprodutseerida, siis andmed on all joonistel. Jõudu!

Allikas: mitmed tublid ajakirjanikud, kes fotod ja videod on on-line väljaannetesse pannud

Pilte ametniku punasest Ferrarist

EE „Ratase maksureform: jaht kulakutele alaku!“ seletab maksu- ja õiglushuvilistele üpris tõetruult lahti, mis kedagi meist, sõltuvalt meie sotsiaalsest klassist, ees ootab. Muu hulgas …

„Töötajatel miinus 36 eurot kuus

Järgmine kulakute klass, kellele hoop antakse, on töötajad, kelle brutopalk kuus on 2100 eurot ja enam. Nemad saavad 36 eurot kuus ehk aastas 432 eurot miinust. Kallis Ekspressi lugeja, kui sina ei tea ühtegi inimest, kelle brutopalk on 2100 eurot, siis võid neid tänaval ära tunda punaste Ferraride ja Hermese käekottide järgi. Kui märkad, sisista ja näita nende peale näpuga. Neid inimesi oli 2016. aastal Eestis 58 965 tükki. Kaks korda enam, kui Stalin suutis omal ajal Eestist kulakuid leida, aga samas nii vähe (kui uskuda Kersti Sarapuud ja Keskerakonna pressiteadet), et nendega ei pea arvestama. Teie – 58 965 inimest – võite kõngeda oma hanemaksapasteedi ja teemantsõrmuste kätte!“

Põgus pilk avalike teenistujate avalikele palgaandmetele kinnitab, et tubli kolmandik rahandusministeeriumi ametnikke langeb just sellesse kulakute kategooriasse. Samas pole tähele pannud, et mõni punane auto ministeeriumi parklas oleks Ferrari. Kuid küsimus on huvitav – kuidas tuvastada kulakut Eestis AD 2017, palga või auto järgi? Vastusest sellele küsimusele võiks sõltuda maksupoliitika kujundamine!

Ferrarisid on praegusel aastaajal vähe liikvel, nii et võtame „Eesti asja ajamise auto“ Toyota Land Cruiseri (iseseisvuse algusaastail kutsuti niimoodi hellitavalt 20 a vanust Ford Transitit), mis näeb välja niimoodi ja maksab ümmarguselt 50 000 €. Ka seda laialt levinud keskklassi autot pole ministeeriumi parklas näha, kuid maanteeameti andmeil liigub selliseid kuni viie aasta vanuseid isendeid meie teedel umbes tuhat, kõik värvid kokku. Ja valdavas enamuses kuuluvad need firmadele, mida uuem eksemplar, seda kindlamini.

Nüüd huvitume, kui suurt (maksustatud) palka teenivad inimesed, kes võivad istuda sellise auto roolis (on maanteeametei andmeil auto kasutajad), kui me peaksime kõik need tuhat masinat tänaval kinni pidama. E-riigile kohaselt ei vii me seda sotsiaalset eksperimenti läbi talvise ilmaga tänaval. Maksu- ja maanteeameti koostöös joonistub sellise uuepoolse firmaauto roolis istuvate inimeste palgajaotusest allolev pilt. (Andmed on veidi vananenud – see peegeldab 2013 aasta seisu.)

Lugeda tuleb seda nii, et enim juhte (mood 100 euroste palgavahemike puhul) teenis tolle aja alampalka (300–400 eurot kuus), kaks kolmandikku alla 1200 (ei tea, kas õpetajad aastal 2013 sellisest palgast või autost juba unistasid) ja valdav enamus alla 2100 €. Laiendades veidi automarkide valikut on pilt järgmine. (Endiselt 2013 a palgad, mis praeguseks on mõnevõrra tõusnud.) Ei-ei, see ei ole kriitika ettevõtjate aadressil – (tööjõu)maksude maksmine ei pea olema vabatahtlik!

Seda pilti tuleb tõlgendada nii, et kõigist firmale kuuluva Toyota Yarise roolis olijatest pooled teenivad üle 1000 € kuus maksustatud brutopalka ja pooled alla. Ja ükskõik millise uuepoolse Porsche või Land Cruiseri roolis olijate osas on murdepunkt ehk mediaan veidi üle 900 € bruto.

Mida sellest arvata?

• Miinus 36 eurot kuus ei tapa, kui su palk on 2100, aga …;
• Olulist osa nendest, kelle elustiili ja maitseeelistusi peegeldab uus ilus auto, see maksureform ei puuduta;
• Või tegelikult, puudutab küll – riik kingib suurele osale neist, sarnaselt teiste madalapalgalistega, kuni 62 eurot kuus;
• Valija – taksojuhist ajakirjanikuni – ei kujuta ette, kuidas on Eestis jaotunud maksevõime ja maksukoormus, või see teda ei huvita;
• Ilmselget maksevõimet oleks meie tingimustes efektiivne maksustada autode kaudu (mida on raskem varjata kui tulu), aga see ei näi olevat valijate ja/või ettevõtjate tahe (ehk arusaam õiglusest), mistõttu ei saa see demokraatia tingimustes väljenduda ka poliitilistes otsustes;
• Eelnevast tulenevalt maksustame edaspidigi õllepudelit, mis on aktsiisi mõttes juba praegu kõige kõrgemalt maksustatud ELis ja ilmselt kogu maailmas (sissetulekuid arvestades).

Kas Hegel oleks ka selles kontekstis arvanud, et mis on mõistlik, see on tegelik, ja mis on tegelik, see on mõistlik?

Tere tulemast andmekaeve sajandisse, Statistikaamet! Ehk kuidas uuest andmebaasist R-iga andmeid tõmmata

Viimaks on jõudnud ka Eesti Statistikaamet nii kaugele, et enam ei pea nende kodulehelt andmete alla laadimiseks hiirega kümneid kordi klikkima. Teades andmetabeli nime, on nüüd võimalik uuest andmebaasist andmeid otse alla laadida XML-faili kujul niinimetatud SDMX (Statistical Data and Metadata eXchange) vormingus. Järgnevalt jagan enda kogemusi, kuidas seda teha programmiga R. Soovitan lugeda ka Statistikaameti põgusat juhendit, lk 16) selle teema kohta.

Üks võimalus on andmed automaatselt alla laadida paketiga rsdmx (vt info CRANis või paketi ühe autori Emmanuel Blondeli tutvustavat ettekannet]).

Installeerime ja aktiveerime paketi rsdmx, võtame viimase versiooni githubist, ja teised paketid dplyr ja ggplot2 andmete töötlemiseks ning kuvamiseks.

devtools::install_github("opensdmx/rsdmx")

library(rsdmx)
library(dplyr)
library(ggplot2)

Teades tabeli nime, mida soovime alla laadida, nt IT32: 16-74-aastased arvuti- ja internetikasutajad isikute rühma järgi, võime otse kogu tabeli alla laadida pannes tabeli nime linki GetData nii nagu allosas. Osavamad kasutajad võivad lingi täpsema kuju abil ka kohe piirata tunnuste ja väärtuste hulka, mida alla laadida, aga mulle tundus, et lihtsam on kogu tabel alla laadida ja pärast R-s valida tunnuseid ja väärtuseid.

#Määrame ajutise faili tf, kuhu veebilehe alla laeme juhuks, kui suur fail
tf <- tempfile(tmpdir = tdir <- tempdir())
download.file("http://andmebaas.stat.ee/restsdmx/sdmx.ashx/GetData/IT32/all", tf)
#loeme andmed ajutisest failist sisse
temp_sdmx <- readSDMX(tf, isURL = FALSE)
#teeme andmetabeliks
andmetabel <- as.data.frame(temp_sdmx)
str(andmetabel)

## 'data.frame':    2000 obs. of  6 variables:
##  $ DIM2       : chr  "1" "1" "1" "1" ...
##  $ DIM3       : chr  "1" "1" "1" "1" ...
##  $ DIM4       : chr  "1" "1" "1" "1" ...
##  $ TIME_FORMAT: chr  "P1Y" "P1Y" "P1Y" "P1Y" ...
##  $ obsTime    : chr  "2005" "2006" "2007" "2008" ...
##  $ obsValue   : num  624 641 670 674 726 ...

Esialgne valmis tabel sisaldab lisaks väärtustele tunnuses obsVAlue tunnuste nimede asemel mittemidagiütlevaid dimensioone DIM2, DIM3, …, ja nende väärtuste tähistusi. Vajalikud tähistused saame alla laadida kombineerides tabeli nime teistsuguse lingiga GetDataStructure . Tunnuste DIM2, DIM3, …, selgitused saame allalaetava objekti andmestruktuur väljalt concepts.

andmetestrukt <- readSDMX("http://andmebaas.stat.ee/restsdmx/sdmx.ashx/GetDataStructure/IT32")
#tunnuste selgitused
tunnused <- as.data.frame(andmetestrukt@concepts)
#et täpitähti õigesti näitaks
tunnused$Name.et <- iconv(tunnused$Name.et, "UTF-8")
tunnused

##            id agencyID                 Name.en                     Name.et
## 1        DIM2  STAT_EE    Group of individuals                Isikute rühm
## 2        DIM3  STAT_EE Using computer/Internet Arvuti/Interneti kasutamine
## 3        DIM4  STAT_EE               Indicator                     Näitaja
## 4        TIME  STAT_EE                    Year                       Aasta
## 5   OBS_VALUE  STAT_EE                                            
## 6 TIME_FORMAT  STAT_EE                                            
## 7  OBS_STATUS  STAT_EE                                            
##                   en
## 1               
## 2               
## 3               
## 4               
## 5  Observation Value
## 6        Time Format
## 7 Observation Status

Nüüd näeme, et DIM2 on pikalt kirjutades “Isikute rühm”  DIM3 on “Arvuti/Interneti kasutamine” jne.

Tunnuste kodeeringute väärtused saame sama andmetestrukt objekti väljadelt codelists. Nende tähised on kujul, mis sisaldab tabeli ja dimensiooni nime, nt “CL_IT32_DIM2” jne

#tunnuste väärtuste nimed, lähme kaks korda sisse välja codelists ja võtame sealt id-d
margenditelist <- slot(andmetestrukt, "codelists")
margenditenimed <- sapply(slot(margenditelist, "codelists"), function(x) slot(x, "id"))
margenditenimed

## [1] "CL_IT32_DIM2"        "CL_IT32_DIM3"        "CL_IT32_DIM4"       
## [4] "CL_IT32_TIME"        "CL_IT32_OBS_STATUS"  "CL_IT32_TIME_FORMAT"

Konkreetsed väärtused saame andmetabelitesse või lisatunnustena soovi korral juurde siduda. Pakun all ühe variandi, eks nupukamad R-i kasutajad ilmselt oskavad seda lihtsamalt teha.

#kirjutame märgendid tabelitesse
for(i in margenditenimed) {
  #paneme tabelite nimedeks muutujatega samad nimed
  assign(substr(i,nchar("CL_IT32_")+1, nchar(i)), as.data.frame(slot(andmetestrukt, "codelists"), codelistId = i))
}
#ja seome kõigile DIM-algusega tunnustele juurde eesti- ja inglise keelsed märgendid
for(i in grep("DIM", names(andmetabel), value=TRUE)) {
  andmetabel <- merge(andmetabel, get(i), by.x=i, by.y = "id" )
  names(andmetabel)[names(andmetabel) == 'label.en'] <- paste0(i,'label.en')
  names(andmetabel)[names(andmetabel) == 'label.et'] <- paste0(i,'label.et')
  }

Ongi tabel valmis koos kõigi märgenditega ja võime andmeid kasutama hakata. Näiteks joonistame internetikasutajate osakaalu vanusrühmade lõikes:

andmetabel %>% filter(DIM3label.et=="Interneti kasutamine",
                      DIM4label.et=="Osatähtsus vastavas rühmas,%",
                      grepl("aastased", DIM2label.et)) %>%
  ggplot(aes(x=as.numeric(obsTime), y=obsValue, color=DIM2label.et)) +
  geom_line() +
  scale_x_continuous(breaks = seq(2005, 2017, 2)) +
  labs(y = "Interneti kasutajate osakaal, %",
       x = "",
       colour = "",
       title = "Interneti kasutajate osakaal eri vanusrühmades",
       subtitle = "Andmed: Statistikaamet, tabel IT32") +
  theme(text=element_text(size=14))

Analoogselt võime võtta mõne muu tabeli, näiteks leibkondade kulutuste struktuuri tabelist LE207: Leibkonnaliikme kulutused aastas leibkonnapea haridustaseme järgi. Dubleerides ülal tehtut, kuvame alkoholi- ja tubakatoodete kulutuste keskmise osakaalu leibkonna kulude struktuur leibkonna pea haridusetaseme järgi.

tf <- tempfile(tmpdir = tdir <- tempdir())
download.file("http://andmebaas.stat.ee/restsdmx/sdmx.ashx/GetData/LE207/all", tf)
andmetabel <- readSDMX(tf, isURL = FALSE) %>% as.data.frame()
andmetestrukt <- readSDMX("http://andmebaas.stat.ee/restsdmx/sdmx.ashx/GetDataStructure/LE207")
tunnused <- as.data.frame(andmetestrukt@concepts)
tunnused$Name.et <- iconv(tunnused$Name.et, "UTF-8")
margenditelist <- slot(andmetestrukt, "codelists")
margenditenimed <- sapply(slot(margenditelist, "codelists"), function(x) slot(x, "id"))
for(i in margenditenimed) {
  assign(substr(i,nchar("CL_LE207_")+1, nchar(i)), as.data.frame(slot(andmetestrukt, "codelists"), codelistId = i))
}
for(i in grep("DIM", names(andmetabel), value=TRUE)) {
  andmetabel <- merge(andmetabel, get(i), by.x=i, by.y = "id" )
  names(andmetabel)[names(andmetabel) == 'label.en'] <- paste0(i,'label.en')
  names(andmetabel)[names(andmetabel) == 'label.et'] <- paste0(i,'label.et')   
} 

andmetabel %>% filter(DIM2label.et=="Osatähtsus (kulutused kokku = 100), %", DIM3label.et=="Alkoholjoogid ja tubakas") %>%
  ggplot(aes(x=as.numeric(obsTime), y=obsValue, color=DIM4label.et)) +
  geom_line() +
  geom_point()+
  #scale_x_continuous(breaks = seq(2010, 2017, 1)) +
  labs(y = "Osakaal, %",
       x = "",
       colour = "",
       title = "Alkoholi ja tubakatoodete kulutuste osakaal leibkondade eelarves \n leibkonnapea hariduse järgi",
       subtitle = "Andmed: Statistikaamet, tabel LE207") +
  theme(text=element_text(size=14))

Jõudu tööle!

Kuidas TAIst kergesti andmed Ri saada?

Eestis kogub tervise valdkonna statistikat TAI – Tervise Arengu Instituut, kes avaldab info enda tervisestatistika ja terviseuuringute andmebaasis. Nende andmebaas kasutab programmi PC-Axis, mis on välja töötatud Rootsi Statistikaameti poolt.
Sellele tarkvarale on kirjutatud R-is valmis pakett pxweb (vt info CRANis), mis võimaldavad neid andmeid otse lugeda R-i kas otse käsuga või pool-interaktiivsel moel.

Installeerime ja aktiveerime paketi pxweb ning aktiveerime ka paketid dplyr ja ggplot2 andmete töötlemiseks ning kuvamiseks.

install.packages("pxweb")

library(pxweb)
library(dplyr)
library(ggplot2)

Paketiga tuleb kaasa fail tekstifail api.json, mis asub paketi kataloogis /pxweb/extdata/api.json. Failis on kirjas asutused, mis teadaolevalt samuti kasutavad tarkvara PC-AXIS ja võimaldavad andmeid alla laadida. Neid saab loetleda funktsiooniga api_catalogue. Vaatame näiteks api kataloogi kolme esimest elementi, mis on Rootsi ja Soome statistikaametite lehed.

api_catalogue()[1:3]

## [[1]]
## Api: api.scb.se ('scb')
##      Statistics Sweden 
## Version(s)   : v1 
## Language(s)  : en, sv 
## Limit(s)     : 30 calls per 10 sec.
##                1e+05  values per call.
## Url template :
##  http://api.scb.se/OV0104/[version]/doris/[lang] 
## 
## [[2]]
## Api: pxnet2.stat.fi ('statfi')
##      Statistics Finland 
## Version(s)   : v1 
## Language(s)  : fi, sv, en 
## Limit(s)     : 30 calls per 10 sec.
##                1e+05  values per call.
## Url template :
##  http://pxnet2.stat.fi/PXWeb/api/[version]/[lang] 
## 
## [[3]]
## Api: pxwebapi2.stat.fi
##      Statistics Finland (old version) 
## Version(s)   : v1 
## Language(s)  : fi 
## Limit(s)     : 30 calls per 10 sec.
##                1e+05  values per call.
## Url template :
##  http://pxwebapi2.stat.fi/PXWeb/api/[version]/[lang]

Selleks, et pakett töötaks ka TAI andmetega, tuleb lisada faili api.json TAI andmete asukoht ja vajalikud parameetrid. Selleks on nüüd kaks võimalust:

a) lisada käsitsi juurde faili api.json TAI andmed suvalise tekstiredaktoriga.
b) või kirjutame faili api.json üle, kui on kindel, et me teiste riikide andmeid kasutada ei soovi.

Alljärgnevalt kasutan teist varianti: nimetan orignaalfaili ümber ja kirjuta uue faili vaid TAI andmetega.

file.rename("~/R/win-library/3.4/pxweb/extdata/api.json", 
"~/R/win-library/3.4/pxweb/extdata/api_orig.json")

tai_api <- 
  '{
  "apis": {
     "pxweb.tai.ee": {
      "alias" : ["tai_et"],
      "description" : "Estonian National Institute for Health Development",
      "url" : "http://pxweb.tai.ee/PXWeb2015/api/[version]/[lang]/",
      "version": ["v1"], 
      "lang": ["en", "et"],
      "calls_per_period": 10,
      "period_in_seconds": 10,
      "max_values_to_download" : 100000
     }
  },

  "local_apis": {}
}'
write(tai_api, file = "~/R/win-library/3.4/pxweb/extdata/api.json")

Andmete esmakordseks allalaadimiseks, kui veel süntaksi kuju ei ole tuttav, siis on soovitatav kasutada interaktiivset käsku interactive_pxweb().

d <- interactive_pxweb()

See juhendab samm-haaval kuni lõpliku andmepäringuni, mis tehakse funktsiooniga get_pxweb_data. Lõpliku andmepäringu kuju saab seejärel ise kergesti edasi töödelda.

Kui tabeli nimetus ja soovitud tunnuste väärtused on varasemalt teada, siis saab kirjutada andmepäringu juba otse. Näiteks laeme alla tervishoiutöötajate palgad kõikide aastate, ametialade, teenuse osutaja liigi ja palga liigi lõikes. Ma tean, et see asub TAI andmebaasi alamkataloogis 04THressursid/06THTootajatePalk, tabelis TT06.px.

df_long <- 
  get_pxweb_data(
    url = "http://pxweb.tai.ee/PXWeb2015/api/v1/et/04THressursid/06THTootajatePalk/TT06.px",
             dims = list(Aasta = c('*'),
                         Ametiala = c('*'),
                         "Tervishoiuteenuse osutaja liik" = c('*'),
                         "Brutokuupalga liik" = c('*')),
             clean = TRUE)

Ülal päringus tähistab tärn *, et soovime kõiki väärtuseid antud tunnuse puhul. Teine võimalus on kirjutada sinna sisse konkreetne väärtuste loetelu. Valik clean võimaldab valida, kas soovitakse andmeid pikas formaadis (long), kus iga väärtus on eraldi reas, või laias formaadis (wide), kus iga väärtus on eraldi veerus. clean = TRUE annab andmed pikas formaadis.
Kui sooviksin andmeid laias vormis, siis paneksin lihtsalt clean = FALSE

df_wide <- 
  get_pxweb_data(
    url = "http://pxweb.tai.ee/PXWeb2015/api/v1/et/04THressursid/06THTootajatePalk/TT06.px",
             dims = list(Aasta = c('*'),
                         Ametiala = c('*'),
                         "Tervishoiuteenuse osutaja liik" = c('*'),
                         "Brutokuupalga liik" = c('*')),
             clean = FALSE)

Ja ongi andmed alla laetud ja valmis töötlemiseks. Näiteks esitame allalaetud andmete põhjal tervishoiutöötajate mediaanpalkade arengu 2011-2017.

df_long %>%  
  filter(`Tervishoiuteenuse osutaja liik` == "Teenuseosutajad kokku", 
         `Brutokuupalga liik` == "Kuupalk koos kõigi lisatasudega, mediaan") %>% 
  ggplot(aes(y=values, x=Aasta, group=Ametiala, color=Ametiala)) +
  geom_point()  +
  geom_line() + 
  labs(title = "Mediaanpalga areng", color="", x="", y="Euro", 
       caption = "Allikas: TAI, tabel TT06")

Allikad:

Mans Magnusson, Leo Lahti and Love Hansson (rOpenGov 2014). pxweb: R tools for PX-WEB API. URL: http://github.com/ropengov/pxweb