smagexp: MetaRNASeq.R annotate

annotate MetaRNASeq.R @ 3:63355796cdc8 draft

planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit ecb5d3451f1623e0cc649ac5197e66995ce19eef

author	sblanck
date	Thu, 22 Feb 2018 10:10:34 -0500
parents	1024245abc70
children	58052f8bc987

rev	line source
0 1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	1 #!/usr/bin/env Rscript
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	2 # setup R error handling to go to stderr
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	3 options( show.error.messages=F, error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	4
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	5 # we need that to not crash galaxy with an UTF8 error on German LC settings.
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	6 loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	7
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	8 library("optparse")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	9
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	10 ##### Read options
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	11 option_list=list(
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	12 make_option("--input",type="character",default="NULL",help="list of rdata objects containing eset objects"),
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	13 make_option("--result",type="character",default=NULL,help="text file containing result of the meta-analysis"),
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	14 make_option("--htmloutput",type="character",default=NULL,help="Output html report"),
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	15 make_option("--htmloutputpath",type="character",default="NULL",help="Path of output html report"),
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	16 make_option("--htmltemplate",type="character",default=NULL,help="html template)")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	17 );
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	18
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	19 opt_parser = OptionParser(option_list=option_list);
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	20 opt = parse_args(opt_parser);
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	21
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	22 if(is.null(opt$input)){
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	23 print_help(opt_parser)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	24 stop("input required.", call.=FALSE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	25 }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	26
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	27 #loading libraries
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	28
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	29 suppressPackageStartupMessages(require(affy))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	30 suppressPackageStartupMessages(require(annaffy))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	31 suppressPackageStartupMessages(require(VennDiagram))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	32 suppressPackageStartupMessages(require(GEOquery))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	33
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	34 listInput <- trimws( unlist( strsplit(trimws(opt$input), ",") ) )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	35
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	36 listfiles=vector()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	37 listfilenames=vector()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	38
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	39 for (i in 1:length(listInput))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	40 {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	41 inputFileInfo <- unlist( strsplit( listInput[i], ';' ) )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	42 listfiles=c(listfiles,inputFileInfo[1])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	43 listfilenames=c(listfilenames,inputFileInfo[2])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	44 }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	45
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	46 cargs <- commandArgs()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	47 cargs <- cargs[(which(cargs == "--args")+1):length(cargs)]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	48 nbargs=length(cargs)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	49 listfiles=vector()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	50 listfilenames=vector()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	51 for (i in seq(1,nbargs-6,2)) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	52 listfiles=c(listfiles,cargs[[i]])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	53 listfilenames=c(listfilenames,cargs[[i+1]])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	54 }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	55 #mod<-cargs[[length(cargs) - 6]]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	56 outputfile <- opt$result
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	57 result.html = opt$htmloutput
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	58 html.files.path=opt$htmloutputpath
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	59 result.template=opt$htmltemplate
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	60
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	61 alpha=0.05
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	62
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	63 #print(comparison)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	64
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	65 listData=lapply(listfiles,read.table)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	66 orderData=lapply(listData, function(x) x[order(x[1]), ])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	67 rawpval=lapply(orderData,function(x) x[6])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	68 rawpval=lapply(rawpval, function(x) as.numeric(unlist(x)))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	69
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	70 DE=list()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	71 DE=lapply(orderData, function(x) ifelse(x[7]<=0.05,1,0))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	72
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	73 FC=list()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	74 FC=lapply(orderData, function(x) x[3])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	75
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	76 DE=as.data.frame(DE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	77 colnames(DE)=listfilenames
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	78 FC=as.data.frame(FC)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	79 colnames(FC)=listfilenames
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	80 # the comparison must only have two values and the conds must
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	81 # be a vector from those values, at least one of each.
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	82
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	83 #if (length(comparison) != 2) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	84 # stop("Comparison type must be a tuple: ", cargs[length(cargs) - 8])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	85 #}
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	86
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	87 sink("/dev/null")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	88 dir.create(html.files.path, recursive=TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	89 #library(DESeq)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	90 #library(HTSFilter)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	91
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	92 #DE=list()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	93 #FC=list()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	94 #i=1
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	95
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	96 # Open the html output file
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	97 #file.conn = file(diag.html, open="w")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	98
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	99 #writeLines( c("<html><body>"), file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	100
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	101 # Perform deseq analysis on each study
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	102 #for(i in 1:length(listfiles))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	103 #{
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	104 # f=listfiles[i]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	105 # fname=listfilenames[i]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	106 # study_name=unlist(strsplit(fname,"[.]"))[1]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	107 # print(paste0("study.name ",study_name))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	108 # d <- read.table(f, sep=" ", header=TRUE, row.names=1)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	109 # conds<-sapply(strsplit(colnames(d),"[.]"),FUN=function(x) x[1])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	110 # if (length(unique(conds)) != 2) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	111 # warning(as.data.frame(strsplit(colnames(d),"[.]")))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	112 # stop("You can only have two columns types: ", paste(conds,collapse=" "))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	113 # }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	114 # if (!identical(sort(comparison), sort(unique(conds)))) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	115 # stop("Column types must use the two names from Comparison type, and vice versa. Must have at least one of each in the Column types.\nColumn types: ", cargs[2], "\n", "Comparison type: ", cargs[3])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	116 # }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	117 # if (length(d) != length(conds)) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	118 # stop("Number of total sample columns in counts file must correspond to the columns types field. E.g. if column types is 'kidney,kidney,liver,liver' then number of sample columns in counts file must be 4 as well.")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	119 # }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	120 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	121 # cds <- newCountDataSet(d, conds)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	122 # cds <- estimateSizeFactors(cds)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	123 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	124 # cdsBlind <- estimateDispersions( cds, method="blind" )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	125 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	126 # if (length(conds) != 2) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	127 # cds <- estimateDispersions( cds )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	128 # norep = FALSE
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	129 # }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	130 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	131 # if (length(conds) == 2) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	132 # cds <- estimateDispersions( cds, method=method, sharingMode=mod, fitType="parametric" )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	133 # norep = TRUE
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	134 # }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	135 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	136 # filter<-HTSFilter(cds, plot=FALSE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	137 # cds.filter<-filter$filteredData
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	138 # on.index<-which(filter$on==1)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	139 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	140 # res<-as.data.frame(matrix(NA,nrow=nrow(cds),ncol=ncol(cds)))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	141 # nbT <- nbinomTest(cds.filter, comparison[1], comparison[2])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	142 # colnames(res)<-colnames(nbT)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	143 # res[on.index,]<-nbT
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	144 # #write.table(res[order(res$padj), ], file=outputfile, quote=FALSE, row.names=FALSE, sep="\t")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	145 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	146 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	147 # temp.pval.plot = file.path( html.files.path, paste("PvalHist",i,".png",sep=""))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	148 # png( temp.pval.plot, width=500, height=500 )
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	149 # hist(res$pval, breaks=100, col="skyblue", border="slateblue", main="")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	150 # dev.off()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	151 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	152 # writeLines( c("<h2>P-value histogram for ",study_name,"</h2>"), file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	153 # writeLines( c("<img src='PvalHist",i,".png'><br/><br/>"), file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	154 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	155 # #on enregistre la p-value
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	156 # rawpval[[study_name]]<-res$pval
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	157 # DE[[study_name]]<-ifelse(res$padj<=alpha,1,0)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	158 # FC[[study_name]]<-res$log2FoldChange
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	159 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	160 # i=i+1
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	161 #}
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	162
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	163
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	164 # combinations
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	165 library(metaRNASeq)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	166 fishcomb<-fishercomb(rawpval, BHth=alpha)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	167 warning(length(rawpval))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	168 invnormcomb<-invnorm(rawpval, nrep=c(8,8), BHth=alpha)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	169 #DE[["fishercomb"]]<-ifelse(fishcomb$adjpval<=alpha,1,0)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	170 #DE[["invnormcomb"]]<-ifelse(invnormcomb$adjpval<=alpha,1,0)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	171
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	172 signsFC<-mapply(FC,FUN=function(x) sign(x))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	173 sumsigns<-apply(signsFC,1,sum)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	174 commonsgnFC<-ifelse(abs(sumsigns)==dim(signsFC)[2],sign(sumsigns),0)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	175
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	176 DEresults <- data.frame(DE=DE,"DE.fishercomb"=ifelse(fishcomb$adjpval<=alpha,1,0),"DE.invnorm"=ifelse(invnormcomb$adjpval<=alpha,1,0))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	177
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	178 unionDE <- unique(c(fishcomb$DEindices,invnormcomb$DEindices))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	179 FC.selecDE <- data.frame(DEresults[unionDE,],FC[unionDE,],signFC=commonsgnFC[unionDE])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	180 keepDE <- FC.selecDE[which(abs(FC.selecDE$signFC)==1),]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	181
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	182 fishcomb_de <- rownames(keepDE)[which(keepDE[,"DE.fishercomb"]==1)]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	183 invnorm_de <- rownames(keepDE)[which(keepDE[,"DE.invnorm"]==1)]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	184 indstudy_de = list()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	185 for (i in 1:length(listfiles)) {
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	186 currentIndstudy_de = rownames(keepDE)[which(keepDE[,i]==1)]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	187 indstudy_de[[listfilenames[i]]]=currentIndstudy_de
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	188 }
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	189
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	190 IDDIRRfishcomb=IDD.IRR(fishcomb_de,indstudy_de)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	191 IDDIRRinvnorm=IDD.IRR(invnorm_de,indstudy_de)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	192
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	193 #conflits<-data.frame(ID=listData[[1]][rownames(DEresults),1],Fishercomb=DEresults[["DE.fishercomb"]],Invnormcomb=DEresults[["DE.invnorm"]],sign=commonsgnFC)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	194 conflits<-data.frame(ID=listData[[1]][rownames(DEresults),1],DE=DEresults,FC=FC,signFC=commonsgnFC)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	195 #write DE outputfile
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	196 write.table(conflits, outputfile,sep="\t",,row.names=FALSE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	197 library(VennDiagram)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	198 DE_num=apply(DEresults, 2, FUN=function(x) which(x==1))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	199 venn.plot<-venn.diagram(x=as.list(DE_num),filename=NULL, col="black", fill=1:length(DE_num)+1,alpha=0.6)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	200 temp.venn.plot = file.path( html.files.path, paste("venn.png"))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	201 png(temp.venn.plot,width=500,height=500)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	202 grid.draw(venn.plot)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	203 dev.off()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	204
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	205 library(jsonlite)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	206 matrixConflits=as.matrix(conflits)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	207 datajson=toJSON(matrixConflits,pretty = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	208 summaryFishcombjson=toJSON(as.matrix(t(IDDIRRfishcomb)),pretty = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	209 summaryinvnormjson=toJSON(as.matrix(t(IDDIRRinvnorm)),pretty = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	210
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	211
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	212 #vennsplit=strsplit(result.venn,split="/")[[1]]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	213 #venn=paste0("./",vennsplit[length(vennsplit)])
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	214
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	215
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	216 vennFilename="venn.png"
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	217 vennFile=file.path(html.files.path,vennFilename)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	218 htmlfile=readChar(result.template, file.info(result.template)$size)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	219 htmlfile=gsub(x=htmlfile,pattern = "###DATAJSON###",replacement = datajson, fixed = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	220 htmlfile=gsub(x=htmlfile,pattern = "###FISHSUMMARYJSON###",replacement = summaryFishcombjson, fixed = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	221 htmlfile=gsub(x=htmlfile,pattern = "###INVSUMMARYJSON###",replacement = summaryinvnormjson, fixed = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	222 htmlfile=gsub(x=htmlfile,pattern = "###VENN###",replacement = vennFilename, fixed = TRUE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	223 write(htmlfile,result.html)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	224
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	225 #library(VennDiagram)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	226 #flog.threshold(ERROR)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	227 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	228 ##venn.plot<-venn.diagram(x = c(res[c(1:(length(res)-3))],meta=list(res$Meta)),filename = v, col = "black", fill = c(1:(length(res)-2)), margin=0.05, alpha = 0.6,imagetype = "png")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	229 #dir.create(result.path, showWarnings = TRUE, recursive = FALSE)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	230 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	231 #showVenn<-function(liste,file)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	232 #{
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	233 # venn.plot<-venn.diagram(x = liste,
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	234 # filename = vennFilename, col = "black",
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	235 # fill = 1:length(liste)+1,
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	236 # margin=0.05, alpha = 0.6,imagetype = "png")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	237 ## png(file);
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	238 ## grid.draw(venn.plot);
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	239 ## dev.off();
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	240 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	241 #}
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	242 #
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	243 #l=list()
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	244 #for(i in 1:length(esets))
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	245 #{
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	246 # l[[paste("study",i,sep="")]]<-res[[i]]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	247 #}
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	248 #l[["Meta"]]=res[[length(res)-1]]
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	249 #showVenn(l,vennFile)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	250 #file.copy(vennFilename,result.path)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	251
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	252
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	253 #writeLines( c("<h2>Venn Plot</h2>"), file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	254 #writeLines( c("<img src='venn.png'><br/><br/>"), file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	255 #writeLines( c("</body></html>"), file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	256 #close(file.conn)
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	257 #print("passe6")
1024245abc70 planemo upload for repository https://github.com/sblanck/smagexp/tree/master/smagexp_tools commit 5974f806f344dbcc384b931492d7f023bfbbe03b sblanck parents: diff changeset	258 #sink(NULL)

Mercurial > repos > sblanck > smagexp

annotate MetaRNASeq.R @ 3:63355796cdc8 draft