Commit 35b51d6a authored by patrick.lambert's avatar patrick.lambert
Browse files

thermal range

parent 91214804
library(tidyverse)
#TODO extract from xml
thermalRange = tibble(process = "Grow", Tmin = 1.7, Topt = 4.5, Tmax = 27.9, type = "Rosso") %>%
bind_rows(tibble(process = "SurvivalSpawnerInRiv", Tmin = 5.4, Topt = 16.7, Tmax = 27.5, type = "Rosso") ) %>%
bind_rows(tibble(process = "Reproduction", Tmin = 5.1, Topt = 13.5, Tmax = 24.5, type = "Rosso") ) %>%
bind_rows(tibble(process = "SurviveAfterReproduction", Tmin = -Inf, Topt = 19.58 - log(19)/0.58, Tmax = 19.58 + log(19)/0.58, type = "logit") )
thermalRange %>% mutate(numero = row_number(),
process = factor(process,
levels = c("Grow","SurvivalSpawnerInRiv", "SurviveAfterReproduction", "Reproduction" ))) %>%
ggplot() +
geom_segment(aes(x = Tmin, y = process, xend = Tmax, yend = process)) +
geom_point(aes(x = Topt, y = process))
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment