Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/env Rscript
- # clargs
- args = commandArgs(trailingOnly=TRUE)
- # microarr messung
- microarr <- read.table(args[1], sep="\t", header=T)
- # phosphoroproteomics
- phospho <- read.table(args[2], sep="\t", header=T)
- # merge
- df <-merge(microarr,phospho, by.x="YORF", by.y="ORF", all=F)
- # viewer zeigt es gibt nur werte für manche... sonst NA, deswegen all=F
- # plots
- outpath <- paste(args[3], "/plot_collection.pdf", sep="")
- pdf(outpath)
- plot(df$X0.min, df$fc)
- plot(df$X15.min, df$fc)
- plot(df$X30.min, df$fc)
- plot(df$X45.min, df$fc)
- plot(df$X60.min, df$fc)
- plot(df$X90.min, df$fc)
- plot(df$X120.min, df$fc)
- dev.off()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement