Mercurial > repos > drosofff > msp_sr_readmap_and_size_histograms

comparison readmap.xml @ 10:71a46afb9ce7 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
author mvdbeek
date Tue, 02 Feb 2016 12:54:00 -0500
parents 44a0b0fa52a9
children edd57052f0bb
comparison
equal deleted inserted replaced
9:44a0b0fa52a9 10:71a46afb9ce7
95 library(gridExtra) 95 library(gridExtra)
96 96
97 ## data frames implementation 97 ## data frames implementation
98 98
99 rm=read.delim("${readmap_dataframe}", header=T, row.names=NULL) 99 rm=read.delim("${readmap_dataframe}", header=T, row.names=NULL)
100 n_samples=length(unique(rm$sample)) 100 n_samples=length(unique(rm\$sample))
101 genes=unique(levels(rm$gene)) 101 genes=unique(levels(rm$gene))
102 per_gene_readmap=lapply(genes, function(x) subset(rm, gene==x)) ####### ? 102 per_gene_readmap=lapply(genes, function(x) subset(rm, gene==x)) ####### ?
103 n_genes=length(per_gene_readmap) 103 n_genes=length(per_gene_readmap)
104 104
105 size=read.delim("${size_distribution_dataframe}", header=T, row.names=NULL) 105 size=read.delim("${size_distribution_dataframe}", header=T, row.names=NULL)