plot_data_and_model.R

#Script to plot data and model

growth_data <- read.csv("experiment2.csv")

logistic_fun <- function(t) {
  
  N <- (N0*K*exp(r*t))/(K-N0+N0*exp(r*t))
  
  return(N)
  
}

#exp(7.53849)
N0 <- 1879 #
  
r <- 0.031 #
  
K <- 999979341 #

ggplot(aes(t,N), data = growth_data) +
  
  geom_function(fun=logistic_fun, colour="red") +
  
  geom_point() #+

  #scale_y_continuous(trans='log10')