Korona – Tartunnat Suomessa ja kuolleisuus koko maailmassa – GeoGebra -kirja

Kiinnitän tämän artikkelin blogini alkuun joksikin ajaksi. Päivitän tätä sivua muutaman päivän välein. Ilmoitan päivityksistä yms. Tarinaan liittyvä GeoGebra-kirja https://www.geogebra.org/m/v6jqtm9w

Tällä hetkellä GeoGebra-kirjassa on taulukot ”Kuolleisuus koko maailmassa” ja ”Tartunnat Suomessa”.

THL:n luvut julkaistaan puolen päivän jälkeen, joten päivitän appeja iltapäivällä.


31.3.

Päivitin maailman luvut. Kasvukerroin ei jurikaan muuttunut. Loin uuden Suomen kuolleista, en saa sitä vielä lisättyä GG-kirjaan, pitää odottaa, että palvelin päivittyy. En laittanut sinne nollia, näin saan päivitettyä sovitukset helpommin (nollasta ei voi ottaa logaritmia). https://www.geogebra.org/m/yjqrpgpu En päivitä enää Suomen tartuntalukuja, niistä kun on niin ristiriitaisia lukuja.

28.3.

Päivitin koko maailman kuolinluvut. Nyt kuolleita on 3271. Kasvukerroin on noin 1,17. Pohdin vielä mitä teen Suomen luvuilla. THL:n sivuilta on haastavaa löytää tietoa aiheesta. Nytpä he ovatkin muuttaneen avoimen datansa raja pinnan, pitänee opiskella tuota.

17.3.

Lisäsin eilisen luvut taulukkoihin, maailmalla 641 kuollutta Suomessa 47 uutta todettua tapausta. Nyt minun summani (318) poikkeaa yhdellä THL:n luvuista (319).

16.3.

Eilinen tartuntamäärä Suomessa oli 31. Näissä luvuissa näkyy testaamisessa tapahtunut muutos. Tällä hetkellä eksponentiaalisen kasvun päivittäinen kasvukerroin näyttäisi olevan noin 1,3. Maailman laajuinen kuolleisuus nousi tosi paljon. Eilen kuolleita oli 686.

15.3.

Päivitin luvut. Suomessa tartuntojen määrä pieneni 30:een. Syy taitaa olla muutoksissa miten testejä tehdään esim Helsingissä ja Uudellamaalla.

14.3.

Tein Suomen tartunnat appletin kokonaan uusiksi. Luvut otin THL:n Ajankohtaista-sivulta. Valitsin pisteitä alkaen päivämäärästä 3.3. Silloin oli ensimmäisen kerran tartuntoja yli kaksi. Tein näihin pisteihin kasvusovituksen ja suoran sovituksen linlog-pisteihin. Pisteet näkyvät kuvissa sinisinä. Tällä hetkellä eksponentiaalisen kasvun päivittäinen kasvukerroin näyttäisi olevan noin 1,4.

Vihdoin sain omat lukuni täsmäämään THL:n lukuihin. Lisäsin sen ensimmäisen tapauksen Lapista taulukkooni ja korjasin pari virhettä. Tosin THL:n Ajankohtaista uutisten yksi luku ei vastaa saman päivän arvoa heidän omassa kuvaajassaan.

13.3.

Päivitin tämän päivän luvut. Käytin THL:n ilmoittamia lukuja. Suomen tilanteesta tulen käyttämään THL:n sivuilla olevia lukuja. Kumma, kun he eivät jaa taulukoitua dataa omilla sivuillaan. Noissa GitHub-luvuissa ja Suomessa esitetyissä luvuissa on pieniä poikkeamia.

12.3.

Loin uuden taulukon aiheesta koronatartunnat Suomessa. Loin samalla uuden GeoGebra-kirjan. Siinä on tällä hetkellä tuo appletti sekä aiemmin tekemäni materiaali koronakuolleista koko Maailmassa.

https://www.geogebra.org/m/v6jqtm9w

Pyrin päivittämään näitä vähintään viikottain, joskus useamminkin.

Lue aiheeseen liittyvät aiemmat tarinat

Koronavirus – mallintamista GeoGebralla

Miksi kuvaajissa esitetään tietoa logaritmisella asteikolla?

Miksi kuvaajissa esitetään tietoa logaritmisella asteikolla?

