library(lubridate)
library(data.table)
library(ggplot2)

market <- as.data.frame(fread("~/Project/FHT_Salford/Results/NRR_INVOICEDATE_NRR1ONLY0322.csv"))

market$INVOICEDATEMKEY <- ymd(market$INVOICEDATEMKEY)
market$MAIL_DT <- ymd_hms(market$MAIL_DT)

market$DAY_TO_ACT <- as.numeric(difftime(market$INVOICEDATEMKEY, market$MAIL_DT, units = 'days'))


ggplot(market[market$PRIMARYRESPONSEMODEL %in% c("TU16", "ZENTU16") & market$MAIL_DT %in% market$MAIL_DT[4],],aes(x=DAY_TO_ACT,color=PRIMARYRESPONSEMODEL)) +
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("TU16") & market$MAIL_DT %in% market$MAIL_DT[4],],aes(y=cumsum(..count..)),geom="step")+
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("ZENTU16") & market$MAIL_DT %in% market$MAIL_DT[4],],aes(y=cumsum(..count..)),geom="step")+
    ggtitle("Net responders from TU (Jan campaign)")+
    theme(plot.title = element_text(hjust = 0.5,face = "bold",color="green",size=25))

ggplot(market[market$PRIMARYRESPONSEMODEL %in% c("TU16", "ZENTU16") & market$MAIL_DT %in% market$MAIL_DT[1],],aes(x=DAY_TO_ACT,color=PRIMARYRESPONSEMODEL)) +
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("TU16") & market$MAIL_DT %in% market$MAIL_DT[1],],aes(y=cumsum(..count..)),geom="step")+
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("ZENTU16") & market$MAIL_DT %in% market$MAIL_DT[1],],aes(y=cumsum(..count..)),geom="step")+
    ggtitle("Net responders from TU (Feb campaign)")+
    theme(plot.title = element_text(hjust = 0.5,face = "bold",color="green",size=25))


ggplot(market[market$PRIMARYRESPONSEMODEL %in% c("EQ16V2", "ZENEQ16") & market$MAIL_DT %in% market$MAIL_DT[4],],aes(x=DAY_TO_ACT,color=PRIMARYRESPONSEMODEL)) +
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("EQ16V2") & market$MAIL_DT %in% market$MAIL_DT[4],],aes(y=cumsum(..count..)),geom="step")+
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("ZENEQ16") & market$MAIL_DT %in% market$MAIL_DT[4],],aes(y=cumsum(..count..)),geom="step")+
    ggtitle("Net responders from EQ (Jan campaign)")+
    theme(plot.title = element_text(hjust = 0.5,face = "bold",color="green",size=25))

ggplot(market[market$PRIMARYRESPONSEMODEL %in% c("EQ16V2", "ZENEQ16") & market$MAIL_DT %in% market$MAIL_DT[1],],aes(x=DAY_TO_ACT,color=PRIMARYRESPONSEMODEL)) +
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("EQ16V2") & market$MAIL_DT %in% market$MAIL_DT[1],],aes(y=cumsum(..count..)),geom="step")+
    stat_bin(data=market[market$PRIMARYRESPONSEMODEL %in% c("ZENEQ16") & market$MAIL_DT %in% market$MAIL_DT[1],],aes(y=cumsum(..count..)),geom="step")+
    ggtitle("Net responders from EQ (Feb campaign)")+
    theme(plot.title = element_text(hjust = 0.5,face = "bold",color="green",size=25))