#Survival analysis - example.
#Control:
#======================
a0 <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/guinea_control.txt", header=TRUE)

q0 <- seq(1,nrow(a0),1)
qq0 <- q0/nrow(a0)
a00 <- as.data.frame(cbind(a0$x, qq0))
names(a00) <- c("x", "ecdf0")

b0 <- which(duplicated(a0$x))

b00 <- a00[-(b0-1),]
#=====================

#Dose 1:
#======================
a1 <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/guinea1.txt", header=TRUE)

q1 <- seq(1,nrow(a1),1)
qq1 <- q1/nrow(a1)
a11 <- as.data.frame(cbind(a1$x1, qq1))
names(a11) <- c("x1", "ecdf1")

b1 <- which(duplicated(a1$x1))

b11 <- a11[-(b1-1),]
#=====================

#Dose 2:
#======================
a2 <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/guinea2.txt", header=TRUE)

q2 <- seq(1,nrow(a2),1)
qq2 <- q2/nrow(a2)
a22 <- as.data.frame(cbind(a2$x2, qq2))
names(a22) <- c("x2", "ecdf2")

b2 <- which(duplicated(a2$x2))

b22 <- a22[-(b2-1),]
#=====================

#Dose 3:
#======================
a3 <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/guinea3.txt", header=TRUE)

q3 <- seq(1,nrow(a3),1)
qq3 <- q3/nrow(a3)
a33 <- as.data.frame(cbind(a3$x3, qq3))
names(a33) <- c("x3", "ecdf3")

a333 <- a33[-c(15,17,18,27,39,60),]
#=====================


#Dose 4:
#======================
a4 <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/guinea4.txt", header=TRUE)

q4 <- seq(1,nrow(a4),1)
qq4 <- q4/nrow(a4)
a44 <- as.data.frame(cbind(a4$x4, qq4))
names(a44) <- c("x4", "ecdf4")

a444 <- a44[-c(4,13,15,16,19,24,26,29,31,34,35),]

#========================


Dose 5:
======================
a5 <- read.table("http://www.stat.ucla.edu/~nchristo/statistics13/guinea5.txt", header=TRUE)

q5 <- seq(1,nrow(a5),1)
qq5 <- q5/nrow(a5)
a55 <- as.data.frame(cbind(a5$x5, qq5))
names(a55) <- c("x5", "ecdf5")

a555 <- a55[-c(4,6,10,17,22,25,26,27,32,36,41,60,62,66,70),]

========================


plot(b00$x, 1-b00$ecdf0, type="l", xlab="Days", ylab="Proportion of live animals")
points(b11$x1, 1-b11$ecdf1, type="l",lty=2)
points(b22$x2, 1-b22$ecdf2, type="l",lty=3)


points(a333$x3, 1-a333$ecdf3, type="l",lty=4)
points(a444$x4, 1-a444$ecdf4, type="l",lty=5)
points(a555$x5, 1-a555$ecdf5, type="l",lty=6)


legend('topright', lty=1:6, c("control", "I", "II", "III", "IV", "V"))