source: trunk/brs2010p35/src/jag/multiple_qq_plot_for_sets.R @ 2

Last change on this file since 2 was 2, checked in by tim.te.beek@…, 8 years ago

Add revision 338 snapshot of previous SVN instance, omitting PLINK binaries

File size: 504 bytes
Line 
1
2#read all arguments given by python
3#DO NOT REMOVE!
4args=(commandArgs(TRUE))
5if(length(args)==0){
6    print("No arguments supplied.")
7}else{
8    for(i in 1:length(args)){
9         eval(parse(text=args[[i]]))
10    }
11}
12#loading R file wit all functions
13source(paste(r_path,"global_functions.R",sep=""))
14
15
16raw_p_table<-read.table(p_values_file)
17pdf(file=image_filename)
18for (set_name in levels(raw_p_table[,1])){
19  p_per_set=raw_p_table[raw_p_table[,1]==set_name,2]
20  plot_qq(p_per_set,set_name)
21}
22dev.off()
Note: See TracBrowser for help on using the repository browser.