Advertisement
zsewa

Untitled

May 27th, 2019
95
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.62 KB | None | 0 0
  1. #!/usr/bin/env Rscript
  2. # clargs
  3. args = commandArgs(trailingOnly=TRUE)
  4.  
  5. # microarr messung
  6. microarr <- read.table(args[1], sep="\t", header=T)
  7. # phosphoroproteomics
  8. phospho <- read.table(args[2], sep="\t", header=T)
  9.  
  10. # merge
  11. df <-merge(microarr,phospho, by.x="YORF", by.y="ORF", all=F)
  12. # viewer zeigt es gibt nur werte für manche... sonst NA, deswegen all=F
  13.  
  14. # plots
  15. outpath <- paste(args[3], "/plot_collection.pdf", sep="")
  16. pdf(outpath)
  17. plot(df$X0.min, df$fc)
  18. plot(df$X15.min, df$fc)
  19. plot(df$X30.min, df$fc)
  20. plot(df$X45.min, df$fc)
  21. plot(df$X60.min, df$fc)
  22. plot(df$X90.min, df$fc)
  23. plot(df$X120.min, df$fc)
  24. dev.off()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement