Eksponentialfunksjonen som modell. Regresjon
Talet på aure i eit vatn har auka kraftig etter at det i 1998 vart sett i gang med kalking av vatnet. Tabellen viser talet på aure i vatnet nokre år etter 1998.
Årstal | 1998 | 2000 | 2002 | 2004 |
---|---|---|---|---|
Talet på år etter 1998, | 0 | 2 | 4 | 6 |
Talet på aure i tusen, | 4,0 | 6,7 | 10,9 | 17,4 |
Tabellen viser at aurebestanden aukar meir og meir. Vi vil finne ein modell for utviklinga av aurebestanden.
Vi legg dataa frå tabellen inn i reknearket i GeoGebra. Så merker vi cellene og klikkar på knappen for regresjonsanalyse.
Punkta viser at vi må finne ein funksjon som veks raskare og raskare ettersom -verdiane aukar, noko vi òg ser direkte av tabellen. Derfor kan ein eksponentiell modell passe godt med dei observerte verdiane.
Ved å velje "Eksponentiell 2" som regresjonsmodell får vi talet som grunntal i potensen.
Vi får funksjonen gitt ved
som modell for utviklinga av aurebestanden.
For å få grafen og punkta over i det vanlege grafikkfeltet kopierer vi han over ved å trykke på knappen til høgre for innstillingshjulet, og vi vel "Kopier til grafikkfeltet".
Av grafane kan vi til dømes sjå at ut frå denne modellen vil aurebestanden ha passert 46 000 individ i 2008 og 76 000 individ i 2010.
Regresjon med Python kan gjerast på fleire måtar. Her har vi valt å bruke metoden curvefit
frå scipy.optimize
sjølv om han gir oss meir data enn vi skal bruke her.
Metoden er basert på at vi skriv i koden kva slags type funksjon som skal brukast i den matematiske modellen. Dette blir gjort ved å definere modellen som ein eigendefinert pythonfunksjon i koden. I vårt tilfelle ønsker vi å tilpasse målingane til ein eksponentialfunksjon av typen
Parameterane til pythonfunksjonen modell
er først den frie variabelen x
, deretter konstantane som skal brukast i modellfunksjonen, modell
blir brukt av metoden curvefit
saman med ei liste med curvefit
gir tilbake to lister. Den eine inneheld dei to konstantane a
og b
(i vårt tilfelle) som gjer at den matematiske funksjonen vi definerte i pythonfunksjonen modell
, passar mest mogleg med måledataa. Den andre lista inneheld dei såkalla kovariansane, som vi ikkje skal bruke her.
Koden kan vidare sjå slik ut:
Køyring av koden gir denne utskrifta:
"Funksjonen blir f(x) = 4.126e(0.240x)."
Vi får nesten den same modellfunksjonen som med GeoGebra. Årsaka til at det ikkje blir likt, kan vere at ulike program kan bruke litt ulike måtar å rekne seg fram til resultatet på.
Etter modellen vil aurebestanden halde fram med å vekse raskare og raskare. Er det slik det vil utvikle seg vidare?
Registreringa av aurebestanden heldt fram kvart andre år etter 2004. Resultata ser du i tabellen nedanfor.
Årstal | 1998 | 2000 | 2002 | 2004 | 2006 | 2008 | 2010 |
---|---|---|---|---|---|---|---|
Talet på år etter 1998, | 0 | 2 | 4 | 6 | 8 | 10 | 12 |
Talet på aure i tusen, | 4,0 | 6,7 | 10,9 | 17,4 | 21,5 | 24,5 | 26,0 |
I koordinatsystemet nedanfor ser du grafen til
Det viser seg altså at modellen ikkje var eigna til å seie noko om utviklinga i perioden etter 2004.
Eksponentialfunksjonar passar ofte godt til å beskrive korleis populasjonar forandrar seg i eit avgrensa tidsrom mens det er rik tilgang på mat og lite eller ingenting som avgrensar veksten.
Etter kvart blir utviklinga annleis. Når populasjonen blir stor nok, blir veksten bremsa fordi det til dømes blir for lite mat.
På teorisida "Logistisk vekst" ser vi på ein matematisk modell som passar betre til utviklinga i aurebestanden.