# NOT RUN {
library("tidyverse")
time_profile <- seq(0, 50, length = 20)
Temp_profile <- seq(52.5,60, length = 3)
exp_design <- expand.grid(time_profile,Temp_profile) %>%
rename(times = Var1, temperature = Var2)
pars <- list(z = 4.2, D_R = 3.9, temp_ref = 55)
my_sensitivities <- isothermal_sensitivities("Bigelow", exp_design, pars)
plot(my_sensitivities)
plot(my_sensitivities, limit = 6)
# }
Run the code above in your browser using DataLab