Mercurial > repos > davidvanzessen > phenotype_gene_relations
changeset 24:14afe3c8a3d5 draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Mon, 16 Nov 2015 07:06:56 -0500 |
| parents | 0fa1b346ab52 |
| children | bce9fdf89e48 |
| files | phenotype_gene_relations_plot.r |
| diffstat | 1 files changed, 5 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/phenotype_gene_relations_plot.r Wed Nov 11 08:36:42 2015 -0500 +++ b/phenotype_gene_relations_plot.r Mon Nov 16 07:06:56 2015 -0500 @@ -11,6 +11,7 @@ omt5=strsplit(args[6], ",")[[1]] inheritance=strsplit(args[7], ",")[[1]] orderby.list=as.character(args[8]) +orderby = ifelse(nchar(orderby) == 0, "1", orderby) outdir=args[9] dat = read.table(input, header=T, sep="\t", comment.char="#", quote = "") @@ -62,11 +63,13 @@ filter = dat.sub2$GROUP.CODE == as.numeric(orderby) & dat.sub2$ratio != 0 print(paste("orderby:",orderby, "+1:", sum(filter))) if(sum(filter) > 0){ - dat.sub2[filter,"disease.gene.order"] = dat.sub2[filter,"disease.gene.order"] + max(dat.sub2$ratio) + #dat.sub2[filter,"disease.gene.order"] = dat.sub2[filter,"disease.gene.order"] + max(dat.sub2$ratio) + dat.sub2[filter,"disease.gene.order"] = dat.sub2[filter,"disease.gene.order"] + 1 } } -dat.sub2[dat.sub2$disease.gene.order < max(dat.sub2$ratio),"disease.gene.order"] = 0 +#dat.sub2[dat.sub2$disease.gene.order < max(dat.sub2$ratio),"disease.gene.order"] = 0 +dat.sub2[dat.sub2$disease.gene.order < 1,"disease.gene.order"] = 0 print(paste("min ratio:", min(dat.sub2$ratio)))