Mercurial > repos > sblanck > mpagenomics_normalize

annotate segcall.R @ 6:7dc6ce39fb89 default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
add selection tool
author blanck
date Wed, 29 Apr 2015 10:08:52 +0200
parents b7f3854e08f8
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
5
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
1 args<-commandArgs(TRUE)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
2
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
3 chrom=args[1]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
4 dataset=args[2]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
5 output=args[3]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
6 tmp_dir=args[4]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
7 nbcall=as.numeric(args[5])
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
8 input=args[6]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
9 outputfigures=type.convert(args[7])
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
10 snp=type.convert(args[8])
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
11 tumorcsv=args[9]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
12 cellularity=as.numeric(args[10])
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
13 user=args[11]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
14 method=args[12]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
15
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
16 library(MPAgenomics)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
17 workdir=file.path(tmp_dir, "mpagenomics",user)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
18 setwd(workdir)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
19
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
20 if (grepl("all",tolower(chrom)) | chrom=="None") {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
21 chrom_vec=c(1:25)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
22 } else {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
23 chrom_tmp <- strsplit(chrom,",")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
24 chrom_vecstring <-unlist(chrom_tmp)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
25 chrom_vec <- as.numeric(chrom_vecstring)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
26 }
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
27
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
28 input_tmp <- strsplit(input,",")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
29 input_tmp_vecstring <-unlist(input_tmp)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
30
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
31
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
32 input_vecstring = sub("^([^.]*).*", "\\1", input_tmp_vecstring)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
33
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
34 if (dataset == input) {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
35 if (tumorcsv== "none")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
36 {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
37 segcall=cnSegCallingProcess(dataset,chromosome=chrom_vec, nclass=nbcall, savePlot=outputfigures,onlySNP=snp, cellularity=cellularity, method=method)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
38 } else {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
39 segcall=cnSegCallingProcess(dataset,chromosome=chrom_vec, normalTumorArray=tumorcsv, nclass=nbcall, savePlot=outputfigures,onlySNP=snp, cellularity=cellularity, method=method)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
40 }
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
41 } else {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
42 if (tumorcsv== "none")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
43 {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
44 segcall=cnSegCallingProcess(dataset,chromosome=chrom_vec, listOfFiles=input_vecstring, nclass=nbcall, savePlot=outputfigures, onlySNP=snp, cellularity=cellularity, method=method)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
45 } else {
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
46 segcall=cnSegCallingProcess(dataset,chromosome=chrom_vec, normalTumorArray=tumorcsv, listOfFiles=input_vecstring, nclass=nbcall, savePlot=outputfigures, onlySNP=snp, cellularity=cellularity, method=method)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
47 }
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
48 }
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
49
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
50 sink(output)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
51 print(format(segcall))
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
52 sink()
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
53 #write.table(format(segcall),output,row.names = FALSE, quote=FALSE, sep = "\t")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
54 #write.fwf(segcall,output,rownames = FALSE, quote=FALSE, sep = "\t")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
55 quit()