source: trunk/brs2010p35/src/jag/generate_single_qq.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: 676 bytes
Line 
1
2generate_qq <- function(pvalues, group, filename, adjusted)
3{
4        pdf.title <- paste(filename)
5        pdf(file=pdf.title)
6       
7        plot_qq(pvalues,group)
8        dev.off()
9}
10
11#read all arguments given by python
12#DO NOT REMOVE!
13args=(commandArgs(TRUE))
14if(length(args)==0){
15    print("No arguments supplied.")
16}else{
17    for(i in 1:length(args)){
18         eval(parse(text=args[[i]]))
19    }
20}
21#loading R file wit all functions
22source(paste(r_path,"global_functions.R",sep=""))
23
24input <- read.delim(assocfile, sep = '', header = TRUE, colClasses = "character")
25
26if(adjusted == 'True')
27{
28        pvalues <- input$GC
29}
30
31if(adjusted == 'False')
32{
33        pvalues <- input$P
34}
35
36
37generate_qq(pvalues, header, filename)
Note: See TracBrowser for help on using the repository browser.