annotate markersSelection.R @ 6:7dc6ce39fb89 default tip

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 input=args[1]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
4 response=args[2]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
5 tmp_dir=args[3]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
6 nbFolds=as.numeric(args[4])
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
7 loss=args[5]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
8 output=args[6]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
9
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
10 library(MPAgenomics)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
11 workdir=file.path(tmp_dir, "mpagenomics")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
12 setwd(workdir)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
13
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
14 CN=read.table(input,header=TRUE,check.names=FALSE)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
15 drops=c("chromosome","position","probeName")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
16 CNsignal=CN[,!(names(CN)%in% drops)]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
17 samples=names(CNsignal)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
18 CNsignalMatrix=t(data.matrix(CNsignal))
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
19 resp=read.table(response,header=TRUE,sep=",")
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
20 listOfFile=resp[[1]]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
21 responseValue=resp[[2]]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
22 index = match(listOfFile,rownames(CNsignalMatrix))
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
23 responseValueOrder=responseValue[index]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
24
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
25 result=variableSelection(CNsignalMatrix,responseValueOrder,nbFolds=nbFolds,loss=loss,plot=TRUE)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
26
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
27 CNsignalResult=CN[result$markers.index,(names(CN)%in% drops)]
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
28
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
29 CNsignalResult["coefficient"]=result$coefficient
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
30 CNsignalResult["index"]=result$markers.index
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
31
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
32 sink(output)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
33 print(format(CNsignalResult),row.names=FALSE)
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
34 sink()
b7f3854e08f8 Adding all tools
blanck
parents:
diff changeset
35 #write.table(CNsignalResult,output,row.names = FALSE, quote=FALSE, sep = "\t")