[edit. 21.2.20 Korjasin virheellisen kuvaajan.]
[edit. 22.2.20 Korjasin muutaman fysikaalisen virheen, kiitos Harri. Vaihdoin Worldometer-kuvan.]

Kun menet jollekin sivulle, joka esittää tietoa kuvaajina, niin usein törmäät tilanteeseen, että kuvaajan asteikot ovatkin logaritmisia. Miksi? Mitä koronaviruksen etenemistä esittävillä sivustoilla kuolleiden määrän logaritmi kertoo?

Minusta yksi kauneimmista tilastollisista esityksistä on Hertzsprung–Russell -diagrammi. Tähtitieteilijöillä on siitä monia eri versioita, mutta perinteisesti vaaka-akselilla on tähden pinnan väri sinisestä punaiseen ja pystyakselilla tähden absoluuttinen luminositeetti eli säteilyteho. Koska tähden pintalämpötila on sen pintalämpötilan mitta, niin loppujen lopuksi vaaka-akselilla onkin pintalämpötilan logaritmi oikealta vasemmalle eli mitä kuumempi tähti, niin sitä enemmän vasemmalla se on. Pystyakseli luminositeetti eli säteilyteho kertoo, kuinka moninkertainen tähden säteilyteho on Auringon tehoon (3,9·10^26 W) verrattuna. HR-diagrammissa pystyakseli on suhteellisen tehon logaritmi. Näin Aurinkomme sattuu mukavasti diagrammin keskialueelle.

HR-diagrammi. Tähtien luminositeetti pintalämpötilan funktiona. ESA:n kuva
In the Hertzprung-Russell diagram the temperatures of stars are plotted against their luminosities. The position of a star in the diagram provides information about its present stage and its mass. Stars that burn hydrogen into helium lie on the diagonal branch, the so-called main sequence. Red dwarfs like AB Doradus C lie in the cool and faint corner. AB Dor C has itself a temperature of about 3,000 degrees and a luminosity which is 0.2% that of the Sun. When a star exhausts all the hydrogen, it leaves the main sequence and becomes a red giant or a supergiant, depending on its mass (AB Doradus C will never leave the main sequence since it burns so little hydrogen). Stars with the mass of the Sun which have burnt all their fuel evolve finally into a white dwarf (left low corner).
lähde ESO: https://www.eso.org/public/images/eso0728c/

Kuuloalueen kuvaaja on fysiikan opiskelijoille tuttu kuvaaja MAOL taulukot -kirjassa. Siinä vaaka-akselilla on taajuus logaritmisella asteikolla ja pystyakseli on äänen intensiteettitaso desibeleinä, joka on suoraan verrannollinen äänen tehosuhteen logaritmiin.

kuuloalue kuvaaja
lähde: MAOL taulukot. Otava-MAOL. 2020

Gapminder sivustolla on valtavasti eri valtioihin liittyvää tilastotietoa. Aineistojen aikasarjoja voi tutkia myös logaritmisilla asteikoilla. Kuvassa on vaaka-akselilla BKT/asukas eri maissa ja pystyakselilla odotettavissa oleva elinikä. Itse asiassa myös valtioita kuvaavien ympyröiden säde on logaritminen.

Kuva, joka sisältää kohteen kartta, teksti

Kuvaus luotu automaattisesti

Lähde Gapminder sivu https://www.gapminder.org/tools/#$state$time$value=2018;;&chart-type=bubbles

Worldometer-sivustolla voi seurata korona-viruksen kuolleiden kokonaismäärää logaritmisella asteikolla.

Koronavirukseen kuolleiden ja sairastuneiden määrät
Lähde: https://www.worldometers.info/coronavirus/

polynomiaalinen ja eksponentiaalinen kasvu

Tutkitaan esimerkkien avulla, miten akseleiden muuttaminen näkyy polynomeilla ja eksponenttifunktiolla. Käytän tässä tarinassa 10-kantaista logaritmia, sillä sitä käytetään yleensä kansantajuisissa taulukoissa. Tällöin näkee alkuperäisten lukujen suuruusluokan mukavasti. 10-kantainen logaritmi kun kertoo luvun nollien lukumäärän tyyliin lg(1000) = 3.  

