Commit c697a18b authored by patrick.lambert's avatar patrick.lambert
Browse files

new temperature effect function

parent 4ddfe864
......@@ -185,6 +185,13 @@ modifiedMasterReaction = function(temperature, par){
return(rate)
}
# DeucthEtal
# Deutsch, Curtis A., Joshua J. Tewksbury, Raymond B. Huey, Kimberly S. Sheldon, Cameron K. Ghalambor, David C. Haak, et Paul R. Martin. « Impacts of Climate Warming on Terrestrial Ectotherms across Latitude ». Proceedings of the National Academy of Sciences 105, nᵒ 18 (5 juin 2008): 6668‑72. https://doi.org/10.1073/pnas.0709472105.
# Huey, Raymond B., et R. D. Stevenson. « Integrating Thermal Physiology and Ecology of Ectotherms: A Discussion of Approaches ». American Zoologist 19, nᵒ 1 (1 février 1979): 357‑66. https://doi.org/10.1093/icb/19.1.357.
# Payne, Nicholas L., James A. Smith, Dylan E. van der Meulen, Matthew D. Taylor, Yuuki Y. Watanabe, Akinori Takahashi, Teagan A. Marzullo, Charles A. Gray, Gwenael Cadiou, et Iain M. Suthers. « Temperature Dependence of Fish Performance in the Wild: Links with Species Biogeography and Physiological Thermal Tolerance ». Functional Ecology 30, nᵒ 6 (2016): 903‑12. https://doi.org/10.1111/1365-2435.12618.
temperature = seq(0,100,.1)
plot(temperature, modifiedMasterReaction(temperature, par = c(48.6, 49.7, 388, 1)), type = 'l', ylab = 'rate', main = "modifiedMasterReaction",
......@@ -250,3 +257,17 @@ plot(temperature, modifiedThorntonLessem(temperature, par = c(Tmin, 14, 27, Tmax
temperature[which.max(modifiedThorntonLessem(temperature, par = c(Tmin, 14, 27, Tmax)))]
lines(temperature, modifiedLehman(temperature, par = c(Tmin, 21.3, Tmax)), col = 'red')
StrabakaGnauck = function(temperature, par){
# in Lefrançois, C., et G. Claireaux (2003). Influence of Ambient Oxygenation and Temperature on Metabolic Scope and
# Scope for Heart Rate in the Common Sole Solea Solea. Marine Ecology Progress Series 259 : 273‑84.
Topt = par[1]
Tmax = par[2]
omega = par[3]
rate = ((Tmax - temperature) / (Tmax - Topt))^omega * exp(-omega * (Tmax - temperature) / (Tmax - Topt) )
return(rate)
}
plot(temperature, StrabakaGnauck(temperature, par = c(19.74, 30.48, 1.87)))
\ No newline at end of file
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