Guest User

Untitled

a guest
Oct 18th, 2017
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.44 KB | None | 0 0
  1. library("pryr")
  2. myCumulativePlot <- function (data1,data2){
  3. for(x in 1:10){
  4. a%<a-% {
  5. plot(ecdf(data1[x,]),col="red",lwd = 2,verticals = TRUE, do.points = FALSE)
  6. legend("topleft",c("Enriched","Depleted"),lwd=c(2.5,2.5),col=c("red","blue"))
  7. lines(ecdf(data2[x,]),col="blue",xaxt='n',yaxt='n', ann=FALSE,lwd = 2,verticals = TRUE, do.points = FALSE)
  8. }
  9. myData[[x]]=a
  10. }
  11. return(myData)
  12. }
Add Comment
Please, Sign In to add comment