Mercurial > repos > iuc > raceid_inspecttrajectory
changeset 4:86e2358cf273 draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/raceid3 commit 5d7c4cd914295a81d42bf8baaad4981c23a83059"
| author | iuc |
|---|---|
| date | Mon, 12 Aug 2019 13:03:03 -0400 |
| parents | 4164c0da0a5d |
| children | 69018f285aa3 |
| files | macros.xml macros_cheetah.xml macros_cluster.xml macros_inspectclusters.xml macros_inspecttrajectories.xml macros_trajectory.xml raceid_inspecttrajectories.xml scripts/cluster.R scripts/clusterinspect.R test-data/matrix.genelist test-data/matrix.pdf test-data/matrix.rdat test-data/out_diffgene_multiple.pdf test-data/out_traject_inspect_fateid.tabular |
| diffstat | 13 files changed, 604 insertions(+), 1390 deletions(-) [+] |
line wrap: on
line diff
--- a/macros.xml Mon Apr 15 17:53:28 2019 -0400 +++ b/macros.xml Mon Aug 12 13:03:03 2019 -0400 @@ -1,5 +1,4 @@ <macros> - <token name="@GET_VERSION@"><![CDATA[ 2>&1 | head -1 | cut -d' ' -f 2 ]]></token> <token name="@STRING2VECTOR@"> string2numericvector = function(string){ #' Converts a string to a valid vector @@ -27,7 +26,25 @@ } </token> <token name="@VERSION_RACEID@">3</token> - <token name="@VERSION_PACKAGE@">0.2</token> + <token name="@VERSION_WRAPPER@">1</token> + + <macro name="version_command_config" token_prog="temp" token_cheetah="temp2" token_out="2> '$outlog'"> + <version_command><![CDATA[ +Rscript '$__tool_directory__/scripts/@PROG@' 2>&1 | head -1 | cut -d' ' -f 2 + ]]></version_command> + <command detect_errors="exit_code"><![CDATA[ +Rscript '$__tool_directory__/scripts/@PROG@' '$userconf' @OUT@ +]]> + </command> + <configfiles> + <configfile name="userconf" ><![CDATA[ +@STRING2VECTOR@ + +@@CHEETAH@@ +]]></configfile> + </configfiles> + </macro> + <macro name="requirements" > <requirements> <requirement type="package" version="0.1.1">r-raceid</requirement> @@ -69,36 +86,6 @@ <citation type="doi">10.1038/nmeth.4662</citation> </citations> </macro> - <macro name="select_cells" token_sectionname="@SECTIONNAME@" token_sectiontitle="@SECTIONTITLE@" > - <section name="@SECTIONNAME@" title="@SECTIONTITLE@" expanded="true" > - <param name="name_set" type="text" optional="true" label="Name of Set" > - <!-- Aname, Bname --> - <expand macro="sanitize_title" /> - </param> - <conditional name="meth" > - <param name="type" type="select" label="Selection method" > - <option value="cln" selected="true">Cluster Numbers</option> - <option value="regex" >Regular Expression</option> - <option value="manual" >Manual Selection</option> - </param> - <when value="cln" > - <param name="selector" type="text" value="" label="List of clusters" > - <expand macro="sanitize_numeric_vector" /> - </param> - </when> - <when value="regex" > - <param name="selector" type="text" value="" label="Regular Expression to select gene names" > - <expand macro="sanitize_regex" /> - </param> - </when> - <when value="manual" > - <param name="selector" type="text" value="" label="List of genes" > - <expand macro="sanitize_string_vector" /> - </param> - </when> - </conditional> - </section> - </macro> <!-- Sanitizers --> <macro name="sanitize_title"> <sanitizer invalid_char="">
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/macros_cheetah.xml Mon Aug 12 13:03:03 2019 -0400 @@ -0,0 +1,377 @@ +<macros> + <token name="@FILTNORM_CHEETAH@"><![CDATA[ +## Perform do.filter +use.filtnormconf = TRUE + +## Perform do.cluster, do.outlier, do.clustmap, mkgenelist +use.cluster = FALSE + +in.table = read.table( + '${intable}', + stringsAsFactors = F, + na.strings=c("NA", "-", "?", "."), + sep='\t', + header=TRUE, + row.names=1 +) + +## Hidden flag to use test data instead +## see: test-data/use.intestinal + +use.test.data = (names(in.table)[1] == "test") + +sc = NULL +if (use.test.data) { + sc = SCseq(intestinalData) + message("Loading test data from library") +} else { + sc = SCseq(in.table) +} + + +filt = formals(filterdata) +filt.ccc = formals(CCcorrect) +filt.use.ccorrect = FALSE +filt.lbatch.regexes = NULL + +filt.geqone = as.logical( '$filt.hist_geq_one' ) +filt\$mintotal = as.integer( '$filt.mintotal' ) +filt\$minexpr = as.integer( '$filt.minexpr' ) +filt\$minnumber = as.integer( '$filt.minnumber' ) +#if str($filt.use.def) == "no": +filt\$knn = as.integer( '$filt.use.knn' ) +filt\$ccor = as.numeric( '$filt.use.ccor' ) +filt\$bmode = as.character( '$filt.use.bmode' ) + #if $filt.use.LBatch_regexes: +filt.lbatch.regexes = string2textvector( '$filt.use.LBatch_regexes' ) + #end if + #if $filt.use.CGenes: +filt\$CGenes = string2textvector( '$filt.use.CGenes' ) + #end if + #if $filt.use.FGenes: +filt\$FGenes = string2textvector( '$filt.use.FGenes' ) + #end if + #if str($filt.use.ccc.use) == "yes" +filt.use.ccorrect = TRUE + #if $filt.use.ccc.vset: +filt.ccc\$vset = string2textvector( '$filt.use.ccc.vset' ) + #end if + #if $filt.use.ccc.ncomp: +filt.ccc\$nComp = as.integer( '$filt.use.ccc.ncomp' ) + #end if +filt.ccc\$pvalue = as.numeric( '$filt.use.ccc.pvalue' ) +filt.ccc\$quant = as.numeric( '$filt.use.ccc.quant' ) +filt.ccc\$dimR = as.logical( '$filt.use.ccc.dimr' ) +filt.ccc\$mode = as.character( '$filt.use.ccc.mode.value' ) +filt.ccc\$logscale = as.logical( '$filt.use.ccc.logscale' ) + #end if +#end if + +out.pdf = '${outpdf}' +out.rdat = '${outrdat}' +out.table = '${outtable}' + +]]></token> + <token name="@CLUSTER_CHEETAH@"><![CDATA[ + +in.rdat = readRDS('${inputrds}') + +sc = in.rdat + +## Perform do.filter +use.filtnormconf = FALSE + +## Perform do.cluster, do.outlier, do.clustmap, mkgenelist +use.cluster = TRUE + + +clust.compdist = formals(compdist) +clust.clustexp = formals(clustexp) +clust.compdist\$metric = as.character( '$clust.metric' ) +clust.clustexp\$FUNcluster = as.character( '$clust.funcluster' ) + +#if str($clust.use.def) == "no": + +clust.compdist\$FSelect = as.logical( '$clust.use.fselect' ) + #if $clust.use.knn: +clust.compdist\$knn = as.integer( '$clust.use.knn' ) + #end if +clust.clustexp\$sat = as.logical( '$clust.use.sat' ) + #if $clust.use.samp: +clust.clustexp\$samp = as.integer( '$clust.use.samp' ) + #end if + #if $clust.use.cln: +clust.clustexp\$cln = as.integer( '$clust.use.cln' ) +clust.clustexp\$clustnr = as.integer( '$clust.use.clustnr' ) +clust.clustexp\$bootnr = as.integer( '$clust.use.bootnr' ) +##clust.clustexp\$rseed = as.integer( '$clust.use.rseed' ) + #end if +#end if + +outlier.use.randomforest = FALSE +outlier.findoutliers = formals(findoutliers) +outlier.clustheatmap = formals(clustheatmap) +outlier.rfcorrect = formals(rfcorrect) + +outlier.findoutliers\$outminc = as.integer( '$outlier.outminc' ) +outlier.findoutliers\$outlg = as.integer( '$outlier.outlg' ) +outlier.rfcorrect\$final = as.logical( '$outlier.final' ) + +#if str($outlier.use.def) == "no": + #if $outlier.use.nbtree: +outlier.rfcorrect\$nbtree = as.integer( '$outlier.use.nbtree' ) + #end if +outlier.findoutliers\$probthr = as.numeric( '$outlier.use.probthr' ) +outlier.findoutliers\$outdistquant = as.numeric( '$outlier.use.outdistquant' ) +##outlier.rfcorrect\$rfseed = as.integer( '$outlier.use.rfseed' ) +outlier.rfcorrect\$nbfactor = as.integer( '$outlier.use.nbfactor' ) +#end if + +cluster.comptsne = formals(comptsne) +cluster.compfr = formals(compfr) + +cluster.comptsne\$perplexity = as.integer( '$tsne.perplexity' ) +cluster.compfr\$knn = as.integer( '$tsne.knn' ) +#if str($tsne.use.def) == "no": +cluster.comptsne\$initial_cmd = as.logical( '$tsne.use.initial_cmd' ) +cluster.comptsne\$rseed = as.integer( '$tsne.use.rseed_tsne' ) +cluster.compfr\$rseed = as.integer( '$tsne.use.rseed_fr' ) +#end if + +genelist.tablelim = as.integer( '$extra.tablelim' ) +genelist.plotlim = as.integer( '$extra.plotlim' ) +genelist.foldchange = as.integer( '$extra.foldchange' ) +genelist.pvalue = as.numeric( '$extra.pvalue' ) + +out.pdf = '${outpdf}' +out.rdat = '${outrdat}' +out.genelist = '${outgenelist}' +out.assignments = '${outassignments}' +out.log = '${outlog}' + +]]> + </token> + <token name="@INSPECTCLUSTERS_CHEETAH@"><![CDATA[ +in.rdat = readRDS('${inputrds}') + +perform.plotting = FALSE +perform.symbolmap = FALSE +perform.genesofinterest = FALSE +perform.diffgene = FALSE + +#if str($plotgen.do_opt) == "yes": +perform.plotting = TRUE +plotting.cln = string2numericvector( '$plotgen.clusts_plot' ) +#end if + +#if str($plotsym.do_opt) == "yes": +perform.symbolmap = TRUE +plotsym = formals(plotsymbolsmap) +plotsym.use.typeremoveregex = NULL +plotsym.use.typeremoveregex.subselect = NULL +plotsym.use.typeremoveregex = as.character( '$plotsym.types_regex' ) + #if str($plotsym.use.def) == "no": +plotsym\$fr = as.logical( '$plotsym.use.fr' ) + #if $plotsym.use.subset_regex: +plotsym.use.typeremoveregex.subselect = as.character( '$plotsym.use.subset_regex' ) + #end if + #end if +#end if + +#if str($gois.do_opt) == "yes": +perform.genesofinterest = TRUE +plotexp = formals(plotexpmap) +plotmarkg = formals(plotmarkergenes) + +inspect.goi.use.genes = string2textvector( '$gois.inspect_goi_genes' ) +plotexp\$g = inspect.goi.use.genes +plotmarkg\$genes = inspect.goi.use.genes + + #if $gois.inspect_goi_cells: +inspect.goi.use.cells = as.character( '$gois.inspect_goi_cells' ) +plotexp\$cells = inspect.goi.use.cells +plotmarkg\$cells = inspect.goi.use.cells + #end if + #if str($gois.use.def) == "no": +inspect.goi.use.imputed = as.logical( '$gois.use.inspect_goi_imputed' ) +plotexp\$imputed = inspect.goi.use.imputed +plotmarkg\$imputed = inspect.goi.use.imputed + #if $gois.use.plotmarkg_cl +plotmarkg\$cl = string2numericvector( '$gois.use.plotmarkg_cl' ) + #end if + #if $gois.use.plotmarkg_cap +plotmarkg\$cap = as.integer( '$gois.use.plotmarkg_cap' ) + #end if + #if $gois.use.plotmarkg_flo +plotmarkg\$flo = as.integer( '$gois.use.plotmarkg_flo' ) + #end if + #if $gois.use.plotmarkg_samples +plotmarkg\$samples = as.character( '$gois.use.plotmarkg_samples' ) + #end if + #if $gois.use.plotexp_n +plotexp\$n = as.character( '$gois.use.plotexp_n' ) + #end if +plotmarkg\$cthr = as.integer( '$gois.use.plotmarkg_cthr' ) +plotmarkg\$order.cells = as.logical( '$gois.use.plotmarkg_order_cells' ) +plotmarkg\$aggr = as.logical( '$gois.use.plotmarkg_aggr' ) +plotmarkg\$norm = as.logical( '$gois.use.plotmarkg_norm' ) +plotmarkg\$cluster_cols = as.logical( '$gois.use.plotmarkg_cluster_cols' ) +plotmarkg\$cluster_rows = as.logical( '$gois.use.plotmarkg_cluster_rows' ) +plotmarkg\$cluster_set = as.logical( '$gois.use.plotmarkg_cluster_set' ) + #end if +#end if + +#if str($diffgtest.do_opt) == "yes": +plotdiffg = formals(plotdiffgenesnb) + +perform.diffgene = TRUE +plotdiffg\$Aname = '$diffgtest.set_a.name_set' +plotdiffg\$Bname = '$diffgtest.set_b.name_set' + +gfdat.A.use = list() +gfdat.B.use = list() + +gfdat.A.use\$manual = NULL +gfdat.A.use\$regex = NULL +gfdat.A.use\$cln = NULL +gfdat.B.use\$manual = NULL +gfdat.B.use\$regex = NULL +gfdat.B.use\$cln = NULL + + #if str($diffgtest.set_a.meth.type) == "cln": +gfdat.A.use\$cln = string2numericvector( '$diffgtest.set_a.meth.selector' ) + #else if str($diffgtest.set_a.meth.type) == "regex": +gfdat.A.use\$regex = as.character( '$diffgtest.set_a.meth.selector' ) + #else if str($diffgtest.set_a.meth.type) == "manual": +gfdat.A.use\$manual = string2textvector( '$diffgtest.set_a.meth.selector' ) + #end if + + #if str($diffgtest.set_b.meth.type) == "cln": +gfdat.B.use\$cln = string2numericvector( '$diffgtest.set_b.meth.selector' ) + #else if str($diffgtest.set_b.meth.type) == "regex": +gfdat.B.use\$regex = as.character( '$diffgtest.set_b.meth.selector' ) + #else if str($diffgtest.set_b.meth.type) == "manual": +gfdat.B.use\$manual = string2textvector( '$diffgtest.set_b.meth.selector' ) + #end if + + #if str($diffgtest.use.def) == "no": +plotdiffg\$pthr = as.numeric( '$diffgtest.use.plotdiffg_pthr' ) +plotdiffg\$padj = as.logical( '$diffgtest.use.plotdiffg_padj' ) +plotdiffg\$lthr = as.integer( '$diffgtest.use.plotdiffg_lthr' ) +plotdiffg\$show_names = as.logical( '$diffgtest.use.plotdiffg_show_names' ) + #if $diffgtest.use.plotdiffg_mthr +plotdiffg\$mthr = as.numeric( '$diffgtest.use.plotdiffg_mthr' ) + #end if + #end if +#end if + +out.pdf='${outpdf}' + ]]></token> + <token name="@TRAJECTORY_CHEETAH@"><![CDATA[ + +in.rdat = readRDS('${inputrds}') + +pstc.projc = formals(projcells) +pstc.projc\$knn = as.integer( '$projcell.knn' ) +pstc.projc\$cthr = as.integer( '$projcell.cthr' ) +#if str($projcell.use.def) == "no": +pstc.projc\$nmode = as.logical( '$projcell.use.nmode' ) +pstc.projc\$fr = as.logical( '$projcell.use.fr' ) +#end if + +pstc.projb = formals(projback) +pstc.projb\$pdishuf = as.integer( '$projback.pdishuf' ) +#if str($projback.use.def) == "no": +pstc.projb\$fast = as.logical( '$projback.use.fast' ) +pstc.projb\$rseed = as.integer( '$projback.use.rseed' ) +#end if + +pstc.comppval = formals(comppvalue) +pstc.comppval\$pthr = as.numeric( '$comppval.pthr' ) +pstc.comppval\$sensitive = as.logical( '$comppval.sensitive' ) + +pstc.plotgraph = formals(plotgraph) +pstc.plotgraph\$showCells = as.logical( '$plotgraph.showcells' ) +pstc.plotgraph\$scthr = as.numeric( '$plotgraph.scthr' ) +#if str($plotgraph.use.def) == "no": +pstc.plotgraph\$showTsne = as.logical( '$plotgraph.use.showtsne' ) +pstc.plotgraph\$tp = as.numeric( '$plotgraph.use.tp' ) +#end if + +pstc.compscore = formals(compscore) +pstc.compscore\$nn = as.integer( '$compscore.nn' ) +pstc.compscore\$scthr = as.numeric( '$compscore.scthr' ) + +out.pdf='${outpdf}' +out.rdat='${outrdat}' + ]]></token> + <token name="@INSPECTTRAJECTORIES_CHEETAH@"><![CDATA[ +in.rdat = readRDS('${inputrds}') + +perform.stemID = FALSE +perform.fateID = FALSE +perform.fateID.sominspect = FALSE + +#if str($trjsid.basic.doit) == "yes" +perform.stemID = TRUE +trjsid.getproj = formals(getproj) + +trjsid.numdiffgenes = 10 +trjsid.getproj\$i = as.integer( '$trjsid.basic.i' ) +trjsid.branchcells.ijk = string2numericvector( '$trjsid.basic.br' ) + + #if str($trjsid.basic.use.def) == "no": + +trjsid.getproj\$zscore = as.logical( '$trjsid.basic.use.zscore' ) +trjsid.numdiffgenes = as.integer( '$trjsid.basic.use.ndiffgenes' ) + + #end if +#end if + +#if str($trjfid.basic.doit) == "yes": +perform.fateID = TRUE +trjfid.cellsfrom = formals(cellsfromtree) +trjfid.filterset = formals(filterset) +trjfid.getsom = formals(getsom) +trjfid.procsom = formals(procsom) +trjfid.plotheat = list() + +trjfid.cellsfrom\$z = string2numericvector( '$trjfid.basic.cellsfromz' ) + #if str($trjfid.basic.use.def) == "no": + +trjfid.filterset\$minexpr = as.integer( '$trjfid.basic.use.filterset_minexpr' ) +trjfid.filterset\$minnumber = as.integer( '$trjfid.basic.use.filterset_minnumber' ) +trjfid.getsom\$nb = as.numeric( '$trjfid.basic.use.getsom_nb' ) +trjfid.getsom\$alpha = as.numeric( '$trjfid.basic.use.getsom_alpha' ) +trjfid.procsom\$corthr = as.numeric( '$trjfid.basic.use.procsom_corthr' ) +trjfid.procsom\$minsom = as.integer( '$trjfid.basic.use.procsom_minsom' ) +trjfid.plotheat\$xgrid = as.logical( '$trjfid.basic.use.plotheat_xgrid' ) +trjfid.plotheat\$ygrid = as.logical( '$trjfid.basic.use.plotheat_ygrid' ) +trjfid.plotheat\$xlab = as.logical( '$trjfid.basic.use.plotheat_xlab' ) + #end if + + #if str($trjfid.basic.som.doit) == "yes": +perform.fateID.sominspect = TRUE + +trjfidsomi = list() + #if str($trjfid.basic.som.use_genes.typer) == "genelist": +trjfidsomi.use.genes = string2textvector( '$trjfid.basic.som.use_genes.use_genes' ) + #else if str($trjfid.basic.som.use_genes.typer) == "cln": +trjfidsomi.use.genes = as.integer( '$trjfid.basic.som.use_genes.use_genes' ) + #end if + +trjfidsomi.use.types = '$trjfid.basic.som.use_types' +trjfidsomi\$name = '$trjfid.basic.som.title' + + #if str($trjfid.basic.som.use.def) == "no": +trjfidsomi\$cluster = as.logical( '$trjfid.basic.som.use.cluster' ) +trjfidsomi\$alpha = as.numeric( '$trjfid.basic.som.use.alpha' ) + #end if + #end if +#end if + +out.pdf = '${outpdf}' +out.diffgenes = '${outdiffgenes}' + +]]></token> +</macros>
--- a/macros_cluster.xml Mon Apr 15 17:53:28 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,404 +0,0 @@ -<macros> - <macro name="cluster_inputs" > - <param name="intable" type="data" format="tabular" label="Count Matrix" /> - <section name="filt" title="Filtering" expanded="true" > - <param name="mintotal" type="integer" min="1" value="3000" label="Min Transcripts" help="The minimum total transcripts required. Cells with less than mintotal transcripts are filtered out." /> - <param name="minexpr" type="integer" min="1" value="5" label="Min Expression" help="The minimum required transcript counts of a gene in the minimum number of cells (below)" /> - <param name="minnumber" type="integer" min="1" value="5" label="Min Cells" help="The minumum number of cells for gene expression to be counted" /> - <expand macro="use_defaults_no" > - <param name="knn" type="integer" min="0" value="10" label="K-nearest-neighbours" help="Number of nearest neighbors used to infer corresponding cell types in different batches" /> - <param name="CGenes" type="text" optional="true" label="CGenes" help="Filter out genes with correlated expression for cell type inference" > - <expand macro="sanitize_string_vector" /> - </param> - <param name="FGenes" type="text" optional="true" label="FGenes" help="Explicitly filter out genes for cell type inference" > - <expand macro="sanitize_string_vector" /> - </param> - <param name="LBatch_regexes" type="text" optional="true" label="Batch Regex" help="List of regexes to capture experimental batches for batch effect correction" > - <expand macro="sanitize_string_vector" /> - </param> - <param name="ccor" type="float" value="0.4" label="CCor" help="Correlation coefficient used as a threshold for determining correlated genes" /> - <param name="bmode" type="select" label="Batch Mode" help="Method to regress out batch effects" > - <option value="RaceID" selected="true" >RaceID</option> - <option value="scran">SCRAN</option> - </param> - <conditional name="ccc" > - <param name="use" type="select" label="Perform Cell-cycle correction?" > - <option value="yes" >Yes</option> - <option value="no" selected="true" >No</option> - </param> - <when value="no" /> - <when value="yes" > - <param name="vset" type="text" optional="true" label="List of Gene Sets" > - <expand macro="sanitize_string_vector" /> - </param> - <param name="pvalue" type="float" value="0.01" min="0" max="1" label="P-value Cutoff" help="P-value cutoff for determining enriched components" /> - <param name="quant" type="float" value="0.01" min="0" max="1" label="Quantification Fraction" help="Upper and lower fraction of gene loadings use for determining enriched components" /> - <param name="ncomp" type="integer" min="0" optional="true" label="Number of components to use" help="If left blank, the maximum number of components are used" /><!-- 0 = NULL --> - <param name="dimr" type="boolean" value="true" label="Derive Components from saturation criterion" /> - <param name="mode" type="select" label="Type of Component Analysis" help="If ICA is selected, ensure that the number of components value above is sufficiently high" > - <option value="pca" selected="true">PCA</option> - <option value="ica">ICA</option> - </param> - <param name="logscale" type="boolean" value="false" label="Log-transform data prior to PCA or ICA" help="" /> - </when> - </conditional> - </expand> - </section> - <section name="clust" title="Clustering" expanded="true" > - <!-- CompDist --> - <param name="metric" type="select" label="Distance Metric" > - <option value="pearson" selected="true" >Pearson</option> - <option value="spearman">Spearman</option> - <option value="logpearson">Log Pearson</option> - <option value="euclidean">Euclidean</option> - </param> - <!-- ClustExp --> - <param name="funcluster" type="select" label="Clustering method" > - <option value="kmedoids" selected="true" >K-medoids</option> - <option value="kmeans">K-means</option> - <option value="hclust">H-Clust</option> - </param> - <expand macro="use_defaults_no" > - <!-- CompDist --> - <param name="fselect" type="boolean" value="true" label="Perform feature selection" /> - <param name="knn" type="integer" min="0" optional="true" label="KNN" help="Number of nearest neighbours for imputing gene expression" /><!-- 0: NULL --> - <!-- ClustExp --> - <param name="sat" type="boolean" checked="true" label="Saturation-based clustering?" help="Determine number of clusters on saturation point of the mean within-cluster dispersion as a function of the cluster number." /> - <param name="clustnr" type="integer" min="0" value="30" label="Max number of clusters using Saturation-by-mean" help="Max number of clusters for the derivation of the cluster number by the saturation of mean within-cluster-dispersion." /> - <param name="samp" type="integer" min="0" optional="true" label="Sample random number of cells" help="Number of random sample of cells used for the inference of cluster number and Jaccard similarity" /><!-- 0:NULL --> - <param name="cln" type="integer" min="0" optional="true" label="Number of clusters" /><!-- 0:Null --> - <param name="bootnr" type="integer" min="0" value="50" label="Number of booststrapping runs" /> - <param name="rseed" type="integer" value="17000" label="Random seed" /> - </expand> - </section> - <section name="outlier" title="Outliers" expanded="true" > - <!-- Find Outliers --> - <param name="outminc" type="integer" min="0" value="5" label="Minimum Transcripts" help="minimal transcript count of a gene in a clusters to be tested for being an outlier gene" /> - <param name="outlg" type="integer" min="1" value="2" label="Minimum Genes" help="Minimum number of outlier genes required for being an outlier cell" /> - <!-- RFCorrect --> - <param name="final" type="boolean" value="true" label="Plot Final Clusters?" help="Reclassification of cell types using out-of-bag analysis is performed based on the final clusters after outlier identification. If 'FALSE', then the cluster partition prior to outlier identification is used for reclassification." /> - <expand macro="use_defaults_no" > - <!-- Find Outliers --> - <param name="probthr" type="float" min="0" value="0.001" label="Outlier Probability Threshold" help="Probability threshold for the above specified minimum number of genes to be an outlier cell. This probability is computed from a negative binomial background model of expression in a cluster" /> - <param name="outdistquant" type="float" min="0" max="1" value="0.95" label="Outlier Distance Quantile" help="Outlier cells are merged to outlier clusters if their distance smaller than the outdistquant-quantile of the distance distribution of pairs of cells in the orginal clusters after outlier removal" /> - <!-- RFCorrect --> - <param name="nbtree" type="integer" optional="true" label="Number of trees to be built" /><!-- 0:Null --> - <param name="nbfactor" type="integer" min="0" value="5" label="Tree Factor" help="Number of trees based on the number of cells multiplied by this factor. Effective only if the number of trees parameter is set to 0" /> - <param name="rfseed" type="integer" value="12345" label="Random Seed" /> - </expand> - </section> - <section name="tsne" title="tSNE and FR" expanded="true" > - <!-- CompTSNE --> - <param name="perplexity" type="integer" min="0" value="30" label="Perplexity" help="Perplexity of the t-SNE map" /> - <!-- CompFR --> - <param name="knn" type="integer" min="0" value="10" label="KNN" help="Number of nearest neighbours used for the inference of the Fruchterman-Rheingold layout" /> - <expand macro="use_defaults_no" > - <!-- CompTSNE --> - <param name="initial_cmd" type="boolean" checked="true" label="tSNE map initialised by classical multidimensional scaling" /> - <param name="rseed_tsne" type="integer" value="15555" label="Random Seed (tSNE)" /> - <!-- CompFR --> - <param name="rseed_fr" type="integer" min="0" value="15555" label="Random Seed (FR)" /> - </expand> - </section> - <section name="extra" title="Extra Parameters" expanded="false" > - <param name="tablelim" type="integer" min="1" value="25" label="Table Limit" help="Top N genes to print per cluster" /> - <param name="plotlim" type="integer" min="1" value="10" label="Plot Limit" help="Top N genes to plot. Must be less than or equal to the Table Limit" /> - <param name="foldchange" type="float" min="0" value="1" label="Fold change" /> - <param name="pvalue" type="float" min="0" max="1" value="0.01" label="P-value Cutoff" help="P-value cutoff for the inference of differential gene expression" /> - </section> - </macro> - <macro name="cluster_tests" > - <test> - <!-- default test --> - <conditional name="tool" > - <param name="mode" value="cluster" /> - <!-- This is a file with a single word 'test', which prompts the scripts to use the test intestinalData in the library --> - <param name="intable" value="use.intestinal" /> - </conditional> - <output name="outgenelist" value="intestinal.genelist" /> - <output name="outpdf" value="intestinal.pdf" compare="sim_size" delta="50" /> - </test> - <test> - <!-- defaults, feeding in a matrix with reduced filtering --> - <conditional name="tool" > - <param name="mode" value="cluster" /> - <param name="intable" value="matrix.tabular" /> - <section name="filt" > - <param name="mintotal" value="1000" /> - <param name="minexpr" value="1" /> - <param name="minnumber" value="3" /> - </section> - <param name="use_log" value="true" /> - <output name="outgenelist" value="matrix.genelist" /> - <output name="outrdat" value="matrix.rdat" compare="sim_size" delta="15" /> - <output name="outpdf" value="matrix.pdf" compare="sim_size" delta="10" /> - <output name="outlog" value="matrix.log" /> - </conditional> - </test> - <test> - <!-- defaults, but manually specified. No opts, no CC. Generates identical to above --> - <conditional name="tool" > - <param name="mode" value="cluster" /> - <param name="intable" value="use.intestinal" /> - <section name="filt" > - <param name="mintotal" value="3000" /> - <param name="minexpr" value="5" /> - <param name="minnumber" value="5" /> - <expand macro="test_nondef" > - <param name="knn" value="10" /> - <param name="ccor" value="0.4" /> - <param name="bmode" value="RaceID" /> - </expand> - </section> - <section name="clust" > - <param name="metric" value="pearson" /> - <param name="funcluster" value="kmedoids" /> - <expand macro="test_nondef" > - <param name="fselect" value="true" /> - <param name="sat" value="true" /> - <param name="clustnr" value="30" /> - <param name="bootnr" value="50" /> - <param name="rseed" value="17000" /> - </expand> - </section> - <section name="outlier" > - <param name="outminc" value="5" /> - <param name="outlg" value="2" /> - <param name="final" value="false" /> - <expand macro="test_nondef" section_name="outlier" > - <param name="probthr" value="0.001" /> - <param name="outdistquant" value="0.95" /> - <param name="rfseed" value="12345" /> - <param name="nbfactor" value="5" /> - </expand> - </section> - <section name="tsne" > - <param name="perplexity" value="30" /> - <param name="knn" value="10" /> - <expand macro="test_nondef" section_name="tsne" > - <param name="initial_cmd" value="true" /> - <param name="rseed_tsne" value="15555" /> - <param name="rfseed_fr" value="15555" /> - </expand> - </section> - </conditional> - <output name="outgenelist" value="intestinal.genelist" /> - <output name="outpdf" value="intestinal.pdf" compare="sim_size" delta="50" /> - </test> - <test> - <!-- Advanced. Opts, CC used --> - <conditional name="tool" > - <param name="mode" value="cluster" /> - <param name="intable" value="use.intestinal" /> - <section name="filt" > - <param name="mintotal" value="2000" /> - <param name="minexpr" value="3" /> - <param name="minnumber" value="2" /> - <expand macro="test_nondef" > - <param name="knn" value="5" /> - <param name="ccor" value="0.5" /> - <param name="CGenes" value="Gga3,Ggact,Ggct" /> - <param name="FGenes" value="Zxdc,Zyg11a,Zyg11b,Zyx" /> - <param name="LBatch_regexes" value="^I5,^II5,^III5,^IV5d,^V5d" /> - <param name="bmode" value="scran" /> - <conditional name="ccc" > - <param name="use" value="yes" /> - <param name="pvalue" value="0.05" /> - <param name="quant" value="0.05" /> - <param name="ncomp" value="3" /> - <param name="dimr" value="true" /> - <param name="mode" value="pca" /> - <param name="logscale" value="true" /> - </conditional> - </expand> - </section> - <section name="clust" > - <param name="metric" value="euclidean" /> - <param name="funcluster" value="hclust" /> - <expand macro="test_nondef" > - <param name="fselect" value="false" /> - <param name="knn" value="5" /> - <param name="sat" value="false" /> - <param name="samp" value="10" /> - <param name="cln" value="10" /> - <param name="clustnr" value="10" /> - <param name="bootnr" value="30" /> - <param name="rseed" value="17000" /> - </expand> - </section> - <section name="outlier" > - <param name="outminc" value="3" /> - <param name="outlg" value="5" /> - <param name="final" value="true" /> - <expand macro="test_nondef" > - <param name="probthr" value="0.01" /> - <param name="outdistquant" value="0.5" /> - <param name="rfseed" value="12345" /> - <param name="nbfactor" value="5" /> - <param name="nbtree" value="10" /> - </expand> - </section> - <section name="tsne" > - <param name="perplexity" value="20" /> - <param name="knn" value="6" /> - <expand macro="test_nondef" > - <param name="initial_cmd" value="false" /> - <param name="rseed_tsne" value="15555" /> - <param name="rfseed_fr" value="15555" /> - </expand> - </section> - </conditional> - <output name="outgenelist" value="intestinal_advanced.genelist" /> - <output name="outpdf" value="intestinal_advanced.pdf" compare="sim_size" delta="150" /> - </test> - </macro> - <token name="@FILTNORM_CHEETAH@"><![CDATA[ -## Perform do.filter -use.filtnormconf = TRUE - -## Perform do.cluster, do.outlier, do.clustmap, mkgenelist -use.cluster = FALSE - -in.table = read.table( - '${intable}', - stringsAsFactors = F, - na.strings=c("NA", "-", "?", "."), - sep='\t', - header=TRUE, - row.names=1 -) - -## Hidden flag to use test data instead -## see: test-data/use.intestinal - -use.test.data = (names(in.table)[1] == "test") - -sc = NULL -if (use.test.data) { - sc = SCseq(intestinalData) - message("Loading test data from library") -} else { - sc = SCseq(in.table) -} - - -filt = formals(filterdata) -filt.ccc = formals(CCcorrect) -filt.use.ccorrect = FALSE -filt.lbatch.regexes = NULL - -filt.geqone = as.logical( '$filt.hist_geq_one' ) -filt\$mintotal = as.integer( '$filt.mintotal' ) -filt\$minexpr = as.integer( '$filt.minexpr' ) -filt\$minnumber = as.integer( '$filt.minnumber' ) -#if str($filt.use.def) == "no": -filt\$knn = as.integer( '$filt.use.knn' ) -filt\$ccor = as.numeric( '$filt.use.ccor' ) -filt\$bmode = as.character( '$filt.use.bmode' ) - #if $filt.use.LBatch_regexes: -filt.lbatch.regexes = string2textvector( '$filt.use.LBatch_regexes' ) - #end if - #if $filt.use.CGenes: -filt\$CGenes = string2textvector( '$filt.use.CGenes' ) - #end if - #if $filt.use.FGenes: -filt\$FGenes = string2textvector( '$filt.use.FGenes' ) - #end if - #if str($filt.use.ccc.use) == "yes" -filt.use.ccorrect = TRUE - #if $filt.use.ccc.vset: -filt.ccc\$vset = string2textvector( '$filt.use.ccc.vset' ) - #end if - #if $filt.use.ccc.ncomp: -filt.ccc\$nComp = as.integer( '$filt.use.ccc.ncomp' ) - #end if -filt.ccc\$pvalue = as.numeric( '$filt.use.ccc.pvalue' ) -filt.ccc\$quant = as.numeric( '$filt.use.ccc.quant' ) -filt.ccc\$dimR = as.logical( '$filt.use.ccc.dimr' ) -filt.ccc\$mode = as.character( '$filt.use.ccc.mode.value' ) -filt.ccc\$logscale = as.logical( '$filt.use.ccc.logscale' ) - #end if -#end if - -out.pdf = '${outpdf}' -out.rdat = '${outrdat}' - -]]></token> - <token name="@CLUSTER_CHEETAH@"><![CDATA[ - -in.rdat = readRDS('${inputrds}') - -sc = in.rdat - -## Perform do.filter -use.filtnormconf = FALSE - -## Perform do.cluster, do.outlier, do.clustmap, mkgenelist -use.cluster = TRUE - - -clust.compdist = formals(compdist) -clust.clustexp = formals(clustexp) -clust.compdist\$metric = as.character( '$clust.metric' ) -clust.clustexp\$FUNcluster = as.character( '$clust.funcluster' ) - -#if str($clust.use.def) == "no": - -clust.compdist\$FSelect = as.logical( '$clust.use.fselect' ) - #if $clust.use.knn: -clust.compdist\$knn = as.integer( '$clust.use.knn' ) - #end if -clust.clustexp\$sat = as.logical( '$clust.use.sat' ) - #if $clust.use.samp: -clust.clustexp\$samp = as.integer( '$clust.use.samp' ) - #end if - #if $clust.use.cln: -clust.clustexp\$cln = as.integer( '$clust.use.cln' ) -clust.clustexp\$clustnr = as.integer( '$clust.use.clustnr' ) -clust.clustexp\$bootnr = as.integer( '$clust.use.bootnr' ) -##clust.clustexp\$rseed = as.integer( '$clust.use.rseed' ) - #end if -#end if - -outlier.use.randomforest = FALSE -outlier.findoutliers = formals(findoutliers) -outlier.clustheatmap = formals(clustheatmap) -outlier.rfcorrect = formals(rfcorrect) - -outlier.findoutliers\$outminc = as.integer( '$outlier.outminc' ) -outlier.findoutliers\$outlg = as.integer( '$outlier.outlg' ) -outlier.rfcorrect\$final = as.logical( '$outlier.final' ) - -#if str($outlier.use.def) == "no": - #if $outlier.use.nbtree: -outlier.rfcorrect\$nbtree = as.integer( '$outlier.use.nbtree' ) - #end if -outlier.findoutliers\$probthr = as.numeric( '$outlier.use.probthr' ) -outlier.findoutliers\$outdistquant = as.numeric( '$outlier.use.outdistquant' ) -##outlier.rfcorrect\$rfseed = as.integer( '$outlier.use.rfseed' ) -outlier.rfcorrect\$nbfactor = as.integer( '$outlier.use.nbfactor' ) -#end if - -cluster.comptsne = formals(comptsne) -cluster.compfr = formals(compfr) - -cluster.comptsne\$perplexity = as.integer( '$tsne.perplexity' ) -cluster.compfr\$knn = as.integer( '$tsne.knn' ) -#if str($tsne.use.def) == "no": -cluster.comptsne\$initial_cmd = as.logical( '$tsne.use.initial_cmd' ) -cluster.comptsne\$rseed = as.integer( '$tsne.use.rseed_tsne' ) -cluster.compfr\$rseed = as.integer( '$tsne.use.rseed_fr' ) -#end if - -genelist.tablelim = as.integer( '$extra.tablelim' ) -genelist.plotlim = as.integer( '$extra.plotlim' ) -genelist.foldchange = as.integer( '$extra.foldchange' ) -genelist.pvalue = as.numeric( '$extra.pvalue' ) - -out.pdf = '${outpdf}' -out.rdat = '${outrdat}' -out.genelist = '${outgenelist}' - -]]> - </token> -</macros>
--- a/macros_inspectclusters.xml Mon Apr 15 17:53:28 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,326 +0,0 @@ -<macros> - <macro name="inspectclusters_inputs" > - <param name="inputrds" type="data" format="rdata" label="Input RaceID RDS" help="Requires the RDS output from the cluster analysis" /> - <conditional name="plotgen"> - <param name="do_opt" type="select" label="Plot All Clusters?" help="Generates tSNE and F-R plots" > - <option value="yes" selected="true" >Yes</option> - <option value="no" >No</option> - </param> - <when value="no" /> - <when value="yes" /> - </conditional> - <conditional name="plotsym" > - <param name="do_opt" type="select" label="Perform Subset Analysis?" > - <expand macro="yesno_checkedno" /> - </param> - <when value="no" /> - <when value="yes"> - <param name="types_regex" type="text" optional="true" label="Types Regex (removal)" help="Regular expression to be used to remove portions of the cell names, in order to indicate cell type. e.g. if your cells are labelled as 'cdk5_1, cdk5_2, glow_1, glow_2', then a regex of '_\\d$' would trim off the ends and set the cell types to 'cdk5, cdk5, glow, glow'. " > - <expand macro="sanitize_regex" /> - </param> - <expand macro="use_defaults_no" > - <param name="subset_regex" type="text" optional="true" label="Subset of Types: Regex (select)" help="Regular expression to be used to select for cell types of interest. e.g. if we are only interested in 'cdk5' positive cells, then we would simply put 'cdk5' here. "> - <expand macro="sanitize_regex" /> - </param> - <param name="fr" type="boolean" checked="false" label="Output FR map instead of tSNE" help="Performs a Fruchterman-Rheingold projection instead of tSNE" /> - </expand> - </when> - </conditional> - <conditional name="gois" > - <param name="do_opt" type="select" label="Examine Genes of Interest" > - <expand macro="yesno_checkedno" /> - </param> - <when value="no" /> - <when value="yes"> - <param name="inspect_goi_genes" type="text" value="" label="Genes to Examine" > - <expand macro="sanitize_string_vector" /> - </param> - <param name="inspect_goi_cells" type="text" optional="true" label="Cells to Examine" help="Leave blank to examine all" ><!-- "":NULL --> - <expand macro="sanitize_string_vector" /> - </param> - <expand macro="use_defaults_no" > - <!-- Yes, examine genes of interest, but no do not use defaults... --> - <param name="inspect_goi_imputed" type="boolean" checked="false" label="Impute Genes" help="Only works if k-nearest neighbours was used in the clustering" /> - <param name="plotexp_n" type="text" optional="true" label="Title of plot" help="If blank, the list of genes are used" > - <expand macro="sanitize_title" /> - </param> - <param name="plotmarkg_cthr" type="integer" min="0" value="0" label="Cluster Threshold" help="Include only clusters with cells greater than this" /> - <param name="plotmarkg_cl" type="text" optional="true" label="Cluster List" help="List of cluster numbers to include. Leave blank to include all" > - <expand macro="sanitize_numeric_vector" /> - </param><!-- "":NULL --> - <param name="plotmarkg_order_cells" type="boolean" checked="false" label="Order Cells" help="Order heatmap by cell names, not by cluster" /> - <param name="plotmarkg_aggr" type="boolean" checked="false" label="Average Expression" help="Show only average expression for each cluster" /> - <param name="plotmarkg_norm" type="boolean" checked="false" label="Normalise Gene Expression " help="Normalise gene expression to 1 to depict genes on the same scale" /> - <param name="plotmarkg_flo" type="integer" max="-5" optional="true" label="Lower-bound for Gene Expression" /><!-- negative or null vals only for this seem to work--> - <param name="plotmarkg_cap" type="integer" min="1" optional="true" label="Upper-bound for Gene Expression" /> - <param name="plotmarkg_samples" type="text" optional="true" label="Regex to Select Sample Names For Each cell" > - <expand macro="sanitize_regex" /> - </param> - <param name="plotmarkg_cluster_cols" type="boolean" checked="false" label="Cluster columns" /> - <param name="plotmarkg_cluster_rows" type="boolean" checked="true" label="Cluster rows" /> - <param name="plotmarkg_cluster_set" type="boolean" checked="false" label="Order Clusters by Heirarchal Clustering of Cluster Medoids." /> - </expand> - </when> - </conditional> - <conditional name="diffgtest" > - <param name="do_opt" type="select" label="Differential Gene Testing" > - <expand macro="yesno_checkedno" /> - </param> - <when value="no" /> - <when value="yes" > - <expand macro="select_cells" sectionname="set_a" sectiontitle="Cells in Set A" /> - <expand macro="select_cells" sectionname="set_b" sectiontitle="Cells in Set B" /> - <expand macro="use_defaults_no" > - <param name="plotdiffg_pthr" type="float" min="0" max="1" value="0.05" label="P-value cutoff" help="Cutoff for displaying differentially expressed genes" /> - <param name="plotdiffg_padj" type="boolean" checked="true" label="Apply B-H Correction" help="Display genes with a Benjamini-Hochberg corrected false discovery rate lower than the above P-value threshold" /> - <param name="plotdiffg_lthr" type="float" min="0" value="0" label="Log2-Fold Threshold" help="Differentially expressed genes are displayed only for log2 fold-changes higher than this" /> - <param name="plotdiffg_mthr" type="float" optional="true" label="Log2 Mean Threshold" help="Differentially expressed genes are displayed only for log2 mean expression greater than this" /><!-- opt:-Inf --> - <param name="plotdiffg_show_names" type="boolean" checked="false" label="Display Gene Names" /> - </expand> - </when> - </conditional> - </macro> - <macro name="inspectclusters_tests" > - <test><!-- general plots, plotgen --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <conditional name="plotgen" > - <param name="do_opt" value="yes" /> - </conditional> - <output name="outpdf" value="out_cluster_general_default.pdf" /> - </test> - <test><!-- default test, plotsym --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <conditional name="plotsym" > - <param name="do_opt" value="yes" /> - <param name="types_regex" value="\\_\\d+" /> - <expand macro="test_nondef" > - <param name="subset_regex" value="IV|V" /> - </expand> - </conditional> - <conditional name="gois" > - <param name="do_opt" value="yes" /> - <param name="inspect_goi_genes" value="Apoa1,Apoa1bp,Apoa2,Apoa4,Apoa5" /> - </conditional> - <output name="outpdf" value="out_cluster_inspect_default.pdf" /> - </test> - <test><!-- Advanced 1 --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <conditional name="plotsym" > - <param name="do_opt" value="yes" /> - <param name="types_regex" value="\\_\\d+" /> - <expand macro="test_nondef" > - <param name="subset_regex" value="IV|V" /> - <param name="fr" value="true" /> - </expand> - </conditional> - <conditional name="gois" > - <param name="do_opt" value="yes" /> - <param name="inspect_goi_genes" value="Apoa1,Apoa1bp,Apoa2,Apoa4,Apoa5" /> - <expand macro="test_nondef" > - <param name="inspect_goi_imputed" value="false" /> - <param name="plotexp_n" value="Test Title" /> - <param name="plotmarkg_order_cells" value="true" /> - <param name="plotmarkg_aggr" value="true" /> - <param name="plotmarkg_norm" value="true" /> - <param name="plotmarkg_flo" value="-10" /> - <param name="plotmarkg_cap" value="100" /> - <param name="plotmarkg_samples" value="(\\_\\d+)$" /> - <param name="plotmarkg_cluster_cols" value="true" /> - <param name="plotmarkg_cluster_rows" value="false" /> - <param name="plotmarkg_cluster_set" value="true" /> - </expand> - </conditional> - <conditional name="diffgtest" > - <param name="do_opt" value="yes" /> - <section name="set_a" > - <param name="name_set" value="Test set A" /> - <conditional name="meth" > - <param name="type" value="regex" /> - <param name="selector" value="^V5.*" /> - </conditional> - </section> - <section name="set_b" > - <param name="name_set" value="Test set B" /> - <conditional name="meth" > - <param name="type" value="manual" /> - <param name="selector" value="I5d_3,I5d_4,I5d_6,I5d_8,I5d_9,I5d_10,I5d_11,I5d_12,I5d_13,I5d_14" /> - </conditional> - </section> - <expand macro="test_nondef" > - <param name="plotdiffg_pthr" value="0.1" /> - <param name="plotdiffg_padj" value="false" /> - <param name="plotdiffg_lthr" value="0.2" /> - <param name="plotdiffg_mthr" value="0.2" /> - <param name="plotdiffg_show_names" value="true" /> - </expand> - </conditional> - <output name="outpdf" value="out_cluster_inspect_advanced1.pdf" /> - </test> - <test><!-- Advanced 2 --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <conditional name="plotsym" > - <param name="do_opt" value="yes" /> - <expand macro="test_nondef" > - <param name="fr" value="false" /> - </expand> - </conditional> - <conditional name="gois" > - <param name="do_opt" value="yes" /> - <param name="inspect_goi_genes" value="Sp1,Spc24,Spcs1,Spcs2,Spcs3" /> - <expand macro="test_nondef" > - <param name="plotexp_n" value="Test Title 2" /> - <param name="plotmarkg_order_cells" value="true" /> - <param name="plotmarkg_aggr" value="true" /> - <param name="plotmarkg_norm" value="false" /> - <param name="plotmarkg_flo" value="-10" /> - <param name="plotmarkg_cap" value="10" /> - <param name="plotmarkg_cluster_cols" value="true" /> - <param name="plotmarkg_cluster_rows" value="true" /> - <param name="plotmarkg_cluster_set" value="true" /> - </expand> - </conditional> - <conditional name="diffgtest" > - <param name="do_opt" value="yes" /> - <section name="set_a" > - <param name="name_set" value="Test set A" /> - <conditional name="meth" > - <param name="type" value="regex" /> - <param name="selector" value="^IV.*" /> - </conditional> - </section> - <section name="set_b" > - <param name="name_set" value="Test set B" /> - <conditional name="meth" > - <param name="type" value="cln" /> - <param name="selector" value="3" /> - </conditional> - </section> - <expand macro="test_nondef" > - <param name="plotdiffg_pthr" value="0.8" /> - <param name="plotdiffg_padj" value="true" /> - <param name="plotdiffg_lthr" value="0.8" /> - <param name="plotdiffg_mthr" value="0.8" /> - <param name="plotdiffg_show_names" value="false" /> - </expand> - </conditional> - <output name="outpdf" value="out_cluster_inspect_advanced2.pdf" /> - </test> - </macro> - <token name="@INSPECTCLUSTERS_CHEETAH@"><![CDATA[ - -in.rdat = readRDS('${inputrds}') - -perform.plotting = FALSE -perform.symbolmap = FALSE -perform.genesofinterest = FALSE -perform.diffgene = FALSE - -#if str($plotgen.do_opt) == "yes": -perform.plotting = TRUE -#end if - -#if str($plotsym.do_opt) == "yes": -perform.symbolmap = TRUE -plotsym = formals(plotsymbolsmap) -plotsym.use.typeremoveregex = NULL -plotsym.use.typeremoveregex.subselect = NULL -plotsym.use.typeremoveregex = as.character( '$plotsym.types_regex' ) - #if str($plotsym.use.def) == "no": -plotsym\$fr = as.logical( '$plotsym.use.fr' ) - #if $plotsym.use.subset_regex: -plotsym.use.typeremoveregex.subselect = as.character( '$plotsym.use.subset_regex' ) - #end if - #end if -#end if - -#if str($gois.do_opt) == "yes": -perform.genesofinterest = TRUE -plotexp = formals(plotexpmap) -plotmarkg = formals(plotmarkergenes) - -inspect.goi.use.genes = string2textvector( '$gois.inspect_goi_genes' ) -plotexp\$g = inspect.goi.use.genes -plotmarkg\$genes = inspect.goi.use.genes - - #if $gois.inspect_goi_cells: -inspect.goi.use.cells = as.character( '$gois.inspect_goi_cells' ) -plotexp\$cells = inspect.goi.use.cells -plotmarkg\$cells = inspect.goi.use.cells - #end if - #if str($gois.use.def) == "no": -inspect.goi.use.imputed = as.logical( '$gois.use.inspect_goi_imputed' ) -plotexp\$imputed = inspect.goi.use.imputed -plotmarkg\$imputed = inspect.goi.use.imputed - #if $gois.use.plotmarkg_cl -plotmarkg\$cl = string2numericvector( '$gois.use.plotmarkg_cl' ) - #end if - #if $gois.use.plotmarkg_cap -plotmarkg\$cap = as.integer( '$gois.use.plotmarkg_cap' ) - #end if - #if $gois.use.plotmarkg_flo -plotmarkg\$flo = as.integer( '$gois.use.plotmarkg_flo' ) - #end if - #if $gois.use.plotmarkg_samples -plotmarkg\$samples = as.character( '$gois.use.plotmarkg_samples' ) - #end if - #if $gois.use.plotexp_n -plotexp\$n = as.character( '$gois.use.plotexp_n' ) - #end if -plotmarkg\$cthr = as.integer( '$gois.use.plotmarkg_cthr' ) -plotmarkg\$order.cells = as.logical( '$gois.use.plotmarkg_order_cells' ) -plotmarkg\$aggr = as.logical( '$gois.use.plotmarkg_aggr' ) -plotmarkg\$norm = as.logical( '$gois.use.plotmarkg_norm' ) -plotmarkg\$cluster_cols = as.logical( '$gois.use.plotmarkg_cluster_cols' ) -plotmarkg\$cluster_rows = as.logical( '$gois.use.plotmarkg_cluster_rows' ) -plotmarkg\$cluster_set = as.logical( '$gois.use.plotmarkg_cluster_set' ) - #end if -#end if - -#if str($diffgtest.do_opt) == "yes": -plotdiffg = formals(plotdiffgenesnb) - -perform.diffgene = TRUE -plotdiffg\$Aname = '$diffgtest.set_a.name_set' -plotdiffg\$Bname = '$diffgtest.set_b.name_set' - -gfdat.A.use = list() -gfdat.B.use = list() - -gfdat.A.use\$manual = NULL -gfdat.A.use\$regex = NULL -gfdat.A.use\$cln = NULL -gfdat.B.use\$manual = NULL -gfdat.B.use\$regex = NULL -gfdat.B.use\$cln = NULL - - #if str($diffgtest.set_a.meth.type) == "cln": -gfdat.A.use\$cln = c( '$diffgtest.set_a.meth.selector' ) - #else if str($diffgtest.set_a.meth.type) == "regex": -gfdat.A.use\$regex = as.character( '$diffgtest.set_a.meth.selector' ) - #else if str($diffgtest.set_a.meth.type) == "manual": -gfdat.A.use\$manual = string2textvector( '$diffgtest.set_a.meth.selector' ) - #end if - - #if str($diffgtest.set_b.meth.type) == "cln": -gfdat.B.use\$cln = c( '$diffgtest.set_b.meth.selector' ) - #else if str($diffgtest.set_b.meth.type) == "regex": -gfdat.B.use\$regex = as.character( '$diffgtest.set_b.meth.selector' ) - #else if str($diffgtest.set_b.meth.type) == "manual": -gfdat.B.use\$manual = string2textvector( '$diffgtest.set_b.meth.selector' ) - #end if - - #if str($diffgtest.use.def) == "no": -plotdiffg\$pthr = as.numeric( '$diffgtest.use.plotdiffg_pthr' ) -plotdiffg\$padj = as.logical( '$diffgtest.use.plotdiffg_padj' ) -plotdiffg\$lthr = as.integer( '$diffgtest.use.plotdiffg_lthr' ) -plotdiffg\$show_names = as.logical( '$diffgtest.use.plotdiffg_show_names' ) - #if $diffgtest.use.plotdiffg_mthr -plotdiffg\$mthr = as.numeric( '$diffgtest.use.plotdiffg_mthr' ) - #end if - #end if -#end if - -out.pdf='${outpdf}' - ]]> - </token> -</macros>
--- a/macros_inspecttrajectories.xml Mon Apr 15 17:53:28 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,243 +0,0 @@ -<macros> - <macro name="inspecttrajectories_inputs" > - <param name="inputrds" type="data" format="rdata" label="Input RDS" help="Requires the RDS output from the trajectory analysis" /> - <section name="trjsid" title="StemID Branch Link Examine" expanded="true" help="StemID derives cell lineage trees and predicts multipotent cell identities"> - <conditional name="basic" > - <param name="doit" type="select" label="Perform StemID?" > - <expand macro="yesno_checkedno" /> - </param> - <when value="no" /> - <when value="yes" > - <param name="i" type="integer" min="1" value="1" label="Cluster Number" /> - <param name="br" type="text" label="Trajectory Path i, j, k" help="A path of three clusters starting at i, passing through j, and ending at k." > - <expand macro="sanitize_numeric_vector" /> - </param> - <expand macro="use_defaults_no"> - <param name="zscore" type="boolean" checked="false" label="Plot Z-score Transformed Projections" /> - <param name="ndiffgenes" type="integer" min="1" value="10" label="Number of DE Genes" help="Number of differentially expressed genes to output per cluster" /> - </expand> - </when> - </conditional> - </section> - <section name="trjfid" title="FateID Branch Link Examine" expanded="true" help="FateID infers cell fate bias in multipotent progenitor cells" > - <conditional name="basic" > - <param name="doit" type="select" label="Perform FateID?" > - <expand macro="yesno_checkedno" /> - </param> - <when value="no" /> - <when value="yes" > - <param name="cellsfromz" type="text" value="" label="Cells from Clusters" help="Vector of valid cluster numbers ordered along the trajectory" > - <expand macro="sanitize_numeric_vector" /> - </param> - <expand macro="use_defaults_no"> - <param name="filterset_minexpr" type="integer" min="0" value="2" label="Min Expression" help="Minimum expression required for at least minnumber cells" /> - <param name="filterset_minnumber" type="integer" min="0" value="1" label="Min Number of Cells" help="Minimum number of cells in which a gene needs to be expressed at least a level of minexpr."/> - <param name="getsom_nb" type="integer" min="1" value="1000" label="SOM Nodes" help="Number of nodes of the self-organizing map." /> - <param name="getsom_alpha" type="float" min="0" value="0.5" label="Smoothing parameter" help="Pseudo-temporal expression profiles are derived by a local regression of expression values across the ordered cells using the function 'loess' from the package 'stats'. This is the parameter, which controls the degree of smoothing. Larger values return smoother profiles." /> - <param name="procsom_corthr" type="float" min="0" max="1" value="0.85" label="Correlation threshold" help="The z-score of the average normalized pseudo-temporal expression profiles within each node of the self-organizing map is computed, and the correlation of these z-scores between neighbouring nodes is computed. If the correlation is greater than 'corthr', neighbouring nodes are merged" /> - <param name="procsom_minsom" type="integer" min="0" value="3" label="Min SOM" help="Nodes of the self-organizing map with less than this number of transcripts are discarded" /> - <param name="plotheat_xgrid" type="boolean" checked="false" label="Partitioning along the x-axis" /> - <param name="plotheat_ygrid" type="boolean" checked="false" label="Partitioning along the y-axis" /> - <param name="plotheat_xlab" type="boolean" checked="false" label="Average position is given for each partition along the x-axis" /> - </expand> - <conditional name="som" > - <param name="doit" type="select" label="Perform Additional FateID Analysis with Self-Organised Map?" > - <expand macro="yesno_checkedno" /> - </param> - <when value="no" /> - <when value="yes" > - <conditional name="use_genes" > - <param name="typer" type="select" label="Genes to Inspect" > - <option value="genelist">List of Genes</option> - <option value="cln">Cluster</option> - </param> - <when value="genelist"> - <param name="use_genes" type="text" value="" label="Gene List" > - <expand macro="sanitize_string_vector" /> - </param> - </when> - <when value="cln" > - <param name="use_genes" type="integer" min="1" value="1" label="Cluster Number" /> - </when> - </conditional> - <param name="use_types" type="text" value="\\_\\d+" label="Types Regex (select)" help="Regex to select types across cell names" > - <expand macro="sanitize_regex" /> - </param> - <param name="title" type="text" value="Title" label="Plot title" /> - <expand macro="use_defaults_no"> - <param name="cluster" type="boolean" checked="false" label="Partitioning along the x-axis" /> - <param name="alpha" type="float" min="0" max="1" value="0.5" label="Smoothing parameter" /> - </expand> - </when> - </conditional> - </when> - </conditional> - </section> - </macro> - <macro name="inspecttrajectories_tests" > - <test> - <!-- default, generates a blank report --> - <param name="inputrds" value="out_traject_default.ltree.rdat" /> - <output name="outpdf" value="out_traject_inspect_default.ltree.blank.pdf" compare="sim_size" delta="15" /> - </test> - <test> - <!-- stemID branch inspection: vignette search "getproj" --> - <param name="inputrds" value="out_traject_default.ltree.rdat" /> - <section name="trjsid" > - <conditional name="basic" > - <param name="doit" value="yes" /> - <param name="i" value="3" /> - <param name="br" value="1,3,8" /> - </conditional> - </section> - <output name="outpdf" value="out_traject_inspect_stemid.pdf" compare="sim_size" delta="15" /> - <output name="outdiffgenes" value="out_traject_inspect_stemid.tabular" /> - </test> - <test> - <!-- fateID trajectory inspection: vignette search "cellsfromtree" --> - <param name="inputrds" value="out_traject_default.ltree.rdat" /> - <section name="trjfid" > - <conditional name="basic" > - <param name="doit" value="yes" /> - <param name="cellsfromz" value="2,1,4" /> - <expand macro="test_nondef" > - <param name="filterset_minexpr" value="2" /> - <param name="getsom_nb" value="1000" /> - <param name="getsom_alpha" value="0.5" /> - <param name="plotheat_xlab" value="true" /> - </expand> - </conditional> - </section> - <output name="outpdf" value="out_traject_inspect_fateid.pdf" compare="sim_size" delta="15" /> - </test> - <test> - <!-- fateID trajectory inspection with som: vignette search "SOM" --> - <param name="inputrds" value="out_traject_default.ltree.rdat" /> - <section name="trjfid" > - <conditional name="basic" > - <param name="doit" value="yes" /> - <param name="cellsfromz" value="2,1,4" /> - <conditional name="som" > - <param name="doit" value="yes" /> - <expand macro="test_nondef" > - <param name="typer" value="cln" /> - <param name="use_genes" value="12" /> - </expand> - </conditional> - </conditional> - </section> - <output name="outpdf" value="out_traject_inspect_fateid_som.pdf" compare="sim_size" delta="15" /> - </test> - <test> - <!-- uses all 3 sections with additional non-default params --> - <param name="inputrds" value="out_traject_default.ltree.rdat" /> - <section name="trjsid" > - <conditional name="basic" > - <param name="doit" value="yes" /> - <param name="i" value="5" /> - <param name="br" value="6,5,3" /> - <expand macro="test_nondef" > - <param name="zscore" value="true" /> - <param name="ndiffgenes" value="14" /> - </expand> - </conditional> - </section> - <section name="trjfid" > - <conditional name="basic" > - <param name="doit" value="yes" /> - <param name="cellsfromz" value="3,9,2" /> - <expand macro="test_nondef" > - <param name="filterset_minexpr" value="1" /> - <param name="filterset_minnumber" value="2" /> - <param name="procsom_minsom" value="5" /> - <param name="procsom_corthr" value="0.5" /> - <param name="plotheat_ygrid" value="true" /> - </expand> - <conditional name="som" > - <param name="doit" value="yes" /> - <conditional name="use_genes" > - <param name="typer" value="genelist" /> - <param name="use_genes" value="Clca4,Mki67,Defa24,Ybx1,Vasp,Apoa1" /> - <expand macro="test_nondef" > - <param name="cluster" value="true" /> - <param name="alpha" value="0.1" /> - </expand> - </conditional> - </conditional> - </conditional> - </section> - <output name="outdiffgenes" value="out_traject_inspect_allthree.tabular" /> - <output name="outpdf" value="out_traject_inspect_allthree.pdf" compare="sim_size" delta="15" /> - </test> - </macro> - <token name="@INSPECTTRAJECTORIES_CHEETAH@"><![CDATA[ -in.rdat = readRDS('${inputrds}') - -perform.stemID = FALSE -perform.fateID = FALSE -perform.fateID.sominspect = FALSE - -#if str($trjsid.basic.doit) == "yes" -perform.stemID = TRUE -trjsid.getproj = formals(getproj) - -trjsid.numdiffgenes = 10 -trjsid.getproj\$i = as.integer( '$trjsid.basic.i' ) -trjsid.branchcells.ijk = string2numericvector( '$trjsid.basic.br' ) - - #if str($trjsid.basic.use.def) == "no": - -trjsid.getproj\$zscore = as.logical( '$trjsid.basic.use.zscore' ) -trjsid.numdiffgenes = as.integer( '$trjsid.basic.use.ndiffgenes' ) - - #end if -#end if - -#if str($trjfid.basic.doit) == "yes": -perform.fateID = TRUE -trjfid.cellsfrom = formals(cellsfromtree) -trjfid.filterset = formals(filterset) -trjfid.getsom = formals(getsom) -trjfid.procsom = formals(procsom) -trjfid.plotheat = list() - -trjfid.cellsfrom\$z = string2numericvector( '$trjfid.basic.cellsfromz' ) - #if str($trjfid.basic.use.def) == "no": - -trjfid.filterset\$minexpr = as.integer( '$trjfid.basic.use.filterset_minexpr' ) -trjfid.filterset\$minnumber = as.integer( '$trjfid.basic.use.filterset_minnumber' ) -trjfid.getsom\$nb = as.numeric( '$trjfid.basic.use.getsom_nb' ) -trjfid.getsom\$alpha = as.numeric( '$trjfid.basic.use.getsom_alpha' ) -trjfid.procsom\$corthr = as.numeric( '$trjfid.basic.use.procsom_corthr' ) -trjfid.procsom\$minsom = as.integer( '$trjfid.basic.use.procsom_minsom' ) -trjfid.plotheat\$xgrid = as.logical( '$trjfid.basic.use.plotheat_xgrid' ) -trjfid.plotheat\$ygrid = as.logical( '$trjfid.basic.use.plotheat_ygrid' ) -trjfid.plotheat\$xlab = as.logical( '$trjfid.basic.use.plotheat_xlab' ) - #end if - - #if str($trjfid.basic.som.doit) == "yes": -perform.fateID.sominspect = TRUE - -trjfidsomi = list() - #if str($trjfid.basic.som.use_genes.typer) == "genelist": -trjfidsomi.use.genes = string2textvector( '$trjfid.basic.som.use_genes.use_genes' ) - #else if str($trjfid.basic.som.use_genes.typer) == "cln": -trjfidsomi.use.genes = as.integer( '$trjfid.basic.som.use_genes.use_genes' ) - #end if - -trjfidsomi.use.types = '$trjfid.basic.som.use_types' -trjfidsomi\$name = '$trjfid.basic.som.title' - - #if str($trjfid.basic.som.use.def) == "no": -trjfidsomi\$cluster = as.logical( '$trjfid.basic.som.use.cluster' ) -trjfidsomi\$alpha = as.numeric( '$trjfid.basic.som.use.alpha' ) - #end if - #end if -#end if - -out.pdf = '${outpdf}' -out.diffgenes = '${outdiffgenes}' - -]]> - </token> -</macros>
--- a/macros_trajectory.xml Mon Apr 15 17:53:28 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,156 +0,0 @@ -<macros> - <macro name="trajectory_inputs" > - <param name="inputrds" type="data" format="rdata" label="Input RDS" help="This is the output RData file given by the Cluster stage" /> - <section name="projcell" title="Compute transcriptome entropy of each cell" expanded="true" > - <param name="knn" type="integer" min="0" value="3" label="KNN" help="Number of nearest neighbours to consider." /> - <param name="cthr" type="integer" min="0" value="5" label="Cluster Threshold" help="Clusters to be included in the StemID analysis must have more than this number of cells" /> - <expand macro="use_defaults_no" > - <param name="nmode" type="boolean" checked="true" label="Neighbour Mode" help="Assign a cell of a given cluster to the cluster with the smallest average distance of the knn nearest neighbours within this cluster." /> - <param name="fr" type="boolean" checked="false" label="F-R layout?" help="Use Fruchterman-Rheigold layout instead of t-SNE" /> - </expand> - </section> - <section name="projback" title="Compute Cell Projections for Randomized Background Distribution" expanded="true" > - <param name="pdishuf" type="integer" min="100" value="2000" label="Shuffle Projections" help="Number of randomizations of cell positions for which to compute projections of cells on inter-cluster links" /> - <expand macro="use_defaults_no" > - <param name="fast" type="boolean" checked="false" label="Use a fast approximate background model to infer links" /> - <param name="rseed" type="integer" value="17000" label="Random seed" /> - </expand> - </section> - <section name="comppval" title="Computing P-values for Link Significance" expanded="true" > - <param name="pthr" type="float" value="0.01" label="P-value cutoff" help="Cutoff for link significance. This threshold is applied for the calculation of link scores reflecting how uniformly a link is occupied by cells." /> - <param name="sensitive" type="boolean" checked="false" label="Sensitivity" help="All cells on the most highly significant link are and the link itself are disregard to test significance of the remaining links with a binomial p-value" /> - </section> - <section name="plotgraph" title="StemID2 Lineage Graph" expanded="true" > - <param name="showcells" type="boolean" checked="false" label="Show Cells" help="Projections of cells are shown in the plot" /> - <param name="scthr" type="float" min="0" max="1" value="0" label="Score Threshold" help="Score threshold for links to be shown in the graph" /> - <expand macro="use_defaults_no"> - <param name="showtsne" type="boolean" checked="true" label="Show tSNE" help="Show transparent t-SNE map of cells in the background" /> - <param name="tp" type="float" min="0" max="1" value="0.5" label="Transparency" help="Level of transparency of t-SNE map" /> - </expand> - </section> - <section name="compscore" title="Compute StemID2 score" expanded="true" > - <param name="nn" type="integer" min="0" value="1" label="Number of Neighbours" help="Number of higher order neighbours to be included for the determination of links" /> - <param name="scthr" type="float" min="0" max="1" value="0" label="Score Threshold" help="Score threshold for links to be included in the calculation" /> - </section> - </macro> - <macro name="trajectory_tests" > - <test> - <!-- All defaults --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <output name="outpdf" value="out_traject_default.pdf" compare="sim_size" delta="15" /> - <output name="outrdat" value="out_traject_default.ltree.rdat" compare="sim_size" delta="15" /> - </test> - <test> - <!-- Manually specify the defaults: should ideally give same output as above test --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <section name="projcell" > - <param name="knn" value="3" /> - <param name="cthr" value="5" /> - <expand macro="test_nondef" > - <param name="nmode" value="true" /> - <param name="fr" value="false" /> - </expand> - </section> - <section name="projback" > - <param name="pdishuf" value="2000" /> - <expand macro="test_nondef" > - <param name="fast" value="false" /> - <param name="rseed" value="17000" /> - </expand> - </section> - <section name="comppval" > - <param name="pthr" value="0.01" /> - <param name="sensitive" value="false" /> - </section> - <section name="plotgraph" > - <param name="showcells" value="false" /> - <param name="scthr" value="0" /> - <expand macro="test_nondef" > - <param name="showtsne" value="true" /> - <param name="tp" value="0.5" /> - </expand> - </section> - <section name="compscore" > - <param name="nn" value="1" /> - <param name="scthr" value="0" /> - </section> - <output name="outpdf" value="out_traject_default.pdf" compare="sim_size" delta="15" /> - <output name="outrdat" value="out_traject_default.ltree.rdat" compare="sim_size" delta="15" /> - </test> - <test> - <!-- Non-defaults --> - <param name="inputrds" value="out_cluster_default.rdat" /> - <section name="projcell"> - <param name="knn" value="4" /> - <param name="cthr" value="4" /> - <expand macro="test_nondef" > - <param name="nmode" value="false" /> - <param name="fr" value="true" /> - </expand> - </section> - <section name="projback"> - <param name="pdishuf" value="1500" /> - <expand macro="test_nondef" > - <param name="fast" value="true" /> - <param name="rseed" value="1500" /> - </expand> - </section> - <section name="comppval"> - <param name="pthr" value="0.05" /> - <param name="sensitive" value="true" /> - </section> - <section name="plotgraph"> - <param name="showcells" value="true" /> - <param name="scthr" value="0.5" /> - <expand macro="test_nondef" > - <param name="showtsne" value="false" /> - <param name="tp" value="0.25" /> - </expand> - </section> - <section name="compscore"> - <param name="nn" value="2" /> - <param name="scthr" value="0.5" /> - </section> - <output name="outpdf" value="out_traject_adv_nondef.pdf" compare="sim_size" delta="15" /> - </test> - </macro> - <token name="@TRAJECTORY_CHEETAH@"><![CDATA[ - -in.rdat = readRDS('${inputrds}') - -pstc.projc = formals(projcells) -pstc.projc\$knn = as.integer( '$projcell.knn' ) -pstc.projc\$cthr = as.integer( '$projcell.cthr' ) -#if str($projcell.use.def) == "no": -pstc.projc\$nmode = as.logical( '$projcell.use.nmode' ) -pstc.projc\$fr = as.logical( '$projcell.use.fr' ) -#end if - -pstc.projb = formals(projback) -pstc.projb\$pdishuf = as.integer( '$projback.pdishuf' ) -#if str($projback.use.def) == "no": -pstc.projb\$fast = as.logical( '$projback.use.fast' ) -pstc.projb\$rseed = as.integer( '$projback.use.rseed' ) -#end if - -pstc.comppval = formals(comppvalue) -pstc.comppval\$pthr = as.numeric( '$comppval.pthr' ) -pstc.comppval\$sensitive = as.logical( '$comppval.sensitive' ) - -pstc.plotgraph = formals(plotgraph) -pstc.plotgraph\$showCells = as.logical( '$plotgraph.showcells' ) -pstc.plotgraph\$scthr = as.numeric( '$plotgraph.scthr' ) -#if str($plotgraph.use.def) == "no": -pstc.plotgraph\$showTsne = as.logical( '$plotgraph.use.showtsne' ) -pstc.plotgraph\$tp = as.numeric( '$plotgraph.use.tp' ) -#end if - -pstc.compscore = formals(compscore) -pstc.compscore\$nn = as.integer( '$compscore.nn' ) -pstc.compscore\$scthr = as.numeric( '$compscore.scthr' ) - -out.pdf='${outpdf}' -out.rdat='${outrdat}' - ]]> - </token> -</macros>
--- a/raceid_inspecttrajectories.xml Mon Apr 15 17:53:28 2019 -0400 +++ b/raceid_inspecttrajectories.xml Mon Aug 12 13:03:03 2019 -0400 @@ -1,30 +1,85 @@ -<tool id="raceid_inspecttrajectory" name="Lineage Branch Analysis using StemID" version="@VERSION_RACEID@.@VERSION_PACKAGE@.1" > +<tool id="raceid_inspecttrajectory" name="Lineage Branch Analysis using StemID" version="@VERSION_RACEID@.@VERSION_WRAPPER@" > <description>inspects branches of a lineage tree</description> <macros> <import>macros.xml</import> - <import>macros_inspecttrajectories.xml</import> + <import>macros_cheetah.xml</import> </macros> <expand macro="requirements" /> - <version_command><![CDATA[ -Rscript '$__tool_directory__/scripts/trajectoryinspect.R' @GET_VERSION@ -]]></version_command> - - <command detect_errors="exit_code"><![CDATA[ -#set bin = 'trajectoryinspect.R' -Rscript '$__tool_directory__/scripts/$bin' '$userconf' 2>&1 > '$outlog' - ]]></command> - - <configfiles> - <configfile name="userconf" ><![CDATA[ -@STRING2VECTOR@ - -@INSPECTTRAJECTORIES_CHEETAH@ -]]> - </configfile> - </configfiles> - + <expand macro="version_command_config" prog="trajectoryinspect.R" cheetah="INSPECTTRAJECTORIES_CHEETAH" /> <inputs> - <expand macro="inspecttrajectories_inputs" /> + <param name="inputrds" type="data" format="rdata" label="Input RDS" help="Requires the RDS output from the trajectory analysis" /> + <section name="trjsid" title="StemID Branch Link Examine" expanded="true" help="StemID derives cell lineage trees and predicts multipotent cell identities"> + <conditional name="basic" > + <param name="doit" type="select" label="Perform StemID?" > + <expand macro="yesno_checkedno" /> + </param> + <when value="no" /> + <when value="yes" > + <param name="i" type="integer" min="1" value="1" label="Cluster Number" /> + <param name="br" type="text" label="Trajectory Path i, j, k" help="A path of three clusters starting at i, passing through j, and ending at k." > + <expand macro="sanitize_numeric_vector" /> + </param> + <expand macro="use_defaults_no"> + <param name="zscore" type="boolean" checked="false" label="Plot Z-score Transformed Projections" /> + <param name="ndiffgenes" type="integer" min="1" value="10" label="Number of DE Genes" help="Number of differentially expressed genes to output per cluster" /> + </expand> + </when> + </conditional> + </section> + <section name="trjfid" title="FateID Branch Link Examine" expanded="true" help="FateID infers cell fate bias in multipotent progenitor cells" > + <conditional name="basic" > + <param name="doit" type="select" label="Perform FateID?" > + <expand macro="yesno_checkedno" /> + </param> + <when value="no" /> + <when value="yes" > + <param name="cellsfromz" type="text" value="" label="Cells from Clusters" help="Vector of valid cluster numbers ordered along the trajectory" > + <expand macro="sanitize_numeric_vector" /> + </param> + <expand macro="use_defaults_no"> + <param name="filterset_minexpr" type="integer" min="0" value="2" label="Min Expression" help="Minimum expression required for at least minnumber cells" /> + <param name="filterset_minnumber" type="integer" min="0" value="1" label="Min Number of Cells" help="Minimum number of cells in which a gene needs to be expressed at least a level of minexpr."/> + <param name="getsom_nb" type="integer" min="1" value="1000" label="SOM Nodes" help="Number of nodes of the self-organizing map." /> + <param name="getsom_alpha" type="float" min="0" value="0.5" label="Smoothing parameter" help="Pseudo-temporal expression profiles are derived by a local regression of expression values across the ordered cells using the function 'loess' from the package 'stats'. This is the parameter, which controls the degree of smoothing. Larger values return smoother profiles." /> + <param name="procsom_corthr" type="float" min="0" max="1" value="0.85" label="Correlation threshold" help="The z-score of the average normalized pseudo-temporal expression profiles within each node of the self-organizing map is computed, and the correlation of these z-scores between neighbouring nodes is computed. If the correlation is greater than 'corthr', neighbouring nodes are merged" /> + <param name="procsom_minsom" type="integer" min="0" value="3" label="Min SOM" help="Nodes of the self-organizing map with less than this number of transcripts are discarded" /> + <param name="plotheat_xgrid" type="boolean" checked="false" label="Partitioning along the x-axis" /> + <param name="plotheat_ygrid" type="boolean" checked="false" label="Partitioning along the y-axis" /> + <param name="plotheat_xlab" type="boolean" checked="false" label="Average position is given for each partition along the x-axis" /> + </expand> + <conditional name="som" > + <param name="doit" type="select" label="Perform Additional FateID Analysis with Self-Organised Map?" > + <expand macro="yesno_checkedno" /> + </param> + <when value="no" /> + <when value="yes" > + <conditional name="use_genes" > + <param name="typer" type="select" label="Genes to Inspect" > + <option value="genelist">List of Genes</option> + <option value="cln">Cluster</option> + </param> + <when value="genelist"> + <param name="use_genes" type="text" value="" label="Gene List" > + <expand macro="sanitize_string_vector" /> + </param> + </when> + <when value="cln" > + <param name="use_genes" type="integer" min="1" value="1" label="Cluster Number" /> + </when> + </conditional> + <param name="use_types" type="text" value="\\_\\d+" label="Types Regex (select)" help="Regex to select types across cell names" > + <expand macro="sanitize_regex" /> + </param> + <param name="title" type="text" value="Title" label="Plot title" /> + <expand macro="use_defaults_no"> + <param name="cluster" type="boolean" checked="false" label="Partitioning along the x-axis" /> + <param name="alpha" type="float" min="0" max="1" value="0.5" label="Smoothing parameter" /> + </expand> + </when> + </conditional> + </when> + </conditional> + </section> </inputs> <outputs> @@ -38,7 +93,100 @@ </outputs> <tests> - <expand macro="inspecttrajectories_tests" /> + <test> + <!-- default, generates a blank report --> + <param name="inputrds" value="out_traject_default.ltree.rdat" /> + <output name="outpdf" value="out_traject_inspect_default.ltree.blank.pdf" compare="sim_size" delta="15" /> + </test> + <test> + <!-- stemID branch inspection: vignette search "getproj" --> + <param name="inputrds" value="out_traject_default.ltree.rdat" /> + <section name="trjsid" > + <conditional name="basic" > + <param name="doit" value="yes" /> + <param name="i" value="3" /> + <param name="br" value="1,3,8" /> + </conditional> + </section> + <output name="outpdf" value="out_traject_inspect_stemid.pdf" compare="sim_size" delta="30" /> + <output name="outdiffgenes" value="out_traject_inspect_stemid.tabular" /> + </test> + <test> + <!-- fateID trajectory inspection: vignette search "cellsfromtree" --> + <param name="inputrds" value="out_traject_default.ltree.rdat" /> + <section name="trjfid" > + <conditional name="basic" > + <param name="doit" value="yes" /> + <param name="cellsfromz" value="2,1,4" /> + <expand macro="test_nondef" > + <param name="filterset_minexpr" value="2" /> + <param name="getsom_nb" value="1000" /> + <param name="getsom_alpha" value="0.5" /> + <param name="plotheat_xlab" value="true" /> + </expand> + </conditional> + </section> + <output name="outpdf" value="out_traject_inspect_fateid.pdf" compare="sim_size" delta="15" /> + </test> + <test> + <!-- fateID trajectory inspection with som: vignette search "SOM" --> + <param name="inputrds" value="out_traject_default.ltree.rdat" /> + <section name="trjfid" > + <conditional name="basic" > + <param name="doit" value="yes" /> + <param name="cellsfromz" value="2,1,4" /> + <conditional name="som" > + <param name="doit" value="yes" /> + <expand macro="test_nondef" > + <param name="typer" value="cln" /> + <param name="use_genes" value="12" /> + </expand> + </conditional> + </conditional> + </section> + <output name="outpdf" value="out_traject_inspect_fateid_som.pdf" compare="sim_size" delta="15" /> + </test> + <test> + <!-- uses all 3 sections with additional non-default params --> + <param name="inputrds" value="out_traject_default.ltree.rdat" /> + <section name="trjsid" > + <conditional name="basic" > + <param name="doit" value="yes" /> + <param name="i" value="5" /> + <param name="br" value="6,5,3" /> + <expand macro="test_nondef" > + <param name="zscore" value="true" /> + <param name="ndiffgenes" value="14" /> + </expand> + </conditional> + </section> + <section name="trjfid" > + <conditional name="basic" > + <param name="doit" value="yes" /> + <param name="cellsfromz" value="3,9,2" /> + <expand macro="test_nondef" > + <param name="filterset_minexpr" value="1" /> + <param name="filterset_minnumber" value="2" /> + <param name="procsom_minsom" value="5" /> + <param name="procsom_corthr" value="0.5" /> + <param name="plotheat_ygrid" value="true" /> + </expand> + <conditional name="som" > + <param name="doit" value="yes" /> + <conditional name="use_genes" > + <param name="typer" value="genelist" /> + <param name="use_genes" value="Clca4,Mki67,Defa24,Ybx1,Vasp,Apoa1" /> + <expand macro="test_nondef" > + <param name="cluster" value="true" /> + <param name="alpha" value="0.1" /> + </expand> + </conditional> + </conditional> + </conditional> + </section> + <output name="outdiffgenes" value="out_traject_inspect_allthree.tabular" /> + <output name="outpdf" value="out_traject_inspect_allthree.pdf" compare="sim_size" delta="15" /> + </test> </tests> <help><![CDATA[ StemID2 and FateID
--- a/scripts/cluster.R Mon Apr 15 17:53:28 2019 -0400 +++ b/scripts/cluster.R Mon Aug 12 13:03:03 2019 -0400 @@ -1,5 +1,5 @@ #!/usr/bin/env R -VERSION = "0.4" +VERSION = "0.5" args = commandArgs(trailingOnly = T) @@ -139,6 +139,19 @@ write.table(df, file=out.genelist, sep="\t", quote=F) } + +writecellassignments <- function(sc){ + dat <- sc@cluster$kpart + tab <- data.frame(row.names = NULL, + cells = names(dat), + cluster.initial = dat, + cluster.final = sc@cpart, + is.outlier = names(dat) %in% sc@out$out) + + write.table(tab, file=out.assignments, sep="\t", quote=F, row.names = F) +} + + pdf(out.pdf) if (use.filtnormconf){ @@ -148,6 +161,7 @@ message(paste(" :: ", sprintf("%.1f", 100 * nrow(getfdata(sc))/nrow(sc@expdata)), "% of genes remain,", sprintf("%.1f", 100 * ncol(getfdata(sc))/ncol(sc@expdata)), "% of cells remain")) + write.table(as.matrix(sc@ndata), file=out.table, col.names=NA, row.names=T, sep="\t", quote=F) } if (use.cluster){ @@ -161,6 +175,7 @@ sc <- do.clustmap(sc) mkgenelist(sc) + writecellassignments(sc) } dev.off()
--- a/scripts/clusterinspect.R Mon Apr 15 17:53:28 2019 -0400 +++ b/scripts/clusterinspect.R Mon Aug 12 13:03:03 2019 -0400 @@ -1,5 +1,5 @@ #!/usr/bin/env R -VERSION = "0.2" +VERSION = "0.5" args = commandArgs(trailingOnly = T) @@ -17,13 +17,30 @@ test$line = 3 do.plotting <- function(sc){ - print(plotmap(sc, final = FALSE, fr = FALSE)) + + sc.tmp <- sc + + ## If it's a subset, we need to get clever and subset specific parts + if (!(is.null(plotting.cln) || is.na(plotting.cln))){ + cellstokeep <- names(sc.tmp@cpart[sc.tmp@cpart %in% plotting.cln]) + + ## Subselect partitions for initial and final clusters + sc.tmp@cpart <- sc.tmp@cpart[cellstokeep] + sc.tmp@cluster$kpart <- sc.tmp@cluster$kpart[cellstokeep] + + ## Subselect tSNE and FR data + ## - Note: no names in tsne, so we assume it follows the ndata naming + sc.tmp@tsne <- sc.tmp@tsne[colnames(sc.tmp@ndata) %in% cellstokeep,] + sc.tmp@fr <- sc.tmp@fr[cellstokeep,] + } + + print(plotmap(sc.tmp, final = FALSE, fr = FALSE)) print(do.call(mtext, c("Initial Clustering tSNE", test))) - print(plotmap(sc, final = TRUE, fr = FALSE)) + print(plotmap(sc.tmp, final = TRUE, fr = FALSE)) print(do.call(mtext, c("Final Clustering tSNE", test))) - print(plotmap(sc, final = FALSE, fr = TRUE)) + print(plotmap(sc.tmp, final = FALSE, fr = TRUE)) print(do.call(mtext, c("Initial Clustering Fruchterman-Reingold", test))) - print(plotmap(sc, final = TRUE, fr = TRUE)) + print(plotmap(sc.tmp, final = TRUE, fr = TRUE)) print(do.call(mtext, c("Final Clustering Fruchterman-Reingold", test))) }
--- a/test-data/matrix.genelist Mon Apr 15 17:53:28 2019 -0400 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,201 +0,0 @@ -n mean.ncl mean.cl fc pv padj -ENSDARG00000108798 1 0.13225291542414 22.6067842060994 170.935998904815 5.63439752001819e-20 3.43529216795509e-16 -ENSDARG00000109030 1 0.186711263967539 34.6842724288799 185.764220603798 2.48233981248329e-19 7.56741291835532e-16 -ENSDARG00000108379 1 0.168199368419563 24.3874482842581 144.991319012715 4.11994031997326e-19 8.37309204362565e-16 -ENSDARG00000107599 1 0.194980763882093 20.9729826342942 107.564367975175 3.56470297051077e-18 5.43349850280104e-15 -ENSDARG00000107321 1 0.1 1.94683879586459 19.4683879586459 1.39905888784068e-09 1.70601240783292e-06 -ENSDARG00000105564 1 0.1 1.72137476189456 17.2137476189456 5.11677986670882e-09 5.19950114122061e-06 -ENSDARG00000090447 1 0.1 1.4302236712519 14.302236712519 3.51620742912878e-08 3.06261667077116e-05 -ENSDARG00000079591 1 0.11578454541401 1.47825950354764 12.767330028907 4.83701028926633e-08 3.6864064667071e-05 -ENSDARG00000107138 1 0.118643244515052 1.63372529702772 13.7700659123533 8.59237265783686e-08 5.36793599520202e-05 -ENSDARG00000041060 1 0.11578454541401 1.4382199783801 12.421519411226 8.84781911763433e-08 5.36793599520202e-05 -ENSDARG00000077004 1 0.1 1.29299827494356 12.9299827494356 1.0364347799037e-07 5.36793599520202e-05 -ENSDARG00000100003 1 0.1 1.29775216728225 12.9775216728225 1.05650700249999e-07 5.36793599520202e-05 -ENSDARG00000102310 1 0.129577145917004 1.72085266498654 13.2805260743185 1.20446928962954e-07 5.64896096836256e-05 -ENSDARG00000068650 1 0.115601615820105 1.27001296748153 10.9861177845289 2.26895176903151e-07 9.88128495413223e-05 -ENSDARG00000086678 1 0.128030253978187 1.35647794130946 10.5949796955067 3.16910756640577e-07 0.000122948041506908 -ENSDARG00000105857 1 0.105900327795989 1.37647789967886 12.9978625026598 3.22645344285801e-07 0.000122948041506908 -ENSDARG00000105159 1 0.1 1.16700756453584 11.6700756453584 3.44329675018825e-07 0.000123492825211163 -ENSDARG00000102997 1 0.125930886410326 1.40950930054417 11.1927211879659 3.75314178672153e-07 0.00012712725263134 -ENSDARG00000052615 1 0.109701288024116 1.2495211865536 11.3902143635625 4.65728278602572e-07 0.000149449753402099 -ENSDARG00000103046 1 0.111800655591977 1.27916812008568 11.441508221151 5.24275745356588e-07 0.000158237500570989 -ENSDARG00000105009 1 0.107892272707005 1.29390774186871 11.9925895470056 5.45020093815115e-07 0.000158237500570989 -ENSDARG00000069937 1 0.1 1.08683574177416 10.8683574177416 5.98190849813324e-07 0.000165780436877811 -ENSDARG00000034941 1 0.1 1.14009483992511 11.4009483992511 7.60152472832673e-07 0.000201506505515687 -ENSDARG00000099753 1 0.121501943616094 1.29493126702895 10.6577000210014 1.17807930734258e-06 0.000299281230702821 -ENSDARG00000076005 1 0.1 1.05371686166029 10.5371686166029 1.32743044013891e-06 0.00029931174491888 -ENSDARG00000041750 2 0.118845246250015 4.25427279350112 35.7967434772393 1.39550360463363e-07 0.000850838547745124 -ENSDARG00000103735 2 0.119337350478008 3.54379562135099 29.6956117020886 9.49083511390771e-07 0.00228722429441009 -ENSDARG00000101886 2 0.215092571957637 5.30165157388469 24.6482318084367 1.12541789129576e-06 0.00228722429441009 -ENSDARG00000090914 2 0.122178374276303 3.1027845112115 25.3955295246819 1.65948635889896e-06 0.00252947208255174 -ENSDARG00000090496 2 0.111721506487875 2.96053053983273 26.4991999562235 2.19562316256802e-06 0.00267734288443544 -ENSDARG00000078308 2 0.1 2.64934229048601 26.4934229048601 2.65009460603265e-06 0.00269293780216352 -ENSDARG00000075253 2 0.121220276605515 2.61038833734525 21.534254915456 7.21609407341446e-06 0.00624434487572006 -ENSDARG00000086877 2 0.174901657246515 3.44134351834849 19.6758771330457 8.19333426369698e-06 0.00624434487572006 -ENSDARG00000040487 2 0.132551316615554 2.19579056306671 16.5655884764638 5.8552609335758e-05 0.0389789127774957 -ENSDARG00000052912 2 0.29374667394921 3.98299817184644 13.5592962408661 6.97030598265265e-05 0.0389789127774957 -ENSDARG00000052558 2 0.307266492503596 4.25427279350112 13.8455474231422 7.0324428498024e-05 0.0389789127774957 -ENSDARG00000012071 2 0.137312605449286 2.19567880375095 15.9903658995232 8.47157718952949e-05 0.0419371525043834 -ENSDARG00000052515 2 0.142741634362038 2.26259398496241 15.8509743500887 8.94182356170222e-05 0.0419371525043834 -ENSDARG00000002897 2 0.155152648803008 2.06464107607378 13.3071597037005 0.000108713544348236 0.047344748563657 -ENSDARG00000087601 2 0.220544850783226 3.03609059600962 13.7663182124972 0.000128704633773863 0.0523141434746162 -ENSDARG00000020581 2 0.132645565796613 1.97090942840696 14.8584644844365 0.000140065988406919 0.0533738957073114 -ENSDARG00000005368 2 0.109588039327212 1.71813717094559 15.6781450009842 0.000165794148578122 0.0594615837576948 -ENSDARG00000013006 2 0.139223723099655 1.9643911375366 14.1096006758166 0.000196735880595178 0.0666388146660444 -ENSDARG00000059690 2 0.220277183931972 2.42251643631602 10.9975821965482 0.000320441653294837 0.102828040007296 -ENSDARG00000007906 2 0.161941295286845 1.92988721804511 11.9172025555725 0.000357508658491656 0.108986514541181 -ENSDARG00000099235 2 0.355572431070773 3.47503570154998 9.77307405719064 0.0005137976966319 0.145982086583318 -ENSDARG00000058248 2 0.187973361461239 1.79476640836755 9.54798272699739 0.000556896211693426 0.145982086583318 -ENSDARG00000013168 2 0.122745621089648 1.36391188645772 11.1116948559948 0.000584796437234187 0.145982086583318 -ENSDARG00000089066 2 0.213734450098499 2.20785246938186 10.3298858390137 0.000627746814730416 0.145982086583318 -ENSDARG00000018559 2 0.216649831175799 2.21719732132153 10.2340136121426 0.000632257513971329 0.145982086583318 -ENSDARG00000108459 3 0.1 40.0486298308386 400.486298308386 1.17120842343972e-21 7.140857757712e-18 -ENSDARG00000106631 3 0.1 16.7209600896533 167.209600896533 1.27082371024938e-20 2.65261585580336e-17 -ENSDARG00000106893 3 0.109440335656379 20.3249081306309 185.716792704722 1.30520708010662e-20 2.65261585580336e-17 -ENSDARG00000106328 3 0.109440335656379 13.3207925346078 121.717394731248 1.812771473855e-19 2.76311691902349e-16 -ENSDARG00000107500 3 0.194403356563789 23.5777850278158 121.282808304183 8.16923386605951e-19 9.96156377627297e-16 -ENSDARG00000107098 3 0.407548050942585 20.6628907529452 50.7005019534645 7.91289296843658e-15 8.0408180714263e-12 -ENSDARG00000106626 3 0.322331939034494 11.7938446379915 36.5891281928763 3.19108227960682e-14 2.77943266553754e-11 -ENSDARG00000063254 3 0.1 1.54154831743065 15.4154831743065 6.16683050966547e-09 4.69989570217879e-06 -ENSDARG00000105532 3 0.1 1.2691654791726 12.691654791726 6.59753120473869e-08 4.46946086169909e-05 -ENSDARG00000102291 3 0.1 1.20003533117234 12.0003533117234 1.03313712547645e-07 6.29903705402992e-05 -ENSDARG00000098646 3 0.1 1.22194683967797 12.2194683967797 1.14432499391301e-07 6.3426813526251e-05 -ENSDARG00000096706 3 0.1 1.2471262087964 12.471262087964 1.27992769850814e-07 6.50309931483676e-05 -ENSDARG00000056322 3 0.1 1.15253315794626 11.5253315794626 1.81837361851978e-07 7.7399961295303e-05 -ENSDARG00000087417 3 0.1 1.1540362057603 11.540362057603 1.83177220114043e-07 7.7399961295303e-05 -ENSDARG00000100789 3 0.109440335656379 1.31018291186832 11.9716638660726 1.90421423557413e-07 7.7399961295303e-05 -ENSDARG00000076768 3 0.1 1.18681657849799 11.8681657849799 2.13741147520963e-07 8.1448736027207e-05 -ENSDARG00000091869 3 0.1 1.10007995897679 11.0007995897679 3.18290687070941e-07 0.000107899200006755 -ENSDARG00000101145 3 0.1 1.10059661347904 11.0059661347904 3.19125599577782e-07 0.000107899200006755 -ENSDARG00000091967 3 0.1 1.1236667156373 11.236667156373 3.57604258883062e-07 0.000107899200006755 -ENSDARG00000099207 3 0.1 1.13395237668583 11.3395237668583 3.75510462080806e-07 0.000107899200006755 -ENSDARG00000079762 3 0.1 1.13482318163712 11.3482318163712 3.77047564405683e-07 0.000107899200006755 -ENSDARG00000100342 3 0.118880671312758 1.11971099812246 9.4187809149955 3.89336132548567e-07 0.000107899200006755 -ENSDARG00000091631 3 0.1 1.0336438761938 10.336438761938 5.26923351033268e-07 0.000136411029348096 -ENSDARG00000101776 3 0.1 1.03711906765213 10.3711906765213 5.36963212129621e-07 0.000136411029348096 -ENSDARG00000077068 3 0.1 1.0531104586081 10.531104586081 5.84457768774286e-07 0.000138268351040062 -ENSDARG00000103380 4 0.118732364944267 4.29504915730337 36.1742070859909 1.87281976492061e-10 1.1418582106721e-06 -ENSDARG00000037116 4 0.105401803653084 3.3624106287092 31.9008832123606 6.24221125271845e-10 1.90293810039122e-06 -ENSDARG00000102060 4 0.65335417049578 14.8778451385581 22.7714856817527 8.75163683692952e-09 1.77862432649198e-05 -ENSDARG00000102188 4 0.855834419578633 14.179562474169 16.5681142868153 2.26488443176202e-07 0.000318621197056678 -ENSDARG00000105127 4 1.2186946752299 22.1380023868484 18.1653393887785 2.6129342058117e-07 0.000318621197056678 -ENSDARG00000101413 4 0.163690456191388 2.80861987923491 17.1581162676402 3.73531302510856e-07 0.000379570058568115 -ENSDARG00000104118 4 0.989700274307544 15.0174505397369 15.1737358567916 6.77423754716057e-07 0.000590036090357685 -ENSDARG00000075333 4 0.171748726689293 2.3373595505618 13.6091812476157 1.50828806455175e-06 0.00110856685474393 -ENSDARG00000071113 4 0.128246671314964 1.77801966292135 13.8640609123855 1.68613897165096e-06 0.00110856685474393 -ENSDARG00000098655 4 0.153148025950228 2.3373595505618 15.2620938863516 1.81821691773647e-06 0.00110856685474393 -ENSDARG00000105209 4 1.54534561735873 23.4939455347753 15.2030363116638 2.06788281463169e-06 0.00114617104734631 -ENSDARG00000092124 4 0.167783479025787 2.20080865760876 13.1169568683846 2.51104682731582e-06 0.00117801683762716 -ENSDARG00000102424 4 0.995458642689216 12.8896835165257 12.9484872236424 2.51176298001526e-06 0.00117801683762716 -ENSDARG00000043130 4 0.151115417204371 1.96446629213483 12.9997741360702 2.99111460631241e-06 0.00128615842311759 -ENSDARG00000102305 4 1.00227191359572 12.5079159542882 12.4795634643849 3.16424083102573e-06 0.00128615842311759 -ENSDARG00000104251 4 1.25657639483624 16.1896340694694 12.8839234415026 4.03444006517729e-06 0.00153737381733662 -ENSDARG00000016348 4 0.130718858025808 1.82528654590158 13.9634523546802 7.40568226139146e-06 0.00265602616162963 -ENSDARG00000107668 4 2.0525826270598 26.1651913798814 12.7474485240877 1.23716885008222e-05 0.00419056582163961 -ENSDARG00000102171 4 0.209616501874572 2.24413623595506 10.7059139709233 1.41333189895541e-05 0.0045353076778585 -ENSDARG00000102300 4 1.16111393148624 12.0981196442604 10.4194078773774 1.50789368610256e-05 0.00459681390208365 -ENSDARG00000096457 4 0.122651614608118 1.49305555555556 12.1731422804827 1.58593667087189e-05 0.00460450280109806 -ENSDARG00000102445 4 1.3115800466299 12.8389700502745 9.78893364782744 3.22055454104447e-05 0.00892532774397644 -ENSDARG00000089213 4 0.12779725761052 1.38681171572623 10.8516547354461 3.54676572069977e-05 0.00940201330395934 -ENSDARG00000105147 4 1.98689642082857 20.3540415496165 10.2441382128659 4.3971555566058e-05 0.0111706072619273 -ENSDARG00000055075 4 0.169741281332608 1.59157303370787 9.37646411770147 4.82148454826373e-05 0.0117586365163056 -ENSDARG00000102705 5 0.209972550128478 13.5265663218744 64.4206412390466 5.51324807068058e-17 3.36142734869395e-13 -ENSDARG00000102125 5 0.233431470953978 10.4713938457086 44.8585351534415 1.77193932024865e-15 5.401757017778e-12 -ENSDARG00000105858 5 0.462925527809972 31.4875823414393 68.0186778430692 3.07338903352449e-15 6.24615097913294e-12 -ENSDARG00000104031 5 0.417738843341455 16.1643635536032 38.6949018777042 4.63887116279502e-14 7.07079936989031e-11 -ENSDARG00000103505 5 0.455814556718891 13.0057940562389 28.5330818521002 8.00401069394782e-13 8.56538424116237e-10 -ENSDARG00000107953 5 0.768435533248605 31.7530045869607 41.3216245385258 8.42911357175237e-13 8.56538424116237e-10 -ENSDARG00000104624 5 0.586363931713224 18.662951819579 31.8282738930583 1.03929690687218e-12 9.05227605885665e-10 -ENSDARG00000107433 5 0.860111367092678 34.1815671825659 39.7408620445342 2.15526546944614e-12 1.64258169590164e-09 -ENSDARG00000105003 5 0.659817827005108 18.3691699000608 27.8397599280362 4.05465742511909e-12 2.74680514677235e-09 -ENSDARG00000102147 5 0.498991534801074 11.6870331272134 23.4213053972398 5.94420491301204e-12 3.62418173546344e-09 -ENSDARG00000104149 5 0.591410894158957 12.7563563319083 21.5693631245145 2.16849130014449e-11 1.20193558699827e-08 -ENSDARG00000105370 5 0.845960751451064 21.1606791594491 25.0137835864755 2.88314145770954e-11 1.46487612230459e-08 -ENSDARG00000102471 5 0.591016786322552 11.9787944006235 20.2681119687961 3.97757415291587e-11 1.86548227771754e-08 -ENSDARG00000104225 5 0.677500552858849 13.028628755702 19.2304326553315 8.48219367396283e-11 3.69399534501081e-08 -ENSDARG00000102765 5 0.710420955035912 7.69305581025659 10.8288694973358 2.636812754442e-08 1.07177649092219e-05 -ENSDARG00000103900 5 0.926463189872816 8.83493365876684 9.53619502139064 1.4891189989945e-07 5.67447408554341e-05 -ENSDARG00000103207 5 0.809155048654974 6.95078624248591 8.5901784262978 3.17107715066467e-07 0.000113729749338838 -ENSDARG00000100850 5 0.10266571617904 1.06794343733293 10.4021427705286 3.65574205111742e-07 0.000123828107142572 -ENSDARG00000100222 5 0.117813108572049 1.02376729835175 8.6897571141302 1.53996760366899e-06 0.000376422800056221 -ENSDARG00000078327 5 0.105063152860082 0.901542752536163 8.58096038424425 1.59754149266682e-06 0.000376422800056221 -ENSDARG00000076659 5 0.107997148537119 0.902203678729055 8.35395833084392 1.63475485205548e-06 0.000376422800056221 -ENSDARG00000078658 5 0.103730294770562 0.912796947051807 8.79971419218265 1.70133206806343e-06 0.000376422800056221 -ENSDARG00000056152 5 0.109363711084075 0.821293559013241 7.50974478528675 2.30482270129059e-06 0.000411281146237518 -ENSDARG00000042974 5 0.111761147765118 0.823872029464418 7.37172126395753 2.39063753034635e-06 0.000411281146237518 -ENSDARG00000071240 5 0.110394585218162 0.833280258058883 7.54819864046915 2.53282560651987e-06 0.000411281146237518 -ENSDARG00000094010 6 0.111452978443105 4.44006630603854 39.8380228869803 6.38525950158321e-09 3.89309271811528e-05 -ENSDARG00000103594 6 1.13093354481739 15.5296459037369 13.7317050810836 2.76457859470178e-05 0.0842781784594839 -ENSDARG00000062868 6 0.148310913360651 2.35945460562356 15.9088400992179 4.74734577891717e-05 0.0964818907135266 -ENSDARG00000101986 6 1.06278096568171 11.921270245687 11.2170528365083 0.000104992972874659 0.160035538904199 -ENSDARG00000054302 6 0.114167149192999 1.69385093295239 14.8365878006548 0.000159146295603994 0.194062992859511 -ENSDARG00000004931 6 0.146147835707832 1.50738986001055 10.3141442547532 0.000273967771010734 0.271604620436783 -ENSDARG00000099280 6 0.174438502158969 1.76393397177854 10.1120678631546 0.000341416753578639 0.271604620436783 -ENSDARG00000104949 6 1.84012138595949 18.1422310489148 9.85925775731088 0.000356378048793549 0.271604620436783 -ENSDARG00000062192 6 0.219250133212514 1.93124420721189 8.8084060835666 0.000490459989622772 0.332259395192227 -ENSDARG00000014366 6 0.135611988044807 1.00529217065682 7.41300371118129 0.00101799998551309 0.595386680453464 -ENSDARG00000099776 6 0.165584080167415 1.27670227417667 7.7102960193145 0.00113607374054731 0.595386680453464 -ENSDARG00000014017 6 0.130935742022126 1.13118855425646 8.63926485455216 0.00124672027506227 0.595386680453464 -ENSDARG00000004527 6 0.161768890546368 1.40411105178857 8.67973469463903 0.00130573477419644 0.595386680453464 -ENSDARG00000015851 6 0.211160123828113 1.42542261170536 6.75043462687903 0.00136713359461186 0.595386680453464 -ENSDARG00000017037 6 0.11090453956233 0.860975851160853 7.76321559566973 0.0016197152098211 0.658360242285282 -ENSDARG00000018146 6 0.124837526779085 0.822444876257476 6.58812215747347 0.00208359331018621 0.755769031195207 -ENSDARG00000094732 6 0.201066488547157 1.22521720805992 6.09359230826055 0.00228107117057515 0.755769031195207 -ENSDARG00000021855 6 0.109168767482819 1.07499270061704 9.84707188149048 0.00231989733905826 0.755769031195207 -ENSDARG00000062058 6 0.122552380041629 0.881708382142295 7.19454311570935 0.00235519297895833 0.755769031195207 -ENSDARG00000089032 6 0.150988259365316 1.23158448554809 8.15682285977133 0.00250441767427358 0.7634717280023 -ENSDARG000001032071 6 1.66376642491526 10.8791567395556 6.53887263057956 0.00264170112721662 0.766973893935226 -ENSDARG00000101652 6 0.117661697666601 1.00798229432208 8.5667835354393 0.00288382859301584 0.79921376961898 -ENSDARG00000036809 6 0.178132306667668 1.05062387043439 5.89799733742001 0.0033324167246799 0.875163147449564 -ENSDARG00000037018 6 0.13196794390221 0.929979569662069 7.04701113136379 0.00344495908459727 0.875163147449564 -ENSDARG00000077737 6 0.229953264257691 1.41304633493481 6.1449283596658 0.00500781997902645 1 -ENSDARG00000104242 7 0.447206731128574 14.3756294075549 32.1453779805067 2.30113131457065e-13 1.40299976249373e-09 -ENSDARG00000102056 7 0.463880684093747 12.1589657830654 26.2114078037536 2.08698751307313e-12 6.36218143360344e-09 -ENSDARG00000102470 7 0.825035359668146 13.4793918534294 16.3379565438882 1.10737149331994e-09 2.25054799825723e-06 -ENSDARG00000104836 7 1.0201341649907 17.1513950198419 16.8128816860067 1.79596375972667e-09 2.73749776076337e-06 -ENSDARG00000104197 7 0.868813330133413 13.2077072533743 15.2020080669644 2.27590161725937e-09 2.77523443208607e-06 -ENSDARG00000102167 7 0.911623267267556 12.5281341487954 13.7426660755889 6.03157787269791e-09 6.12908838163986e-06 -ENSDARG00000105195 7 1.35837446108669 20.6382812121625 15.193366632977 1.20671859497662e-08 1.05105189622463e-05 -ENSDARG00000105669 7 1.55320812115904 19.4772872937536 12.5400369908053 8.3653757257818e-08 6.37546197501146e-05 -ENSDARG00000102820 7 1.12752840077096 11.8419301998565 10.5025560258699 1.31966552783421e-07 8.94000080356132e-05 -ENSDARG00000092989 7 0.133042434280628 1.21637077531351 9.14272789648302 2.82608004083509e-07 0.000172306100089716 -ENSDARG00000061806 7 0.13278988077121 1.12194957012091 8.44905924762421 7.25529240779189e-07 0.000402141071002793 -ENSDARG00000095512 7 0.134733887854507 1.08703031970933 8.06798005326739 1.27119959332606e-06 0.000645875326709081 -ENSDARG00000090634 7 0.14895440642195 1.13152352961684 7.5964421382173 1.67131715307221e-06 0.000689864250270158 -ENSDARG00000087555 7 0.106854974051793 0.876076138239269 8.19873989033614 1.8103703467808e-06 0.000689864250270158 -ENSDARG00000069269 7 0.131991994835001 1.03283452343888 7.82497851274993 2.01967640768832e-06 0.000724351003392688 -ENSDARG00000059558 7 0.112699471961133 0.967793901919788 8.58738630340302 3.00038333622434e-06 0.000888688802910736 -ENSDARG00000060566 7 0.127846770606366 1.01035836252462 7.90288528785342 3.78083155119976e-06 0.000992894397782283 -ENSDARG00000100206 7 0.139832328948304 1.04998427596811 7.5088806992286 4.54325865418129e-06 0.000992894397782283 -ENSDARG00000057206 7 0.121885624201679 0.88338782626096 7.24767856789453 4.55979057534917e-06 0.000992894397782283 -ENSDARG00000070623 7 0.12670584582223 0.943766658477006 7.44848552450473 5.87432370429556e-06 0.00110819712277626 -ENSDARG00000100514 7 0.121454028840368 0.820406006163069 6.75486860334098 7.45177975757828e-06 0.00116897636495374 -ENSDARG00000039265 7 0.131635113592875 0.87126351316452 6.6187773868543 9.18440046329954e-06 0.00133326880058898 -ENSDARG00000095904 7 0.168129332968718 1.08581712830248 6.45822539785194 1.03525068635208e-05 0.00139894259529663 -ENSDARG00000039350 7 0.133011711683714 0.897490439418404 6.74745425088979 1.06942592337219e-05 0.00139894259529663 -ENSDARG00000070512 7 0.12670239660567 0.910072943888678 7.18276029711621 1.10515908767776e-05 0.00139894259529663 -ENSDARG00000103551 8 1.61835767565622 32.2654393890171 19.93714978732 1.97522608849965e-06 0.0120429534615824 -ENSDARG00000078383 8 0.182610070366458 2.60621132144157 14.2720021749703 1.6366673021414e-05 0.0447464228757837 -ENSDARG00000093799 8 0.115321017732354 1.62941087246693 14.1293487042283 2.20172656433248e-05 0.0447464228757837 -ENSDARG00000036625 8 0.119337255991682 1.48171990557204 12.4162391137543 8.39663841121717e-05 0.127985760982978 -ENSDARG00000016706 8 0.180025642511508 1.82486224560436 10.1366795315712 0.000126325160000158 0.154040900104192 -ENSDARG00000040251 8 0.121540951107174 1.29918699186992 10.6892942669529 0.000198034398721646 0.201235954834313 -ENSDARG00000098674 8 0.253812764860557 2.21906493621457 8.74292093793513 0.000293366486612692 0.255522209839655 -ENSDARG00000098856 8 0.194513509828675 1.64663383642527 8.46539573459756 0.000451769293400039 0.296721325111977 -ENSDARG00000062226 8 0.144813536626749 1.38335068305012 9.55263378875725 0.000474939776632199 0.296721325111977 -ENSDARG00000075123 8 0.142439266402834 1.40532056857682 9.86610366696508 0.00048666774661633 0.296721325111977 -ENSDARG00000070239 8 0.154515334275144 1.27289156626506 8.23796273836765 0.000570475520527726 0.316199022605232 -ENSDARG00000026180 8 0.151758427938041 1.41564289946672 9.32826544595394 0.000683459166751383 0.347254211640265 -ENSDARG00000104786 8 0.182483274181232 1.61787042996763 8.8658559927024 0.00078762411588287 0.369395710349066 -ENSDARG00000101061 8 0.255483815147335 1.62613690189592 6.36493118344166 0.00150635451453941 0.61883047685698 -ENSDARG00000016584 8 0.135272462199637 0.880366460381022 6.50809814551733 0.00158197684628705 0.61883047685698 -ENSDARG00000070256 8 0.203213690445973 1.26976802854965 6.24843742448165 0.00178015471505539 0.61883047685698 -ENSDARG00000062745 8 0.130428963634169 0.955091918456498 7.32269805604966 0.00182445112094621 0.61883047685698 -ENSDARG00000086037 8 0.19632007891506 1.29918699186992 6.61769799120764 0.00182695564760138 0.61883047685698 -ENSDARG00000027310 8 0.193675103177481 1.44901164174846 7.48166190685131 0.00211840042255122 0.678671506770642 -ENSDARG00000103109 8 0.14293938597264 0.936026560054929 6.54841598545899 0.00240745117064372 0.678671506770642 -ENSDARG00000078352 8 0.177346649454539 1.15788746142844 6.52895030715114 0.0024876207141355 0.678671506770642 -ENSDARG00000035389 8 0.143574445102194 1.00421593036345 6.99439186164828 0.0026995869399452 0.678671506770642 -ENSDARG00000103106 8 0.175526172559138 1.23041862190788 7.00988692437495 0.0027017171947583 0.678671506770642 -ENSDARG00000007257 8 0.157933075856061 0.866671316548008 5.48758587680447 0.0028384639169755 0.678671506770642 -ENSDARG00000010596 8 0.110516369140454 0.836156159426777 7.56590327686312 0.00297745919336174 0.678671506770642