Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- library("pryr")
- myCumulativePlot <- function (data1,data2){
- for(x in 1:10){
- a%<a-% {
- plot(ecdf(data1[x,]),col="red",lwd = 2,verticals = TRUE, do.points = FALSE)
- legend("topleft",c("Enriched","Depleted"),lwd=c(2.5,2.5),col=c("red","blue"))
- lines(ecdf(data2[x,]),col="blue",xaxt='n',yaxt='n', ann=FALSE,lwd = 2,verticals = TRUE, do.points = FALSE)
- }
- myData[[x]]=a
- }
- return(myData)
- }
Add Comment
Please, Sign In to add comment