# Load Packages ----------------------------------------------------------------
library(tidyverse)
library(plyr)
library(MASS)

# Data importation -------------------------------------------------------------
data <- read_csv("ltdata.csv")
data
str(data)
summary(data)

# Data recording --------------------------------------------------------

# Recording combination, treatment, replicate, day and death variables as factors --------

data <- read_csv(
  "ltdata.csv",
  col_types = cols(
    combination   = col_factor(),
    treatment       = col_factor(),
    replicate = col_factor(),
    day     = col_factor(),
    death = col_factor()
  )
)
data


# Let's calculate LT50 values according treatment and replicate ------------------------------------------------

# Calculer les valeurs predites de LT50 pour chaque traitement en fonction des repliques

LT50_val <- ddply(data, .(treatment, replicate), summarize, 
                  LT50=as.numeric(dose.p(glm(cbind(cumul_alive,cumul_death)~day1,binomial),p=0.5))) #p=0.5 pour LT50
View(LT50_val)

# Let's get means LT50 for each level of treatment ------------------------
LT50_Moy <- ddply(LT50_val, .(treatment), summarize,    
                  LT50_moy=mean(LT50),
                  se=sd(LT50)/sqrt(length(LT50)),
                  Nb_replique=length(LT50))
View(LT50_Moy)