Mercurial > repos > davidvanzessen > phenotype_gene_relations
changeset 22:04d8e7add1ba draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Tue, 10 Nov 2015 10:04:46 -0500 |
| parents | 5806aed6ddec |
| children | 0fa1b346ab52 |
| files | phenotype_gene_relations_plot.r |
| diffstat | 1 files changed, 3 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/phenotype_gene_relations_plot.r Tue Nov 10 09:40:09 2015 -0500 +++ b/phenotype_gene_relations_plot.r Tue Nov 10 10:04:46 2015 -0500 @@ -59,7 +59,7 @@ dat.sub2$disease.gene.order = dat.sub2$ratio for(orderby in strsplit(orderby.list, ",")[[1]]){ - filter = dat.sub2$GROUP.CODE == as.numeric(orderby) + 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) @@ -68,6 +68,8 @@ dat.sub2[dat.sub2$disease.gene.order < max(dat.sub2$ratio),"disease.gene.order"] = 0 +print(paste("min ratio:", min(dat.sub2$ratio))) + write.table(dat.sub2[,c("disease.gene", "GROUP.CODE", "GROUP.NAME", "ratio", "disease.gene.order")], "debug.txt", quote=F, sep="\t", row.names=F, col.names=T) #dat.sub2.order = dat.sub2[dat.sub2$GROUP.CODE == orderby,c("disease.gene", "ratio")]