Lukiomatikassa tutkitaan yleensä aritmeettista tai geometrista lukujonoa. Aritmeettisessa lukujonossa kahden lukujonon peräkkäisen erotus pysyy vakiona ja geometrisessä suhde. Jos lukujonoja ajattelee funktioina siten, että järjestysluku on muuttuja x, niin aritmeettinen lukujono on 1. asteen polynomi f(x) = k x + b ja geometrinen lukujono on eksponenttifunktio g(x) = a·q^x.

On suhteellisen helppoa osoittaa, että jos k, n, a ja q > 1, niin a·q^x > k x^n kunhan x on riittävän iso. Eli eksponentiaalinen kasvu ohittaa aina jossain vaiheessa polynomiaalisen kasvun.

Kun pohditaan kasvunopeutta, niin ei yleensä välitetä polynomien a….. muista kuin korkeimman asteen termistä. Yleensä myös yleisen potenssifunktion tyyliin 2.5  x^1,23 kasvua kutsutaan polynomiseksi, vaikka 1,23 ei olekaan luonnollinen luku (0, 1, 2, 3, …).

Otetaan esimerkiksi funktiot f(x) = 3*x^2 ja g(x) = 3*2^x. Funktio f on toisen asteen polynomi ja g eksponenttifunktio, jonka kasvukerroin (taulukkokirjassa talousmatematiikan kohdalla korkotekijä) q = 2.

loglog

Luon GeoGebrassa taulukon, jossa kaksi ensimmäistä saraketta ovat x ja f(x), seuraavat kaksi x ja lg(f(x)) ja kolmannessa parissa lg(x) ja lg(f(x)) Luon noista kolme pistelistaa nimetä linlin, linlog ja loglog. Punainen pistelista lililin on normaali ”koordinaatisto”, vihreässä on normaali x ja y-koordinaatti on lg(f(x)). Sinisessä loglog-pistelistassa molemmista on otettu logaritmi.

Kuvaajasta nähdään, että vihreä lukujono kasvaa todella hitaasti, mutta sininen näkyisi kasvavan suoraviivaisesti. Kun sinisiin pisteisiin sovitetaan suora komennolla

p(x) = SovitaPolynomi(loglog, 1)

havaitaan, että suoran kulmakerroin on 2 eli alkuperäisen polynomimme asteluku. Vakiotermin avulla päästään kertoimeen laskemalla CASissa

10^0,47712
-> 2,99999

Tai jos olet lukenut edellisen artikkelini ”Arvojen poimiminen listoista ja yhtälön ratkaisuista GeoGebrassa

10^Alkio(Kertoimet(p), 2)
-> 3

Jätän lukijalle todistettavaksi, että jos f(x) = k x^n, niin loglog-koordinaatistossa syntyy suora, jonka kulmakerroin on n ja k = 10^(suoran vakiotermi).

linlog

Tutkitaan samalla tavalla eksponenttifunktiota

g(x) = 3*2^x

Nyt vihreä loglog-käyrä ei ole suora, mutta linlog on. Sen yhtälöstä saadaan

10^0.30103
-> 2
10^0.47712
-> 3

Jätän taas lukijalle pohdittavaksi miksi eksponenttifunktiolla a q^x loglog-sovituksen suoran yhtälöllä q = 10^kulmakerroin ja a = 10^vakiotermi.

Lisäksi voit pohtia, minkä tyyppinen käyrä sininen loglog lukujono on.

lopuksi

Toki lukusarjoja olisi voinut alun perikin tutkia käyttämällä GeoGebran sovita-komentoja. Silti monessa tapauksessa, kun tutkitaan lukusarjoja, joiden kasvunopeuksia ei tunneta, on logaritmisten akseleiden käyttö aika kätevää.

Vielä muistutukseksi.

Jos linlog- koordinaatistossa on suora, niin se kertoo eksponentiaalisesta  a·q^x kasvusta.

Jos loglog -koordinaatistossa on suora, niin se kertoo polynomiaalisesta (potenssifunktio) k·x^n kasvusta.

Mitä Worldometer-kuvaaja kertoo? Koronavirukseen sairastuneiden ja kuolleiden kokonaismäärä ei tämän datan perusteella kasva eksponentiaalisesti. Palannen tähän aiheeseen, kun käytössä on pidempi aikasarja ja useammasta eri valtiosta.