Please follow the instructions at this link:

https://vectorbyteorg.github.io/bayesTPC-tutorials/intro.html

You can copy-paste the scripts below to run them.

In [ ]:

library(nimble)
library(bayesTPC)

Sys.setenv(PATH = paste("/opt/conda/envs/RGIS/bin",
                        Sys.getenv("PATH"), sep = ":"))

In [ ]:

# Load example data
data <- read.csv("data/AeaegyptiTraitData.csv")
data_list <- list(Trait = data$trait, Temp = data$Temp)

# Quick plot
plot(data_list$Trait ~ data_list$Temp, ylab = "Development Rate", xlab = "Temperature (°C)", pch = 16)

In [ ]:

set.seed(1)
fit <- b_TPC(data = data_list, model = "briere")

In [ ]:

print(fit)

In [ ]:

plot(fit, temp_interval = seq(0, 35, length.out = 1000))