Prdaragadoz interakcik modellezsnek lehetsgei Rben Hogyan hasznljuk az
Préda-ragadozó interakciók modellezésének lehetőségei R-ben
Hogyan használjuk az R-t? Mint más emberek: - adatelemzés - statisztika - grafika n Modellezés R-ben: – alkalmazott modellek – „magyarázó” modellek – demostráló modellek és modellezés n
Ökológiai modellek n n - n Egy ökológiai rendszer vagy alrendszer matematikai modellje Folyamat (process) modellek: rendszer megismerés, hipotézis tesztelés (préda-predátor modell, sejtautomaták) Szintézis (synthetic) modellek: ismeretek összefoglalása, döntéshozás (vízminőség, globális folyamatok)
A feladat(ok)…. n n n 1. csomagok letöltése, telepítése - A csomagok (simecol, odesolve) 2. A scriptek letöltése, értelmezése 3. Néhány példa kiszámítása 4. jegyzőkönyv elkészítése 5. jegyzőkönyvek elküldése (*. pdf-ben) az elek. zoltan@aotk. szie. hu címre 6. elindulni lassan hazafelé…
1. modell n n n n library(odesolve) library(simecol) #####basic Lotka-Volterra predator prey-model # k 1 -growth rate of the prey population. #k 2 -encounter rate of predator and prey. #k 3 -death rate of the predator population. data(lv) print(lv) parms(lv) <- list(k 1=0. 8, k 2=0. 1, k 3=0. 9) times (lv)<- c(from=0, to=300, by=0. 6) init(lv) <- c(prey=1, predator=5) print(lv) plot(sim(lv))
Példa - Lotka-Volterra-alapmodell: init(lv) <- c(prey=10, predator=5) L 1 L 2 L 3 L 4 L 5 L 6 L 7 L 8 Préda b. (k 1) Nagy (0. 9) Kicsi (0. 1) Kicsi(0. 1) Interakt. (k 2) Kicsi (0. 1) Nagy (0. 9) Nagy(0. 9) Ragadozó d. (k 3) Nagy (0. 9) Kicsi (0. 1) Nagy (0. 9) Kicsi(0. 1)
Példa - Lotka-Volterra-alapmodell: init(lv) <- c(prey=5, predator=10) L 1 L 2 L 3 L 4 L 5 L 6 L 7 L 8 Préda b. (k 1) Nagy (0. 9) Kicsi (0. 1) Kicsi(0. 1) Interakt. (k 2) Kicsi (0. 1) Nagy (0. 9) Nagy(0. 9) Ragadozó d. (k 3) Nagy (0. 9) Kicsi (0. 1) Nagy (0. 9) Kicsi(0. 1)
2. modell n n n #####predator prey-model with three equations: predator, prey and resource (e. g. nutriens, grassland) library(odesolve) library(simecol) data(lv 3) plot(sim(lv 3)) solver(lv 3) <- "lsoda" parms(lv 3)<- list (b=0. 1, c=0. 5, d=0. 1, e=0. 4, f= 0. 4, g=0. 1) init(lv 3)<- list (s=0. 8, p=0. 2, k=0. 0) times(lv 3)<- c(from=0, to=100, by=0. 1), print(sim(lv 3)) plot(sim(lv 3)) # wrong! integration overlooks internal inputs plot(sim(lv 3, hmax=1)) # integration with correct maximum time step
n http: //www. univet. hu/users/zelek/okologia-III/ n www. simecol. de n T. J. Case (2000): An Illustrated Guide to Theoretical Ecology. Oxford University Press, New York.
- Slides: 9