WT vs. L3 Empirical count CDF (cumulative sum)

d <- read.csv("/tmp/clusters_list.csv", sep="\t")

levels(d[,"source"]) <- c("L3", "WT")

densitybysource <- with(d, split(density, source))

densitycumsum <- lapply(densitybysource, function (x) cumsum(rev(x)))

y.wt <- cumsum(table(densitybysource$WT))

y.l3 <- cumsum(table(densitybysource$L3))

x.wt <- as.numeric(names(y.wt))

x.l3 <- as.numeric(names(y.l3))

wt <- data.frame(x=x.wt, y=y.wt, d="WT")

l3 <- data.frame(x=x.l3, y=y.l3, d="L3")

d2 <- rbind(wt, l3)

pdf("/tmp/WT_L3_cumsum.pdf")

xyplot(y ~ x, groups=d, data=d2, auto.key=T, type="p")

dev.off()

Done