Repository 'baredsc_1d'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/baredsc_1d

Changeset 0:02b1fe7aed76 (2023-10-02)
Next changeset 1:4fff5a293013 (2023-10-16)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/baredsc commit dad42d788f7407976a67ed50cc886ebe79740e32
added:
baredsc_1d.xml
macros.xml
test-data/2d_small_1gauss.npz
test-data/2d_small_1gauss.png
test-data/2d_small_2gauss.npz
test-data/2d_small_2gauss.png
test-data/2d_small_2gauss_convergence.png
test-data/2d_small_2gauss_corr.txt
test-data/2d_small_2gauss_individuals.png
test-data/2d_small_2gauss_median.png
test-data/2d_small_2gauss_neff.txt
test-data/2d_small_2gauss_p.png
test-data/2d_small_2gauss_p.txt
test-data/2d_small_2gauss_pdf2d.txt
test-data/2d_small_2gauss_pdf2d_flat.txt
test-data/2d_test2.pdf
test-data/2d_test2_split-10.5.pdf
test-data/2d_test2_split-10.5.txt
test-data/2d_test2_split-10.5_renorm.pdf
test-data/2d_test2_split-9.5.pdf
test-data/2d_test2_split-9.5.txt
test-data/2d_test2_split-9.5_renorm.pdf
test-data/generate_small_inputs_for_combine.sh
test-data/nih3t3_generated_2d_2.txt
test-data/small_1gauss.npz
test-data/small_1gauss.png
test-data/small_2gauss.npz
test-data/small_2gauss.png
test-data/small_2gauss_convergence.png
test-data/small_2gauss_corner.png
test-data/small_2gauss_individuals.png
test-data/small_2gauss_means.txt.gz
test-data/small_2gauss_neff.txt
test-data/small_2gauss_p.png
test-data/small_2gauss_p.txt
test-data/small_2gauss_pdf.txt
test-data/small_2gauss_posterior_andco.png
test-data/small_2gauss_posterior_individuals.png
test-data/small_2gauss_posterior_per_cell.txt
test-data/small_2gauss_with_posterior.png
test-data/test1_res.png
test-data/test2_res.pdf
b
diff -r 000000000000 -r 02b1fe7aed76 baredsc_1d.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/baredsc_1d.xml Mon Oct 02 13:24:52 2023 +0000
[
b'@@ -0,0 +1,303 @@\n+<tool id="baredsc_1d" name="baredSC 1d" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="20.01" license="MIT">\n+    <description>Compute distribution for a single gene</description>\n+    <macros>\n+        <import>macros.xml</import>\n+    </macros>\n+    <expand macro="edam_topics"/>\n+    <edam_operations>\n+        <edam_operation>operation_2495</edam_operation>\n+    </edam_operations>\n+    <xrefs>\n+        <xref type="bio.tools">baredsc</xref>\n+    </xrefs>\n+    <expand macro="requirements"/>\n+    <command detect_errors="exit_code"><![CDATA[\n+    baredSC_1d\n+    ## Required inputs:\n+    @REQUIRED_INPUTS_1D@\n+\n+    ## Filter cells\n+    @FILTER_CELLS@\n+    ## MCMC\n+    @MCMC_1D@\n+    @BAREDSC_COMMON@\n+    ## Plots\n+    @PLOTS@\n+    @PRETTYBINS_1D@\n+    ## Advanced\n+    @ADVANCED_COMMON_X@\n+    @ADVANCED_BAREDSC_COMMON@\n+    ## Outputs\n+    --output output\n+    --figure baredSC.$plots.image_file_format\n+    --logevidence logevidence.txt &&\n+    mkdir QC output &&\n+    mv baredSC_convergence.* QC &&\n+    mv baredSC_p.$plots.image_file_format QC &&\n+    mv baredSC_corner.* QC &&\n+    mv baredSC_neff.txt output &&\n+    @ORDER_OUTPUTS_1D@\n+    ]]></command>\n+    <inputs>\n+        <expand macro="macro_input_counts"/>\n+        <expand macro="macro_single_gene"/>\n+        <expand macro="macro_filter_cells"/>\n+        <section name="MCMC" title="MCMC parameters">\n+            <expand macro="macro_MCMC_params_common_axis" axis="x" />\n+            <expand macro="macro_scale_seed"/>\n+            <expand macro="macro_MCMC_common_baredSC"/>\n+        </section>\n+        <section name="plots" title="Customize plots">\n+            <expand macro="macro_plots"/>\n+            <expand macro="macro_prettybins_1d"/>\n+        </section>\n+        <section name="advanced" title="Advanced parameters" expanded="false">\n+            <expand macro="macro_advanced_common_axis" axis="x"/>\n+            <expand macro="macro_advanced_evidence"/>\n+            <expand macro="macro_advanced_common_baredSC"/>\n+        </section>\n+    </inputs>\n+    <outputs>\n+        <data format="npz" name="output" label="${tool.name} on $getVar(\'input_counts.input\', \'input_counts.inputAnnData\') for ${geneColName}: Numpy archive" from_work_dir="output.npz"/>\n+        <data format="txt" name="neff" label="${tool.name} on $getVar(\'input_counts.input\', \'input_counts.inputAnnData\') for ${geneColName}: Neff" from_work_dir="output/baredSC_neff.txt"/>\n+        <collection type="list" name="qc_plots" label="${tool.name} on $getVar(\'input_counts.input\', \'input_counts.inputAnnData\') for ${geneColName}: QC plots">\n+            <discover_datasets pattern="baredSC_(?P&lt;identifier_0&gt;\\S+)\\.(?P&lt;ext&gt;.*)" directory="QC"/>\n+        </collection>\n+        <data format="tabular" name="pdf" label="${tool.name} on $getVar(\'input_counts.input\', \'input_counts.inputAnnData\') for ${geneColName}: PDF with error bar" from_work_dir="output/baredSC_pdf.txt">\n+            <actions>\n+                <action name="column_names" type="metadata" default="x,low,mean,high,median" />\n+            </actions>\n+        </data>\n+        <data format="png" name="plot" label="${tool.name} on $getVar(\'input_counts.input\', \'input_counts.inputAnnData\') for ${geneColName}: plot" from_work_dir="baredSC">\n+            <change_format>\n+                <when input="plots.image_file_format" value="png" format="png" />\n+                <when input="plots.image_file_format" value="svg" format="svg" />\n+                <when input="plots.image_file_format" value="pdf" format="pdf" />\n+            </change_format>\n+        </data>\n+        <collection type="list" name="other_outputs" label="${tool.name} on $getVar(\'input_counts.input\', \'input_counts.inputAnnData\') for ${geneColName}: plots and txt">\n+            <discover_datasets pattern="baredSC_(?P&lt;identifier_0&gt;\\S+)\\.(?P&lt;ext&gt;.*)"/>\n+        </collection>\n+        <data format="txt" name="logevidence" label="${tool.name} on $getVar(\'input'..b'ntents>\n+            </output>\n+            <output name="neff">\n+                <assert_contents>\n+                    <has_text_matching expression="^317[56]\\."/>\n+                </assert_contents>\n+            </output>\n+            <output_collection name="qc_plots" count="3" type="list">\n+                <element name="convergence" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="25976" delta="2000"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="p" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="479327" delta="20000"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="corner" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="223275" delta="20000"/>\n+                    </assert_contents>\n+                </element>\n+            </output_collection>\n+            <!--In the next output pdf means probability density function-->\n+            <output name="pdf" ftype="tabular">\n+                <assert_contents>\n+                    <has_text_matching expression="x\\s+low\\s+mean\\s+high\\s+median"/>\n+                    <has_n_lines n="201"/>\n+                    <has_line_matching expression="-8\\.5\\s+0\\.4[0-9]+\\s+0\\.5[0-9]+\\s+0\\.6[0-9]+\\s+0\\.5[0-9]+"/>\n+                </assert_contents>\n+            </output>\n+            <output name="plot" ftype="pdf">\n+                <assert_contents>\n+                    <has_size value="111517" delta="10000"/>\n+                </assert_contents>\n+            </output>\n+            <output_collection name="other_outputs" count="7" type="list">\n+                <element name="individuals" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="95638" delta="10000"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="p" ftype="txt">\n+                    <assert_contents>\n+                       <has_line_matching expression="name\\s+low\\s+median\\s+high"/>\n+                       <has_line_matching expression="mu[01]\\s+-8\\.8[0-9]+\\s+-8\\.7[0-9]+\\s+-8\\.6[0-9]+"/>\n+                       <has_text_matching expression="scale[01]\\s+0\\.1[0-9]+\\s+0\\.2[0-9]+\\s+0\\.2[0-9]+" />\n+                       <has_n_lines n="6"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="with_posterior" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="112356" delta="10000"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="posterior_per_cell" ftype="txt">\n+                    <assert_contents>\n+                        <has_n_lines n="1164"/>\n+                        <has_line_matching expression="mu\\s+sd"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="posterior_andco" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="66922" delta="6000"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="posterior_individuals" ftype="pdf">\n+                    <assert_contents>\n+                        <has_size value="77887" delta="7000"/>\n+                    </assert_contents>\n+                </element>\n+                <element name="means" ftype="txt">\n+                    <assert_contents>\n+                        <has_n_lines n="100000"/>\n+                    </assert_contents>\n+                </element>\n+            </output_collection>\n+            <output name="logevidence">\n+                <assert_contents>\n+                    <has_text_matching expression="^-20[78][0-9]\\."/>\n+                </assert_contents>\n+            </output>\n+        </test>\n+    </tests>\n+    <expand macro="helpcitations"/>\n+</tool>\n'
b
diff -r 000000000000 -r 02b1fe7aed76 macros.xml
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/macros.xml Mon Oct 02 13:24:52 2023 +0000
[
b'@@ -0,0 +1,313 @@\n+<macros>\n+    <token name="@TOOL_VERSION@">1.1.1</token>\n+    <token name="@VERSION_SUFFIX@">0</token>\n+    <xml name="requirements">\n+        <requirements>\n+            <requirement type="package" version="@TOOL_VERSION@">baredsc</requirement>\n+            <requirement type="package" version="1.12">gzip</requirement>\n+        </requirements>\n+        <version_command><![CDATA[baredSC_1d --version]]></version_command>\n+    </xml>\n+    <xml name="helpcitations">\n+        <help><![CDATA[\n+\n+            .. class:: infomark\n+            \n+            **BARED (Bayesian Approach to Retreive Expression Distribution of) Single Cell**\n+            \n+            baredSC is a tool that uses a Monte-Carlo Markov Chain to estimate a confidence interval on the probability density function (PDF) of expression of one or two genes from single-cell RNA-seq data. It uses the raw counts and the total number of UMI for each cell. The PDF is approximated by a number of 1d or 2d gaussians provided by the user. The likelihood is estimated using the asumption that the raw counts follow a Poisson distribution of parameter equal to the proportion of mRNA for the gene in the cell multiplied by the total number of UMI identified in this cell.\n+            \n+            To get a description of outputs, please read the `Documentation <https://baredsc.readthedocs.io/en/latest/index.html>`_\n+\n+            This is a description of the figure with the results. \n+\n+            - When the 1d version is used, it displays the mean PDF in solid red line, the median in black dashed lines (/!\\backslash the integral of the median is not equal to 1) with the confidence interval of 1 sigma (68%), 2 sigma (95%) and 3 sigma (99.7%) as well as in green, the kernel density estimate of the input values, the detected expression (``log(1 + targetSum * raw / total UMI)``).\n+\n+            - When the 2d version is used, it displays the PDF as a heatmap as well as a projection on the x and y axis. On the projection, the confidence interval 68% is indicated as a shaded area as well as the mean with a solid red line and the median with a dashed black line. On the top right corner, the correlation is indicated with the confidence interval 68% as well as a confidence interval on the one-sided p-value (the probability that the correlation is the opposite sign of the mean, one sigma confidence interval).\n+\n+            Usually you should run baredSC_1d or baredSC_2d with 1 to 4 gaussians. Then you combine the different models with combineMultipleModels_1d or combineMultipleModels_2d.\n+            \n+        ]]></help>\n+        <citations>\n+            <citation type="doi">10.1186/s12859-021-04507-8</citation>\n+        </citations>\n+    </xml>\n+    <xml name="edam_topics">\n+        <edam_topics>\n+            <edam_topic>topic_3170</edam_topic>\n+            <edam_topic>topic_4028</edam_topic>\n+            <edam_topic>topic_2269</edam_topic>\n+        </edam_topics>\n+    </xml>\n+    <xml name="macro_input_counts">\n+        <conditional name="input_counts">\n+            <param name="filetype" type="select" label="Input type">\n+                <option value="tabular">Tabular</option>\n+                <option value="anndata">Anndata (for example from Scanpy)</option>\n+            </param>\n+            <when value="tabular">\n+                <param argument="--input" type="data" format="tabular" label="Input table (with header)" help="Expected format is one line per cell, columns with raw counts and one column \'nCount_RNA\' with the total number of UMI per cell (optionally other meta data to filter)" />\n+            </when>\n+            <when value="anndata">\n+                <param argument="--inputAnnData" type="data" format="anndata" label="AnnData containing raw counts" />\n+            </when>\n+        </conditional>\n+    </xml>\n+    <xml name="macro_single_gene">\n+        <param argument="--geneColName" type="text" value="" label="Name of the column with gene counts."/>\n+'..b'tadata1ColName\'\n+            --metadata1Values \'$filter.metadata1Values\'\n+            --metadata2ColName \'$filter.metadata2ColName\'\n+            --metadata2Values \'$filter.metadata2Values\'\n+            --metadata3ColName \'$filter.metadata3ColName\'\n+            --metadata3Values \'$filter.metadata3Values\'\n+        #end if    \n+    ]]></token>\n+    <token name="@MCMC_1D@" ><![CDATA[\n+        --xmin $MCMC.xmin\n+        --xmax $MCMC.xmax\n+        --xscale \'$MCMC.scale.type\'\n+        #if str( $MCMC.scale.type ) == "Seurat":\n+            --targetSum $MCMC.scale.targetSum\n+        #end if\n+        --nx $MCMC.nx\n+        --minScale $MCMC.minScalex\n+        --seed $MCMC.seed\n+    ]]></token>\n+    <token name="@MCMC_2D@" ><![CDATA[\n+        --xmin $MCMC.xmin\n+        --xmax $MCMC.xmax\n+        --nx $MCMC.nx\n+        --minScalex $MCMC.minScalex\n+        --ymin $MCMC.ymin\n+        --ymax $MCMC.ymax\n+        --ny $MCMC.ny\n+        --minScaley $MCMC.minScaley\n+        --scale \'$MCMC.scale.type\'\n+        #if str( $MCMC.scale.type ) == "Seurat":\n+            --targetSum $MCMC.scale.targetSum\n+        #end if\n+        --seed $MCMC.seed\n+    ]]></token>\n+    <token name="@BAREDSC_COMMON@" ><![CDATA[\n+        --nnorm $MCMC.nnorm\n+        --nsampMCMC $MCMC.nsampMCMC\n+        #if str( $MCMC.automaticRestart.set_minNeff ) == "yes":\n+            --minNeff $MCMC.automaticRestart.minNeff\n+        #end if\n+    ]]></token>\n+    <token name="@PLOTS@" ><![CDATA[\n+        #if str( $plots.title ) != \'\':\n+            --title \'$plots.title\'\n+        #end if\n+        #if $plots.removeFirstSamples != -1:\n+            --removeFirstSamples $plots.removeFirstSamples\n+        #end if\n+        --nsampInPlot $plots.nsampInPlot\n+    ]]></token>\n+    <token name="@PRETTYBINS_1D@" ><![CDATA[\n+        #if $plots.prettyBins != -1:\n+            --prettyBins $plots.prettyBins\n+        #end if    \n+    ]]></token>\n+    <token name="@PRETTYBINS_SPLITY_COLORSCALE_2D@" ><![CDATA[\n+        #if $plots.prettyBinsx != -1:\n+            --prettyBinsx $plots.prettyBinsx\n+        #end if\n+        #if $plots.prettyBinsy != -1:\n+            --prettyBinsy $plots.prettyBinsy\n+        #end if\n+        ## splity is space separated floats\n+        #if str($plots.splity) != \'\':\n+            --splity $plots.splity\n+        #end if\n+        #if str($plots.log1pColorScale) != \'\':\n+            \'$plots.log1pColorScale\'\n+        #end if\n+    ]]></token>\n+    <token name="@ADVANCED_COMMON_X@" ><![CDATA[\n+        --osampx $advanced.osampx\n+        --osampxpdf $advanced.osampxpdf\n+        --coviscale $advanced.coviscale\n+        --nis $advanced.nis\n+    ]]></token>\n+    <token name="@ADVANCED_COMMON_COMPLEMENT_2D@" ><![CDATA[\n+        --osampy $advanced.osampy\n+        --osampypdf $advanced.osampypdf\n+        --scalePrior $advanced.scalePrior\n+    ]]></token>\n+    <token name="@ADVANCED_BAREDSC_COMMON@" ><![CDATA[\n+        #if str( $advanced.burn.custom ) == "yes":\n+            #if str( $advanced.burn.nsampBurnMCMC ) != "-1":\n+                --nsampBurnMCMC $advanced.burn.nsampBurnMCMC\n+            #end if\n+            --T0BurnMCMC $advanced.burn.T0BurnMCMC\n+        #end if\n+    ]]></token>\n+    <token name="@COMBINE_OUTPUTS@"><![CDATA[\n+        --outputs\n+        #for $i, $output in enumerate($MCMC.outputs):\n+            $i\n+        #end for\n+    ]]></token>\n+    <token name="@ORDER_OUTPUTS_1D@"><![CDATA[\n+        mv baredSC_pdf.txt output &&\n+        mv baredSC.$plots.image_file_format baredSC &&\n+        gunzip baredSC_means.txt.gz\n+    ]]></token>\n+    <token name="@ORDER_OUTPUTS_2D@"><![CDATA[\n+        mv baredSC_pdf2d.txt output &&\n+        mv baredSC_pdf2d_flat.txt output &&\n+        mv baredSC.$plots.image_file_format baredSC\n+        #if str($plots.splity) != \'\':\n+            #for $value in str($plots.splity).split(\' \'):\n+                && mv baredSC_split\'$value\'.txt baredSC_split\'$value\'_pdf.txt\n+            #end for\n+        #end if\n+    ]]></token>\n+</macros>\n\\ No newline at end of file\n'
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_1gauss.npz
b
Binary file test-data/2d_small_1gauss.npz has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_1gauss.png
b
Binary file test-data/2d_small_1gauss.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss.npz
b
Binary file test-data/2d_small_2gauss.npz has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss.png
b
Binary file test-data/2d_small_2gauss.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_convergence.png
b
Binary file test-data/2d_small_2gauss_convergence.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_corr.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_small_2gauss_corr.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,2 @@
+mean median low high pval error
+0.902971588834658 0.9030205591240211 0.8808479278110088 0.9249439927709214 0.007859609399117668 0.007859609399117668
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_individuals.png
b
Binary file test-data/2d_small_2gauss_individuals.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_median.png
b
Binary file test-data/2d_small_2gauss_median.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_neff.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_small_2gauss_neff.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,1 @@
+127.23278590819812
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_p.png
b
Binary file test-data/2d_small_2gauss_p.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_p.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_small_2gauss_p.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,12 @@
+name low median high
+xy_mux0 0.3421735749355604 0.36400142067277885 0.38414407467979367
+xy_muy0 0.3602753667616778 0.3797834812619477 0.3997083085113705
+xy_scalex0 0.128373484740148 0.13754306887169315 0.154475100117708
+xy_scaley0 0.11267797849970684 0.13404387324572375 0.1636497454683696
+xy_corr0 -0.056936609709389754 0.20134274953150133 0.4823245136404016
+xy_amp1 0.4683630960583824 0.4935067017747425 0.5184297378623833
+xy_mux1 0.9950898144773097 1.0128387895441175 1.0323483080264237
+xy_muy1 0.9882799141570242 1.0073935336038764 1.0262810747590292
+xy_scalex1 0.12680233831791496 0.1317913482460104 0.1439711554909536
+xy_scaley1 0.10695025120288884 0.11525119819932761 0.13138587521715606
+xy_corr1 0.27825576244514594 0.5311286235242539 0.7541979912185461
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_pdf2d.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_small_2gauss_pdf2d.txt Mon Oct 02 13:24:52 2023 +0000
b
b'@@ -0,0 +1,51 @@\n+xy\t0.025\t0.07500000000000001\t0.125\t0.17500000000000002\t0.225\t0.275\t0.32500000000000007\t0.375\t0.42500000000000004\t0.475\t0.525\t0.5750000000000001\t0.625\t0.675\t0.7250000000000001\t0.775\t0.8250000000000001\t0.875\t0.925\t0.9750000000000001\t1.025\t1.0750000000000002\t1.125\t1.1750000000000003\t1.225\t1.275\t1.3250000000000002\t1.375\t1.4250000000000003\t1.475\t1.525\t1.5750000000000002\t1.625\t1.6750000000000003\t1.725\t1.775\t1.8250000000000002\t1.875\t1.9250000000000003\t1.975\t2.0250000000000004\t2.075\t2.125\t2.175\t2.225\t2.2750000000000004\t2.325\t2.375\t2.4250000000000003\t2.475\n+2.500000000000001527e-02\t3.891025267313832220e-02\t6.619362271397291608e-02\t9.962032325456805781e-02\t1.329150165261198668e-01\t1.585042114885953923e-01\t1.698631835606070062e-01\t1.643453093470616522e-01\t1.440979701449102424e-01\t1.147229982887641442e-01\t8.297593836730705885e-02\t5.453632221673256758e-02\t3.260197304273629304e-02\t1.775351906881531341e-02\t8.824550099856330251e-03\t4.019653022896153491e-03\t1.693722380320050562e-03\t6.723899559109907872e-04\t2.583744938047211970e-04\t9.884221137200835654e-05\t3.829570838672849967e-05\t1.502388652651616482e-05\t5.894092358324372588e-06\t2.277322359194214590e-06\t8.559947513669245676e-07\t3.105071541723489091e-07\t1.082088117846700032e-07\t3.615050191473302866e-08\t1.156966286349393425e-08\t3.547337713663130267e-09\t1.042332291023246327e-09\t2.936144188469572462e-10\t7.930373196462731033e-11\t2.053722404340254589e-11\t5.098316782057429645e-12\t1.212828418944710247e-12\t2.763685365599003245e-13\t6.030215549885882683e-14\t1.259510436836432442e-14\t2.517728221803534345e-15\t4.816310522435080954e-16\t8.817026885803886373e-17\t1.544792003475728112e-17\t2.590716812703897422e-18\t4.159550638707556230e-19\t6.394785204636262863e-20\t9.415173774848427457e-21\t1.327729810923157728e-21\t1.793526816745113172e-22\t2.320809782904078269e-23\t2.876730845976302486e-24\n+7.500000000000001110e-02\t6.643484693575134326e-02\t1.183006981844331729e-01\t1.870334996452306109e-01\t2.602824894000552725e-01\t3.208092962187330532e-01\t3.530296873517885792e-01\t3.485449341260358214e-01\t3.098527431345946548e-01\t2.486508709255841332e-01\t1.803438080812760358e-01\t1.182696247307003096e-01\t7.013448855438601803e-02\t3.761041573495765844e-02\t1.826430708113421092e-02\t8.075934593136612633e-03\t3.293371397591066205e-03\t1.265958974562858769e-03\t4.719598896863413484e-04\t1.754024667681958573e-04\t6.606852765071759824e-05\t2.524704536692301126e-05\t9.686290355684962265e-06\t3.679077132580878796e-06\t1.366629307272457547e-06\t4.922366751994848102e-07\t1.710211775214381808e-07\t5.715672258023800778e-08\t1.835103562301085251e-08\t5.657250573536754573e-09\t1.674237655582708414e-09\t4.755900425933429292e-10\t1.296452060781358147e-10\t3.390370930144212431e-11\t8.502204712729585186e-12\t2.043744190884452216e-12\t4.707199491813082972e-13\t1.038485908907196685e-13\t2.194020592867262314e-14\t4.438371167716276240e-15\t8.596588337515135905e-16\t1.594236871632899937e-16\t2.830912777714656474e-17\t4.813690310992910765e-18\t7.838603461554066503e-19\t1.222452146969727692e-19\t1.825867479967406519e-20\t2.611843420933754214e-21\t3.578038081665623504e-22\t4.693805421077943979e-23\t5.895675498776352740e-24\n+1.250000000000000278e-01\t1.051092459126146961e-01\t1.966928020901290708e-01\t3.268527238483234787e-01\t4.774922711687904497e-01\t6.121451034480128728e-01\t6.940302635696841138e-01\t7.008011953746924316e-01\t6.327642743347238197e-01\t5.122592048982417268e-01\t3.724164539772137905e-01\t2.432313027328303900e-01\t1.426556598345526949e-01\t7.514893726870940094e-02\t3.566505733111000054e-02\t1.537704344262894388e-02\t6.116063486705410925e-03\t2.294925944790045854e-03\t8.345625294036262531e-04\t3.017899410916226061e-04\t1.103951145719827496e-04\t4.101821763562869523e-05\t1.536937012505641887e-05\t5.735203619215257598e-06\t2.105293649026688872e-06\t7.531120904003529380e-07\t2.609126896498768556e-07\t8.721928811276988446e-08\t2.807461767329103996e-08\t8.691693342174174952e-09\t2.586345293985255456e-09\t7.393324391977500756e-10\t2.029380174373139609e-10\t5.346360008139223221e-11\t1.351229484442190466e-11'..b'7040694340274685e-18\t1.352683999664027101e-19\t2.233880727327099937e-20\t6.015364889538180105e-21\t1.522194703739943963e-21\t2.917476454271716322e-22\t5.243974032320023343e-23\n+2.375000000000000000e+00\t4.388195368366521442e-19\t1.671506025382529543e-18\t6.229765325347545717e-18\t2.105336918967188615e-17\t6.278555410017589079e-17\t1.639250631837221005e-16\t3.743789285820746140e-16\t7.496484696967835052e-16\t1.321773825353704208e-15\t2.064500076074844191e-15\t2.877726069056643685e-15\t3.609349866953973942e-15\t4.105728290305388638e-15\t4.262007623618967680e-15\t4.050444851776926856e-15\t3.524028773463642917e-15\t2.799948700412575740e-15\t2.025145871788806529e-15\t1.331523783793975740e-15\t7.994498677026531681e-16\t4.471346086813425646e-16\t2.460015038128927040e-16\t1.459647751812339766e-16\t9.728682956516551228e-17\t6.685661857990702758e-17\t4.432618090894337802e-17\t3.246475200005258129e-17\t3.269817719724199626e-17\t4.075921565023788040e-17\t5.046099433093516435e-17\t5.727527152615365706e-17\t5.996672573574389461e-17\t6.095074740235865311e-17\t6.461790423958404874e-17\t7.321735849601882793e-17\t8.262643983910216582e-17\t8.422325253461393166e-17\t7.512672597931765004e-17\t6.617561101351289782e-17\t6.538445205624774633e-17\t5.111031621943633462e-17\t2.178188384974001937e-17\t4.726424170356313971e-18\t5.943522092520003635e-19\t5.601173801802059831e-20\t5.560237658109037348e-21\t9.492385285542926078e-22\t2.692603329480852931e-22\t7.031463560866556195e-23\t1.455205872253979626e-23\n+2.425000000000000266e+00\t8.571564022902764806e-20\t3.395439739908100516e-19\t1.300848088687509419e-18\t4.492294622840806904e-18\t1.365543666828820425e-17\t3.630156830002293030e-17\t8.437556840174051521e-17\t1.719088623969656934e-16\t3.084028994035439583e-16\t4.901756897508358145e-16\t6.954939245671626415e-16\t8.883658474241139825e-16\t1.029794636777511356e-15\t1.090128984112849601e-15\t1.057168644485305165e-15\t9.389506807514792310e-16\t7.616416845394756049e-16\t5.621394798219496508e-16\t3.765337410882662026e-16\t2.292452277214820660e-16\t1.283642507996977865e-16\t6.853291123593672494e-17\t3.768342026303448032e-17\t2.307803313501243235e-17\t1.519192615928129011e-17\t9.534597858565881257e-18\t5.591172693423777730e-18\t3.803442160832206803e-18\t3.757264283381587666e-18\t4.478361656915204916e-18\t5.174891340753303504e-18\t5.500859347528055462e-18\t5.539309038841850533e-18\t5.658811000252308673e-18\t6.217175604513879417e-18\t7.160835889399820430e-18\t7.833684418986813088e-18\t7.501494813931311113e-18\t6.357012408533341629e-18\t5.738709324323115522e-18\t5.666212930823022056e-18\t3.857326183522140865e-18\t1.350534483427198886e-18\t2.424979002822490245e-19\t2.629368026319186048e-20\t2.244619057368868742e-21\t2.223672272825115933e-22\t4.370650857847124275e-23\t1.344292542151976491e-23\t3.760604778948504676e-24\n+2.475000000000000089e+00\t1.621105492714958480e-20\t6.670586069916486759e-20\t2.624185565246232221e-19\t9.254962246112966850e-19\t2.866737297683849356e-18\t7.758560084042475930e-18\t1.835124971200113702e-17\t3.804196814888727620e-17\t6.943670617385706399e-17\t1.123011521110732181e-16\t1.621880133113295537e-16\t2.109710312926076754e-16\t2.492120913452456150e-16\t2.690276800739296073e-16\t2.662263565714788191e-16\t2.414027798437603098e-16\t1.999498241791764298e-16\t1.506523468799500421e-16\t1.029047893501557564e-16\t6.369930956196674557e-17\t3.597208775848153126e-17\t1.896702384562388817e-17\t9.887308089617070670e-18\t5.562792909840032236e-18\t3.467606600929517228e-18\t2.159183945797982487e-18\t1.200782582161042977e-18\t6.195436213305423752e-19\t3.976685628676712431e-19\t3.857767153251908125e-19\t4.354851688234009867e-19\t4.694646198431953851e-19\t4.731792156614820582e-19\t4.700294740951533729e-19\t4.946985957992040183e-19\t5.631367928923222179e-19\t6.453060686425745162e-19\t6.717708568966487404e-19\t6.009024143169283498e-19\t4.923358493148858372e-19\t4.662008114386063913e-19\t4.428227247828746297e-19\t2.550204992080155717e-19\t7.297036822235744721e-20\t1.093605303531067023e-20\t1.041595429534497928e-21\t8.470178050411832510e-23\t9.965011157153024311e-24\t2.539622475898041699e-24\t8.429237166204116204e-25\n'
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_small_2gauss_pdf2d_flat.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_small_2gauss_pdf2d_flat.txt Mon Oct 02 13:24:52 2023 +0000
b
b'@@ -0,0 +1,2501 @@\n+x\ty\tlow\tmean\thigh\tmedian\n+0.025\t0.025000000000000015\t0.0028275171230156943\t0.03891025267313832\t0.07345381655886798\t0.018883354134479356\n+0.07500000000000001\t0.025000000000000015\t0.006166221033717528\t0.06619362271397292\t0.12567989759768028\t0.0337612674083914\n+0.125\t0.025000000000000015\t0.011978878035529729\t0.09962032325456806\t0.1863350204202721\t0.05496303551549483\n+0.17500000000000002\t0.025000000000000015\t0.01877780272265707\t0.13291501652611987\t0.24891690268353578\t0.07993148214679692\n+0.225\t0.025000000000000015\t0.02410756070471183\t0.1585042114885954\t0.29339138509437995\t0.0995079399448785\n+0.275\t0.025000000000000015\t0.02532984131849599\t0.169863183560607\t0.30805715571507436\t0.10746902259407222\n+0.32500000000000007\t0.025000000000000015\t0.023613393634799477\t0.16434530934706165\t0.3006682022616266\t0.10469285364724716\n+0.375\t0.025000000000000015\t0.019670559863889363\t0.14409797014491024\t0.27049426273317895\t0.08673567821541447\n+0.42500000000000004\t0.025000000000000015\t0.012788182883589922\t0.11472299828876414\t0.2228113815127328\t0.06301151553931411\n+0.475\t0.025000000000000015\t0.006537776975833527\t0.08297593836730706\t0.16142273408061317\t0.04113263922283775\n+0.525\t0.025000000000000015\t0.002895780791935031\t0.05453632221673257\t0.10766389836532349\t0.02398314733204404\n+0.5750000000000001\t0.025000000000000015\t0.0010370259649472387\t0.03260197304273629\t0.06480922628228036\t0.012019613818710188\n+0.625\t0.025000000000000015\t0.00027898259498954217\t0.017753519068815313\t0.033892216663849145\t0.005343071189537988\n+0.675\t0.025000000000000015\t6.724666645975823e-05\t0.00882455009985633\t0.01569156940323104\t0.00202601643085585\n+0.7250000000000001\t0.025000000000000015\t1.2041623854199181e-05\t0.0040196530228961535\t0.00629222212450801\t0.0007094145496271375\n+0.775\t0.025000000000000015\t2.022922615331913e-06\t0.0016937223803200506\t0.00237857816866685\t0.0002165095535088674\n+0.8250000000000001\t0.025000000000000015\t2.9839051283922713e-07\t0.0006723899559109908\t0.000785755321359725\t5.657618073183036e-05\n+0.875\t0.025000000000000015\t3.3008043233790196e-08\t0.0002583744938047212\t0.00022839911336038367\t1.2431516799992795e-05\n+0.925\t0.025000000000000015\t3.3053643500745195e-09\t9.884221137200836e-05\t5.881281859435004e-05\t2.4616229055495425e-06\n+0.9750000000000001\t0.025000000000000015\t2.9428889698194374e-10\t3.82957083867285e-05\t1.3765115238751885e-05\t4.055343653562609e-07\n+1.025\t0.025000000000000015\t2.064387009795184e-11\t1.5023886526516165e-05\t2.6853708648058713e-06\t6.092093434681246e-08\n+1.0750000000000002\t0.025000000000000015\t1.4722370597000242e-12\t5.894092358324373e-06\t4.302731050889479e-07\t7.52321040788387e-09\n+1.125\t0.025000000000000015\t1.2092724480936643e-13\t2.2773223591942146e-06\t6.548462248085567e-08\t8.773904252515071e-10\n+1.1750000000000003\t0.025000000000000015\t6.2348799285296446e-15\t8.559947513669246e-07\t9.351271603479492e-09\t8.53584938162264e-11\n+1.225\t0.025000000000000015\t3.5677381106047285e-16\t3.105071541723489e-07\t1.2107914093422262e-09\t7.436403053238842e-12\n+1.275\t0.025000000000000015\t2.5578616694107807e-17\t1.0820881178467e-07\t1.2784274270372268e-10\t5.071125386915198e-13\n+1.3250000000000002\t0.025000000000000015\t1.0260600281822056e-18\t3.615050191473303e-08\t1.3449313594594362e-11\t2.972123526454833e-14\n+1.375\t0.025000000000000015\t3.7039418594411067e-20\t1.1569662863493934e-08\t1.2183844129360137e-12\t1.5596194741125646e-15\n+1.4250000000000003\t0.025000000000000015\t1.2102313342013408e-21\t3.5473377136631303e-09\t1.0384367959432388e-13\t8.357754248596118e-17\n+1.475\t0.025000000000000015\t3.9303115866406374e-23\t1.0423322910232463e-09\t8.229155280840018e-15\t3.390501069155369e-18\n+1.525\t0.025000000000000015\t1.0047933972630657e-24\t2.9361441884695725e-10\t5.870515400119876e-16\t1.480287378743775e-19\n+1.5750000000000002\t0.025000000000000015\t2.4605178697112705e-26\t7.930373196462731e-11\t3.609255229567178e-17\t4.95598012818425e-21\n+1.625\t0.025000000000000015\t4.248958675849221e-28\t2.0537224043402546e-11\t2.0969002750897167e-18\t1.3910292422317308e-22\n+1.6750000000000003\t0.02500'..b'345098053098809e-48\n+0.675\t2.475\t3.8091766597197725e-65\t2.690276800739296e-16\t6.526448310190745e-35\t1.480195602619133e-47\n+0.7250000000000001\t2.475\t2.3771592969854215e-64\t2.662263565714788e-16\t1.773790062443437e-34\t9.047173703726691e-47\n+0.775\t2.475\t4.988422592144346e-63\t2.414027798437603e-16\t4.816244182076834e-34\t5.3377480646970545e-46\n+0.8250000000000001\t2.475\t2.3363916964867087e-62\t1.9994982417917643e-16\t1.4942695401412177e-33\t3.288609392671216e-45\n+0.875\t2.475\t1.3018520678006217e-61\t1.5065234687995004e-16\t3.619129027798108e-33\t2.324104792108237e-44\n+0.925\t2.475\t4.04588138155899e-61\t1.0290478935015576e-16\t9.172907237511743e-33\t8.666846902619233e-44\n+0.9750000000000001\t2.475\t9.845305967414265e-60\t6.369930956196675e-17\t1.496461074859211e-32\t2.7832865384984236e-43\n+1.025\t2.475\t1.293033349303233e-58\t3.597208775848153e-17\t2.8777759271058696e-32\t1.1275695818540316e-42\n+1.0750000000000002\t2.475\t2.6275576047680974e-57\t1.8967023845623888e-17\t8.273870886884713e-32\t4.3792395015718776e-42\n+1.125\t2.475\t2.853877049497956e-56\t9.88730808961707e-18\t1.6261395817554126e-31\t1.2100504935734452e-41\n+1.1750000000000003\t2.475\t4.5048558129388604e-55\t5.562792909840032e-18\t3.493594330887723e-31\t3.5490820078796784e-41\n+1.225\t2.475\t1.4671619706830216e-53\t3.4676066009295172e-18\t6.379190207519922e-31\t1.1237161343014799e-40\n+1.275\t2.475\t2.5216984116045096e-52\t2.1591839457979825e-18\t1.2928108716950301e-30\t3.2263615338868193e-40\n+1.3250000000000002\t2.475\t1.640986499373436e-51\t1.200782582161043e-18\t2.331796994134262e-30\t1.0864048003008563e-39\n+1.375\t2.475\t1.8947741770675045e-50\t6.195436213305424e-19\t3.765846571495058e-30\t2.0956631034465537e-39\n+1.4250000000000003\t2.475\t2.393775220629594e-49\t3.9766856286767124e-19\t7.002345582757834e-30\t6.432075808649665e-39\n+1.475\t2.475\t2.1740418570739486e-48\t3.857767153251908e-19\t9.015969019000993e-30\t1.7287832009103091e-38\n+1.525\t2.475\t2.0565294553638698e-47\t4.35485168823401e-19\t1.3682142130178459e-29\t2.7790891765849295e-38\n+1.5750000000000002\t2.475\t9.619418651609598e-47\t4.694646198431954e-19\t1.7386869895648372e-29\t5.860586454617413e-38\n+1.625\t2.475\t4.812515095437497e-46\t4.731792156614821e-19\t3.1401494643422897e-29\t1.1777792992318676e-37\n+1.6750000000000003\t2.475\t3.174186309486899e-45\t4.700294740951534e-19\t3.9525595202090304e-29\t1.9823961412719416e-37\n+1.725\t2.475\t7.014538039370068e-45\t4.94698595799204e-19\t6.660650906823616e-29\t1.922611826870485e-37\n+1.775\t2.475\t2.362262232546069e-44\t5.631367928923222e-19\t9.677166097424813e-29\t2.424815923292836e-37\n+1.8250000000000002\t2.475\t1.1318410108685484e-43\t6.453060686425745e-19\t1.1489844829911155e-28\t2.69997301176567e-37\n+1.875\t2.475\t2.1207021249606866e-43\t6.717708568966487e-19\t1.1349152256951729e-28\t2.2483673044948416e-37\n+1.9250000000000003\t2.475\t2.706603668641923e-43\t6.0090241431692835e-19\t9.686537686918842e-29\t1.9703510180181214e-37\n+1.975\t2.475\t2.227125479893856e-43\t4.923358493148858e-19\t8.072274241507064e-29\t1.6158083075112546e-37\n+2.0250000000000004\t2.475\t1.9205193974142984e-43\t4.662008114386064e-19\t5.695883940964188e-29\t1.5141126929192331e-37\n+2.075\t2.475\t1.5550933132379314e-43\t4.428227247828746e-19\t3.7799447826164555e-29\t1.0826246353611623e-37\n+2.125\t2.475\t4.737485095003875e-44\t2.5502049920801557e-19\t2.412227325070249e-29\t8.167967230511536e-38\n+2.175\t2.475\t1.7052343983381615e-44\t7.297036822235745e-20\t1.328950798724004e-29\t5.257740759745071e-38\n+2.225\t2.475\t6.358928264061567e-45\t1.093605303531067e-20\t6.852710171786249e-30\t2.594231701517555e-38\n+2.2750000000000004\t2.475\t1.593783598311689e-45\t1.041595429534498e-21\t2.4982608900807184e-30\t1.1361328515105745e-38\n+2.325\t2.475\t2.1115715928687407e-46\t8.470178050411833e-23\t8.030121053113561e-31\t3.320931549850192e-39\n+2.375\t2.475\t2.4483066342670293e-47\t9.965011157153024e-24\t1.7034399210663177e-31\t1.0991883306410254e-39\n+2.4250000000000003\t2.475\t2.2205465453156153e-48\t2.5396224758980417e-24\t4.3762986802076344e-32\t2.3092663689539515e-40\n+2.475\t2.475\t1.646645981074985e-49\t8.429237166204116e-25\t8.324766201430514e-33\t5.202229622039887e-41\n'
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2.pdf
b
Binary file test-data/2d_test2.pdf has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2_split-10.5.pdf
b
Binary file test-data/2d_test2_split-10.5.pdf has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2_split-10.5.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_test2_split-10.5.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,26 @@
+x low1 mean1 high1 low2 mean2 high2
+-1.190000000000000036e+01 2.598287065893160544e-44 1.585274254574060360e-05 1.835135945978420623e-08 1.442773755860150222e-38 2.052048995737900338e-04 1.050694500596239781e-05
+-1.169999999999999929e+01 1.349513053924564785e-35 3.977560456761572108e-05 2.486609810217813625e-07 2.466748504170547658e-30 5.624401881034365698e-04 1.120402940282793763e-04
+-1.150000000000000000e+01 4.531319753132766899e-28 9.594545582756318521e-05 2.787707049014204804e-06 5.097778958715304920e-23 1.533003143316528580e-03 8.720320356902651034e-04
+-1.130000000000000071e+01 1.690733658986313566e-21 2.235265749322897790e-04 2.374482969585529565e-05 1.053795742793207507e-16 4.139358802028571260e-03 5.220510767781757697e-03
+-1.109999999999999964e+01 3.361301828246254505e-16 5.089872288221513556e-04 1.614789875614967732e-04 1.633361785534608728e-11 1.106356465725347356e-02 2.302314040549096816e-02
+-1.090000000000000036e+01 2.058448339772207212e-12 1.156445332511210925e-03 8.883596776243168165e-04 2.781665822380890478e-07 2.944798800469059646e-02 7.659334522399686984e-02
+-1.069999999999999929e+01 7.180708905730886339e-10 2.724782550339245341e-03 3.839707915997163899e-03 3.570567554242068224e-04 7.972997990899827192e-02 1.924991874783543999e-01
+-1.050000000000000000e+01 2.118059392717893225e-08 6.994343260081681973e-03 1.408657995225120162e-02 4.255802869340978567e-02 2.291654660250422582e-01 3.743848094595008669e-01
+-1.030000000000000071e+01 8.941234004590808332e-08 1.617978372272826232e-02 3.455995943978115020e-02 4.052222164743428667e-01 6.149713032591206385e-01 8.406626287071065073e-01
+-1.009999999999999964e+01 1.156514375977928323e-07 1.946062853407351870e-02 4.251300851815696669e-02 5.021047445344608562e-01 8.530517003625158745e-01 1.247500090307065657e+00
+-9.900000000000000355e+00 4.021960876298153002e-08 9.981651013719877036e-03 2.096963599255254398e-02 2.553671179752886045e-01 4.539091434697714589e-01 6.492176009838344086e-01
+-9.699999999999999289e+00 4.119443957387565366e-09 3.215048894368956797e-03 5.534192606720483160e-03 1.038235034416979015e-02 1.511081221419017240e-01 2.887429717738879131e-01
+-9.500000000000000000e+00 4.783275649973367636e-11 1.044277882174386631e-03 1.073064916336143775e-03 9.317854808959093269e-05 5.229663615911780972e-02 1.306600310909605545e-01
+-9.300000000000000711e+00 3.003636268035327315e-14 3.573556103687697510e-04 1.824673131679265003e-04 1.884166311931892410e-05 2.171438747310733969e-02 5.122512949820227923e-02
+-9.099999999999999645e+00 1.687298123156504498e-18 1.240063631822725231e-04 2.417333162285666118e-05 1.055362682369262450e-03 3.598596370719481241e-02 7.373352487026842150e-02
+-8.900000000000000355e+00 1.179297473663470110e-23 4.278942610602769317e-05 2.544003818644774816e-06 9.707895017342857424e-02 2.705574890605504934e-01 4.568789886800128097e-01
+-8.699999999999999289e+00 6.043793715267850584e-29 1.445209254315350252e-05 1.949999051004347995e-07 8.915363695650204701e-01 1.115099752827883428e+00 1.341746322728876839e+00
+-8.500000000000000000e+00 8.345770499217446651e-35 4.741729528548845399e-06 1.211121106256519643e-08 6.866577676678010667e-01 8.766412038384161498e-01 1.071256494170844009e+00
+-8.300000000000000711e+00 7.576091502559685812e-41 1.511795234588014025e-06 6.176640004137890864e-10 4.370807590349612415e-02 1.274936086556556492e-01 2.185021944787013948e-01
+-8.099999999999999645e+00 2.595627804484775466e-47 4.699609966205265639e-07 2.267703208666181153e-11 1.300302155302189978e-04 8.544747678265863719e-03 1.697085463151461887e-02
+-7.900000000000000355e+00 2.196438072005507691e-54 1.429573987255398505e-07 6.318022541491044091e-13 2.167103791447842661e-08 5.524355961968038320e-04 4.479041476095592866e-04
+-7.699999999999999289e+00 2.989352204500754866e-62 4.269383561617732233e-08 1.190794994533467480e-14 2.202011642212695045e-13 3.632927073407426230e-05 6.302616359194336635e-06
+-7.500000000000000000e+00 1.914571918074225196e-70 1.255873286072209397e-08 2.249668958609558916e-16 2.023934352661581583e-19 2.994250821347781221e-06 6.300052419582949432e-08
+-7.299999999999998934e+00 6.093188640531772828e-79 3.649617069822763767e-09 2.805171728965951326e-18 1.501234095664622088e-26 4.724760092182685135e-07 4.382145227679361189e-10
+-7.099999999999999645e+00 1.043764176529380021e-88 1.049813103734197142e-09 2.349399036931189675e-20 1.586457343007895702e-34 1.254599427302130363e-07 2.783718903754569113e-12
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2_split-10.5_renorm.pdf
b
Binary file test-data/2d_test2_split-10.5_renorm.pdf has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2_split-9.5.pdf
b
Binary file test-data/2d_test2_split-9.5.pdf has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2_split-9.5.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/2d_test2_split-9.5.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,26 @@
+x low1 mean1 high1 low2 mean2 high2
+-1.190000000000000036e+01 4.803297546172825515e-39 2.017872702441017119e-04 1.066605775712725591e-05 2.003100457871771228e-42 1.927037187543402806e-05 1.688426799957086159e-07
+-1.169999999999999929e+01 1.161303246704032504e-30 5.496006148108198994e-04 1.099098544141658773e-04 8.889522294287208610e-34 5.261517786022843417e-05 2.222227168005508637e-06
+-1.150000000000000000e+01 3.393936980747613856e-23 1.486302373782895370e-03 8.528887702054521475e-04 3.836472170264625068e-26 1.426462253612124009e-04 2.305306783356143470e-05
+-1.130000000000000071e+01 6.806582203557863020e-17 3.980164398130769249e-03 4.981473167096432328e-03 1.546718412604563423e-19 3.827209788301017381e-04 1.737920350474747751e-04
+-1.109999999999999964e+01 1.258563200681326123e-11 1.055446125729465542e-02 2.186949458038200955e-02 5.790600436075149673e-14 1.018090628780808377e-03 1.017696831184716766e-03
+-1.090000000000000036e+01 2.188119976221788961e-07 2.789667220339803883e-02 7.282719532245499061e-02 1.549518814910123141e-09 2.707761133803449592e-03 4.489710080455356842e-03
+-1.069999999999999929e+01 3.206966582101956377e-04 7.512427747863843763e-02 1.806309476541008152e-01 2.610843065738889603e-06 7.330484980702736859e-03 1.619434140071621467e-02
+-1.050000000000000000e+01 3.921756689641812954e-02 2.152526083096725495e-01 3.550092074419232135e-01 2.726848290829387840e-04 2.090720097544487149e-02 4.642338911014470826e-02
+-1.030000000000000071e+01 3.699328651384568567e-01 5.753208278272271547e-01 7.977196431283402056e-01 2.248706504763652372e-03 5.583025915460779892e-02 1.121623967111159692e-01
+-1.009999999999999964e+01 4.633465497357979190e-01 7.895864364723708118e-01 1.154317929976509616e+00 3.512875341287357255e-03 8.292589242422288009e-02 1.678930660422758014e-01
+-9.900000000000000355e+00 2.283646993866594033e-01 4.113441340860961581e-01 5.877728380667451846e-01 1.428503109746356772e-03 5.254666039738387778e-02 1.095633358830468695e-01
+-9.699999999999999289e+00 8.887427174977266878e-03 1.329787854913495049e-01 2.540325651776247362e-01 1.037217019384402666e-04 2.134438554492266993e-02 4.641724143930865687e-02
+-9.500000000000000000e+00 3.019962819746257844e-05 4.429685143279186305e-02 1.115543510169677982e-01 3.244303170449063270e-06 9.044062608498257272e-03 1.668866422265784244e-02
+-9.300000000000000711e+00 3.442972583844935983e-08 1.543178997185759918e-02 3.742319453421989739e-02 4.591242740380319285e-06 6.639953111618982808e-03 1.137285386163181845e-02
+-9.099999999999999645e+00 1.471986463507429327e-10 5.396603749266395272e-03 9.406826521285129722e-03 7.471197086760812244e-04 3.071336632111123371e-02 6.127961189665303426e-02
+-8.900000000000000355e+00 2.721226402947793872e-12 1.873780764548195824e-03 1.813321119532720427e-03 9.521215798813713327e-02 2.687264977221154316e-01 4.553186638803994879e-01
+-8.699999999999999289e+00 4.757048390018248319e-14 6.435532137013693678e-04 2.872715585477180812e-04 8.906457999184262952e-01 1.114470651706724658e+00 1.341241343278842812e+00
+-8.500000000000000000e+00 4.137697534466050439e-16 2.143814047382074464e-04 3.587008207002907040e-05 6.865111729039222732e-01 8.764315641632429710e-01 1.071067339992828771e+00
+-8.300000000000000711e+00 9.978530950646850700e-19 6.947949580906451661e-05 3.045879549156452239e-06 4.365996815013298976e-02 1.274256409550834357e-01 2.184863667363852180e-01
+-8.099999999999999645e+00 4.288314974106641885e-22 2.222711549172525899e-05 1.705535415530815804e-07 1.255036386645525665e-04 8.522990523770715421e-03 1.693981843682318411e-02
+-7.900000000000000355e+00 1.632483644593243785e-26 7.056652832482650939e-06 6.950482326105564616e-09 1.820434986660159292e-08 5.455219007630228564e-04 4.285414777976365814e-04
+-7.699999999999999289e+00 5.728955648096030054e-32 2.230623699098163889e-06 2.098712040575874791e-10 1.638727215250439284e-13 3.414134087059198349e-05 4.649309106539978739e-06
+-7.500000000000000000e+00 9.162447402609750420e-39 7.069527669934216702e-07 4.679403536918906080e-12 1.240923798772455686e-19 2.299856787215091993e-06 3.322760873763373091e-08
+-7.299999999999998934e+00 8.312775421822221248e-47 2.262099434286902207e-07 7.406853831987881032e-14 7.367198754646305872e-27 2.499156828594033445e-07 1.844675618561730004e-10
+-7.099999999999999645e+00 9.434110606721632665e-56 7.331074235484486987e-08 9.422420351779523587e-16 4.918119958221378093e-35 5.319901347910126341e-08 7.256150775317469917e-13
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/2d_test2_split-9.5_renorm.pdf
b
Binary file test-data/2d_test2_split-9.5_renorm.pdf has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/generate_small_inputs_for_combine.sh
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/generate_small_inputs_for_combine.sh Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,6 @@
+for nnorm in 1 2; do
+    baredSC_1d --input test-data/nih3t3_generated_2d_2.txt --geneColName 0.5_0_0_0.5_x --nnorm ${nnorm} --output test-data/small_${nnorm}gauss --figure test-data/small_${nnorm}gauss.png --nx 10 --prettyBins 100 --nsampMCMC 20000 --force
+done
+for nnorm in 1 2; do
+    baredSC_2d  --input test-data/nih3t3_generated_2d_2.txt --geneXColName '0.5_0_0_0.5_x' --geneYColName '0.5_0_0_0.5_y' --nnorm ${nnorm} --output test-data/2d_small_${nnorm}gauss --figure test-data/2d_small_${nnorm}gauss.png --nx 10 --ny 12 --prettyBinsx 50 --prettyBinsy 50 --nsampMCMC 20000 --force
+done
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/nih3t3_generated_2d_2.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/nih3t3_generated_2d_2.txt Mon Oct 02 13:24:52 2023 +0000
b
b'@@ -0,0 +1,2362 @@\n+nCount_RNA\tgroup\t0.5_0_0_0.5_x\t0.5_0_0_0.5_y\t0.5_0_0_0.5_x_expression\t0.5_0_0_0.5_y_expression\t0.5_0_0_0.5_group\t0_0.5_0.5_0_x\t0_0.5_0.5_0_y\t0_0.5_0.5_0_x_expression\t0_0.5_0.5_0_y_expression\t0_0.5_0.5_0_group\n+20036\tgroup3\t0\t2\t0.35510019026207057\t0.6049293077078384\t0.0\t1\t0\t1.1142752911813003\t0.09168608678486606\t2.0\n+20682\tgroup2\t1\t2\t1.012767750238586\t0.923452095890812\t3.0\t0\t1\t0.19728620786115628\t1.0743203122055762\t1.0\n+36511\tgroup3\t0\t1\t0.13462720393366712\t0.3900268299651572\t0.0\t8\t1\t1.1859765754934843\t0.37579805653745385\t2.0\n+12534\tgroup3\t0\t1\t0.390793063729482\t0.6152979243684825\t0.0\t0\t1\t0.47161531393189604\t0.7961539475432636\t1.0\n+18220\tgroup3\t0\t1\t0.35488704567606544\t0.4460804418942297\t0.0\t1\t3\t0.35586491373283113\t1.0341531453670259\t1.0\n+7350\tgroup2\t0\t0\t0.35013103204542273\t0.4398811772343231\t0.0\t4\t0\t1.1122591752924653\t0.3738889060220699\t2.0\n+13052\tgroup3\t0\t1\t0.3204381207408692\t0.27726718838397035\t0.0\t1\t0\t1.007742102318972\t0.23711647277930997\t2.0\n+16425\tgroup1\t0\t2\t0.27233072744972736\t0.6228530625447832\t0.0\t1\t4\t0.4492853357933467\t1.0670730491038631\t1.0\n+13753\tgroup3\t1\t0\t0.41311449227036884\t0.07138384054975677\t0.0\t0\t1\t0.626807152311016\t0.9996257626088465\t1.0\n+14036\tgroup2\t3\t3\t1.0547981865540972\t1.1129859604818588\t3.0\t2\t0\t0.9270108741415989\t0.24178169986959502\t2.0\n+12190\tgroup3\t1\t1\t0.3096655708377394\t0.29225200862529965\t0.0\t0\t0\t0.38727676159757574\t0.9243906753918077\t1.0\n+12761\tgroup3\t1\t1\t0.9449256437955227\t0.9182142062450089\t3.0\t1\t0\t0.7723190576595443\t0.2783126172084158\t2.0\n+25426\tgroup2\t0\t0\t0.2975714102053989\t0.35471875401999436\t0.0\t4\t0\t0.9824819284205879\t0.34953723259002023\t2.0\n+23016\tgroup3\t6\t3\t1.0027144114277171\t1.1395788534214049\t3.0\t7\t0\t0.9752734452388224\t0.2903407867021942\t2.0\n+15267\tgroup3\t0\t2\t0.4137255881887701\t0.4713894245776208\t0.0\t4\t0\t1.0056867608396607\t0.2037120518899986\t2.0\n+15680\tgroup3\t2\t3\t1.1242353697500371\t0.997094707155079\t3.0\t2\t0\t0.9945809753715524\t0.32561168890034564\t2.0\n+9677\tgroup1\t0\t2\t0.38433888509789293\t0.26924630090085466\t0.0\t2\t0\t1.0139206398638207\t0.25468842709092343\t2.0\n+17116\tgroup3\t6\t2\t1.1260708689336378\t1.0283867042552086\t3.0\t0\t1\t0.3755620820033411\t0.9151439204012725\t1.0\n+16540\tgroup2\t2\t2\t0.2375471054570028\t0.5118387693140509\t0.0\t1\t1\t1.0108967882175413\t0.521078369374232\t2.0\n+13464\tgroup3\t0\t0\t0.37946109742502904\t0.6412769300919429\t0.0\t5\t1\t1.0709186340652448\t0.5818657252701785\t2.0\n+18591\tgroup2\t2\t1\t0.885248931811368\t1.0725436501197576\t3.0\t0\t4\t0.19354761001903656\t1.0708010759899333\t1.0\n+15630\tgroup3\t2\t2\t0.9183974858275997\t1.0046183461456941\t3.0\t2\t0\t1.0560335156342886\t0.15623088848708816\t2.0\n+35148\tgroup3\t0\t1\t0.2631938450790857\t0.2532211457679464\t0.0\t8\t4\t1.0304645965274597\t0.4403752487309259\t2.0\n+17782\tgroup2\t0\t2\t0.8369208902686032\t0.9466374992793527\t3.0\t0\t1\t0.9813666886006327\t0.3981094025029501\t2.0\n+9322\tgroup3\t1\t1\t0.9851083710406641\t0.911367506715367\t3.0\t1\t1\t0.9443911060791809\t0.5338814108838679\t2.0\n+19299\tgroup3\t2\t3\t0.8117383965258445\t0.8599897202718597\t3.0\t2\t2\t1.1940742508021156\t0.4290754541849754\t2.0\n+17503\tgroup3\t1\t1\t0.47391652712246773\t0.4059821338901858\t0.0\t0\t3\t0.2981814548415668\t0.904231509972303\t1.0\n+10534\tgroup3\t1\t0\t0.5068121948249535\t0.08570442449200488\t0.0\t2\t0\t0.9570484465158777\t0.3902176965410252\t2.0\n+22958\tgroup2\t0\t3\t0.373408112057528\t0.3392370374931661\t0.0\t2\t4\t0.37517471469628294\t0.8648950661232981\t1.0\n+16941\tgroup1\t3\t3\t0.989158493167134\t1.0638544328754207\t3.0\t0\t2\t0.43334458123234876\t1.0482739992395556\t1.0\n+19740\tgroup2\t2\t1\t0.5028340469974669\t0.16366680056120142\t0.0\t0\t1\t0.9223751241036221\t0.4415693556726121\t2.0\n+11104\tgroup1\t1\t3\t0.20010946733914475\t0.5892171263132788\t0.0\t2\t0\t0.8492243173547003\t0.33920635766237656\t2.0\n+14576\tgroup3\t2\t2\t1.0383010599812037\t0.9852349306377165\t3.0\t1\t4\t0.4820012917941298\t1.1678685734556675\t1.0\n+14945\tgroup3\t2\t5\t0.9412658784110647\t0.8998571591368998\t3.0\t4\t0\t1.2117011086514062\t0.4387060710338841\t2.0\n+18929\tgroup3\t7\t6\t1.0457251624613944\t1.117467206998317\t3.0\t4\t2\t0.8863414072909881\t0.49187549540392894\t2.0\n+8625\tgroup3\t0\t0\t0.376853079493063\t0.35514'..b'873758031\t0.0\t5\t0\t0.9880431112507074\t0.44681668942699504\t2.0\n+13577\tgroup3\t0\t0\t0.22028805303091836\t0.3869016265751662\t0.0\t3\t1\t1.0808794218188824\t0.4058835355052489\t2.0\n+24151\tgroup2\t0\t0\t0.3561175039850726\t0.40388367791259805\t0.0\t1\t3\t0.5015529369378173\t0.8123570822161652\t1.0\n+19289\tgroup3\t4\t6\t0.9989512208932608\t1.133481240919195\t3.0\t4\t0\t0.9379042447925486\t0.2782680907382488\t2.0\n+15398\tgroup3\t2\t0\t0.3248179986863013\t0.28158780662914434\t0.0\t3\t1\t0.8794089423732879\t0.5498341568627245\t2.0\n+19236\tgroup3\t4\t4\t0.915300011055525\t0.9699695656272898\t3.0\t0\t5\t0.23380390154246436\t1.1243183385615119\t1.0\n+15455\tgroup3\t1\t1\t0.3017293584017071\t0.48486454938668394\t0.0\t0\t4\t0.290456258804643\t0.7791056774358668\t1.0\n+11675\tgroup3\t2\t2\t0.9956920929753222\t1.0062914052461407\t3.0\t2\t1\t1.0236248904354341\t0.5147278904595354\t2.0\n+22414\tgroup2\t6\t6\t1.166995117137059\t0.9362789495828017\t3.0\t5\t0\t0.9246383340396649\t0.42353998200203613\t2.0\n+14686\tgroup3\t1\t1\t0.33116440192692126\t0.5321620425714314\t0.0\t3\t2\t1.0488950175204788\t0.7005445529986285\t2.0\n+12377\tgroup1\t0\t1\t0.3707241437046496\t0.5782552867911115\t0.0\t0\t0\t0.3498105645725788\t0.7749565315524524\t1.0\n+19841\tgroup2\t3\t3\t0.9751006803371836\t1.0003829042622288\t3.0\t3\t0\t0.9966346104575061\t0.40650031140045695\t2.0\n+13295\tgroup3\t0\t1\t0.7769046103487025\t1.0046830864496405\t3.0\t1\t0\t0.19307271521541822\t0.7358620346863765\t1.0\n+15144\tgroup3\t0\t0\t0.49759907343045384\t0.3330015439110453\t0.0\t1\t0\t1.0500221733049644\t0.14476202751550837\t2.0\n+17480\tgroup3\t2\t2\t1.0297667758582996\t1.1836536481499413\t3.0\t0\t1\t0.46267305970500683\t0.8276249806065282\t1.0\n+16251\tgroup3\t0\t2\t0.4124308677803242\t0.6828302048780339\t0.0\t0\t4\t0.46673405966078196\t0.9710969217818207\t1.0\n+13561\tgroup3\t0\t2\t1.0787246582949683\t0.9841428767980445\t3.0\t1\t0\t0.8570738517467593\t0.3413881658359805\t2.0\n+11390\tgroup3\t2\t2\t0.9604294894508263\t1.2253630529191466\t3.0\t0\t0\t1.012756884988206\t0.44265515147011836\t2.0\n+55708\tgroup1\t8\t12\t0.8733377369334027\t0.9801556173055537\t3.0\t1\t9\t0.3070577902385216\t0.9513703345871919\t1.0\n+16063\tgroup3\t2\t1\t0.9644819197446499\t0.8783361692912572\t3.0\t2\t1\t1.085734241123231\t0.4196712229928285\t2.0\n+16422\tgroup3\t1\t0\t0.5406844787097789\t0.3384544586482739\t0.0\t2\t6\t0.5391804026148392\t0.9567956910514819\t1.0\n+17668\tgroup3\t4\t6\t1.0799805445421096\t1.031096050631378\t3.0\t0\t2\t0.39168343939938344\t0.9868614621641726\t1.0\n+16415\tgroup3\t1\t4\t0.9923551015783658\t1.0198173980788054\t3.0\t1\t7\t0.28043205411160976\t1.2111948741909506\t1.0\n+16653\tgroup3\t3\t1\t0.7949225542667253\t1.0857032601068828\t3.0\t6\t2\t1.0650583476769155\t0.3033123505712213\t2.0\n+6630\tgroup1\t3\t1\t1.052039376110604\t1.0196469054881838\t3.0\t0\t1\t0.9599728103261282\t0.27691535289423685\t2.0\n+7997\tgroup3\t2\t3\t0.9373885579676389\t1.0929122783305962\t3.0\t0\t3\t0.4138438514343339\t1.0972371520000848\t1.0\n+16842\tgroup2\t2\t1\t1.0123686877389662\t1.0138181807538984\t3.0\t4\t1\t1.0386558479472863\t0.29642211452153533\t2.0\n+15418\tgroup3\t0\t0\t0.4977397972254476\t0.29342809208998816\t0.0\t0\t3\t0.43103048686992385\t1.0884920273765328\t1.0\n+16977\tgroup3\t4\t3\t0.9812798602404847\t1.1004306057391895\t3.0\t2\t5\t0.44708787541087025\t1.052007688880996\t1.0\n+8249\tgroup3\t0\t0\t0.24398068905795048\t0.4446066052836349\t0.0\t3\t1\t1.0972131920857282\t0.5432789321063332\t2.0\n+38624\tgroup2\t4\t5\t0.914525673837027\t1.0025201921010993\t3.0\t3\t9\t0.6126885981267824\t0.9966513382000214\t1.0\n+13713\tgroup3\t0\t0\t0.3508954425341463\t0.14396342558201491\t0.0\t1\t0\t0.5455349699532419\t1.0157638758893444\t1.0\n+14578\tgroup2\t0\t0\t0.3855539972318272\t0.27075459163213833\t0.0\t1\t0\t0.9670287407119521\t0.3666579233697094\t2.0\n+12152\tgroup3\t0\t0\t0.2933301214987318\t0.23591015760260986\t0.0\t3\t2\t0.8573550195271602\t0.4908410106271917\t2.0\n+28796\tgroup1\t0\t2\t0.39605190921232225\t0.31741453320872637\t0.0\t1\t0\t0.3478635238179973\t0.8191960715673223\t1.0\n+15427\tgroup3\t2\t1\t0.9733283804040067\t1.0062014237283545\t3.0\t0\t1\t0.3537778336845988\t0.9286706453933524\t1.0\n+16065\tgroup3\t1\t2\t0.9379392206019438\t1.0045767370103011\t3.0\t0\t5\t0.12951255367862535\t0.8870814948670983\t1.0\n+15383\tgroup2\t1\t4\t1.009630850730726\t1.0807864197771637\t3.0\t2\t1\t0.34921474837059063\t1.0488463942447201\t1.0\n'
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_1gauss.npz
b
Binary file test-data/small_1gauss.npz has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_1gauss.png
b
Binary file test-data/small_1gauss.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss.npz
b
Binary file test-data/small_2gauss.npz has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss.png
b
Binary file test-data/small_2gauss.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_convergence.png
b
Binary file test-data/small_2gauss_convergence.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_corner.png
b
Binary file test-data/small_2gauss_corner.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_individuals.png
b
Binary file test-data/small_2gauss_individuals.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_means.txt.gz
b
Binary file test-data/small_2gauss_means.txt.gz has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_neff.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/small_2gauss_neff.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,1 @@
+137.99048317720332
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_p.png
b
Binary file test-data/small_2gauss_p.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_p.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/small_2gauss_p.txt Mon Oct 02 13:24:52 2023 +0000
b
@@ -0,0 +1,6 @@
+name low median high
+mu0 0.1678923686281538 0.3241302687316502 0.39107634005727987
+scale0 0.15226068811267196 0.24583330904194123 0.37859276190581687
+amp1 0.4386233342029042 0.5270427979795759 0.6088863690713495
+mu1 0.9164251678140255 0.964314561582573 1.0028530401005753
+scale1 0.1297819189381704 0.1469905372180511 0.17782311769651776
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_pdf.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/small_2gauss_pdf.txt Mon Oct 02 13:24:52 2023 +0000
b
b'@@ -0,0 +1,101 @@\n+x\tlow\tmean\thigh\tmedian\n+0.0125\t0.12670356263254787\t0.376720679973249\t0.6437415260049979\t0.3317127733119675\n+0.037500000000000006\t0.17582330089312104\t0.4133114389760649\t0.6657809415891787\t0.380153005890229\n+0.0625\t0.23661521822386272\t0.4540541316911476\t0.6832995626346544\t0.43691685880206693\n+0.08750000000000001\t0.30424970353426534\t0.4988702388572368\t0.6983802820219329\t0.4915631177662837\n+0.1125\t0.3836397006423265\t0.5473445748809775\t0.7134332346822652\t0.5461082250352205\n+0.1375\t0.4708266792551106\t0.5986391499219265\t0.7285661973612422\t0.6012880766817457\n+0.16250000000000003\t0.5551375505802325\t0.6514384431366499\t0.7470742094200006\t0.653654387119493\n+0.1875\t0.61502546702653\t0.703948288081012\t0.790363360632909\t0.696695268582255\n+0.21250000000000002\t0.6555642774004675\t0.7539661666433396\t0.8733336733584552\t0.7322890417604377\n+0.2375\t0.6736783299242419\t0.7990301467677411\t0.9553680856288194\t0.7704497049572413\n+0.2625\t0.675393329085913\t0.8366382866905671\t1.03127744143166\t0.7982012185867823\n+0.28750000000000003\t0.6698533344840528\t0.8645133836334012\t1.0996024965966875\t0.823522553234274\n+0.3125\t0.6562280552101921\t0.8808740947486529\t1.142785582751362\t0.8392720843559668\n+0.3375\t0.6382151359869952\t0.8846672294323755\t1.1604043250838079\t0.8466083079757122\n+0.36250000000000004\t0.6173377618864876\t0.875720186981675\t1.1560793721936236\t0.8443153762859131\n+0.3875\t0.5914650262833041\t0.8547869133094897\t1.1240052657432054\t0.8313284636476609\n+0.41250000000000003\t0.5705406062389853\t0.8234820162260114\t1.080570410392243\t0.810222608338741\n+0.4375\t0.5419078747864474\t0.7841202275217672\t1.0128961100558458\t0.7766448456894333\n+0.4625\t0.5138324919487681\t0.7394962235359692\t0.9540304329968603\t0.7293416737260087\n+0.48750000000000004\t0.48571534158557045\t0.6926484148940208\t0.9072009344099945\t0.676063588158053\n+0.5125\t0.45391290335369533\t0.6466480118722264\t0.8523350032740965\t0.6331960635576673\n+0.5375000000000001\t0.41047400944220847\t0.6044428406185313\t0.7965606263035216\t0.5993649676518484\n+0.5625\t0.38123148840420584\t0.5687676573306188\t0.7501323052933836\t0.5690424621949819\n+0.5875000000000001\t0.3543061438080897\t0.542113603462916\t0.7264391975547957\t0.5441105721493111\n+0.6125\t0.32724069597173433\t0.5267330905848387\t0.7173242011337101\t0.5288342316799834\n+0.6375\t0.31804086031047735\t0.5246457981204454\t0.7245401188480379\t0.5273719067243517\n+0.6625000000000001\t0.31624006682928635\t0.5376084156812828\t0.742884886754941\t0.5432558741428829\n+0.6875\t0.3310007734440721\t0.5670161897604825\t0.7818727108791022\t0.5700365662085624\n+0.7125000000000001\t0.36940569539834234\t0.6137185131089978\t0.8439623603678892\t0.6116536275744638\n+0.7375\t0.41677681347413686\t0.6777530885263192\t0.9269493856113262\t0.6743122313800531\n+0.7625\t0.4913668702340747\t0.7580313025405345\t1.0276580053619266\t0.7478932141250783\n+0.7875000000000001\t0.5874265189863948\t0.8520365338113601\t1.1247782980233165\t0.8475468262671909\n+0.8125\t0.685810330771705\t0.9556196599591928\t1.2172606247371305\t0.9502169711633773\n+0.8375000000000001\t0.794304005115456\t1.0629829523318268\t1.3154534880295068\t1.0657040275966423\n+0.8625\t0.9165123222702052\t1.1669270429180645\t1.402016413085469\t1.1699316029316447\n+0.8875\t1.0305435760041086\t1.2593928418224503\t1.4900394994454742\t1.2589406428182985\n+0.9125000000000001\t1.1297900139891905\t1.332266397229696\t1.5469073476436306\t1.32263649268865\n+0.9375\t1.1750831650595561\t1.3783437991977383\t1.5895895283441637\t1.3672328869915744\n+0.9625000000000001\t1.1946000198640179\t1.3922954500262892\t1.6033519094265078\t1.3863729202396793\n+0.9875\t1.173803461531464\t1.3714443596966461\t1.573530369418087\t1.3708376533507989\n+1.0125000000000002\t1.1238855540959782\t1.3161939614945208\t1.5181046504608666\t1.3238038057291357\n+1.0375\t1.0439943131278036\t1.2300062920571133\t1.4171311907374347\t1.2354465179558072\n+1.0625\t0.9450233085659544\t1.1189256752375434\t1.2924977187041686\t1.1191252696927967\n+1.0875\t0.8341013019187482\t0.990740559203444\t1.1580443234371118\t0.9824103106873007\n+1.1125\t0.7197142880443617\t0.8539495910420355\t0.99342078903853'..b'7460716629863\t0.010707922712163104\t0.0027450054760013182\n+1.5375\t0.00025613804688462617\t0.0037427375420525885\t0.00753116916763843\t0.0016013136665934324\n+1.5625\t0.00011765083122848543\t0.0026331109585060537\t0.005352484925194219\t0.0009287938926587001\n+1.5875000000000001\t5.345546572807357e-05\t0.0018797189830566151\t0.003706767651857873\t0.0005400365644494052\n+1.6125\t2.3370937768727924e-05\t0.0013633089408364761\t0.0026440118796045684\t0.00030864159562653503\n+1.6375000000000002\t9.757946635934111e-06\t0.0010049673142406117\t0.0018339507562112334\t0.0001746151500012073\n+1.6625\t4.039741112936451e-06\t0.0007526806894907278\t0.001292680054208138\t9.40478145056667e-05\n+1.6875\t1.6500717782002303e-06\t0.0005721885744573551\t0.0009137107861820259\t5.0820204871803656e-05\n+1.7125000000000001\t6.450152545630558e-07\t0.00044086907427803423\t0.0006708552319324961\t2.7717825334092087e-05\n+1.7375\t2.3846466162823973e-07\t0.00034370468767527806\t0.00047543161078548676\t1.46121737605902e-05\n+1.7625000000000002\t8.243420675795641e-08\t0.0002706441019363624\t0.0003278698630743087\t7.877429298834631e-06\n+1.7875\t2.9567149863899117e-08\t0.00021488642580696894\t0.00024168245961583184\t4.111864900975251e-06\n+1.8125\t9.916098951984631e-09\t0.00017176840016090353\t0.00017766484521680227\t2.0436068952489992e-06\n+1.8375000000000001\t3.2193240187408435e-09\t0.00013804293056988233\t0.00013143522703209913\t1.0600421853937708e-06\n+1.8625\t1.0117680854831239e-09\t0.00011141044642529393\t9.785280159960393e-05\t5.473072069593378e-07\n+1.8875000000000002\t3.0781495161163277e-10\t9.021323230157927e-05\t7.074636710704377e-05\t2.7667775698594076e-07\n+1.9125\t8.784410165459977e-11\t7.32347701720531e-05\t5.0896557913777374e-05\t1.4334140711404513e-07\n+1.9375\t2.4857260499140447e-11\t5.956684435104026e-05\t3.658911076926906e-05\t6.963557021501608e-08\n+1.9625000000000001\t6.662200736417528e-12\t4.852052343408162e-05\t2.6311131093125673e-05\t3.4044224732833354e-08\n+1.9875\t1.9056763152600633e-12\t3.9565711780163867e-05\t1.8886957595195303e-05\t1.6580410293477873e-08\n+2.0125\t4.928560749226767e-13\t3.228945069884334e-05\t1.352099797664341e-05\t7.419712778697193e-09\n+2.0375\t1.2339101599689579e-13\t2.6366651984427003e-05\t9.69445751013165e-06\t3.125964019126105e-09\n+2.0625\t3.2403241079413646e-14\t2.1539179820026018e-05\t7.06369595888429e-06\t1.410577788077171e-09\n+2.0875000000000004\t7.757118288584117e-15\t1.7600622124301562e-05\t5.036216375816057e-06\t6.86379142603198e-10\n+2.1125\t1.9271977144216315e-15\t1.4385003928280595e-05\t3.6376430646191413e-06\t3.2763629544560525e-10\n+2.1375\t4.1523692324292487e-16\t1.175828100593813e-05\t2.594323439969036e-06\t1.6090157616855693e-10\n+2.1625\t8.852963237642496e-17\t9.611830711030293e-06\t1.8184378303147285e-06\t7.84766400208787e-11\n+2.1875\t2.1555442747497686e-17\t7.857404013167623e-06\t1.3003637962897524e-06\t3.5599468027780326e-11\n+2.2125000000000004\t4.851942133719642e-18\t6.423165585622127e-06\t8.996097050192052e-07\t1.6160783796601712e-11\n+2.2375\t1.0194642589233542e-18\t5.250557505953479e-06\t6.435592766105672e-07\t6.816812839020974e-12\n+2.2625\t2.17333630688025e-19\t4.2917957238637915e-06\t4.300763413090321e-07\t2.9131748510317703e-12\n+2.2875\t4.107269839154287e-20\t3.5078590515533083e-06\t3.0031522258324975e-07\t1.1114672687770382e-12\n+2.3125\t8.387300721139591e-21\t2.866865809433882e-06\t2.146600618768069e-07\t5.142422579082654e-13\n+2.3375000000000004\t1.4540286557812869e-21\t2.3427584319071066e-06\t1.4845851974195882e-07\t2.0359990500994472e-13\n+2.3625\t2.7484084683373016e-22\t1.914234572951312e-06\t1.0321046150506166e-07\t7.967482449080433e-14\n+2.3875\t6.275110725947119e-23\t1.5638767001202602e-06\t7.203111322752142e-08\t3.189224431493489e-14\n+2.4125000000000005\t1.1451941049101382e-23\t1.2774422581206118e-06\t4.863788994201108e-08\t1.44482250041291e-14\n+2.4375\t2.7661805033454506e-24\t1.043284181264797e-06\t3.339333022316414e-08\t5.052616401472067e-15\n+2.4625000000000004\t5.705909720093128e-25\t8.518774930918978e-07\t2.257968257644151e-08\t1.9689455218905996e-15\n+2.4875\t8.832750390074374e-26\t6.954324041647874e-07\t1.5528422011773676e-08\t6.965344998268519e-16\n'
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_posterior_andco.png
b
Binary file test-data/small_2gauss_posterior_andco.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_posterior_individuals.png
b
Binary file test-data/small_2gauss_posterior_individuals.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_posterior_per_cell.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/small_2gauss_posterior_per_cell.txt Mon Oct 02 13:24:52 2023 +0000
b
b'@@ -0,0 +1,2362 @@\n+mu\tsd\n+0.33086857263360775\t0.247036503702648\n+0.5491477359946301\t0.28386170259352617\n+0.22467653306651206\t0.16938334730238958\n+0.4261111139938885\t0.2985330354803467\n+0.35012225039459355\t0.25907805436643583\n+0.5190546403018004\t0.3301770026801263\n+0.4181008608434141\t0.2949734570797291\n+0.37139014781670254\t0.27141752164420935\n+0.6634902375868398\t0.2988762749115825\n+0.9316004760046239\t0.20368830150373668\n+0.6912419775213312\t0.29895595016356136\n+0.6810604805580249\t0.299082058616062\n+0.28487181099410125\t0.21515376084048576\n+1.018508233223392\t0.1535046993564713\n+0.38639257821030937\t0.27951081132144695\n+0.8080105375675627\t0.2563078052200935\n+0.4744757791421052\t0.3171621570504075\n+1.0571575189460516\t0.15260288313698803\n+0.795848919163722\t0.2586011367080841\n+0.41189334166703034\t0.29211974235085175\n+0.7664210634124616\t0.2632227975032267\n+0.8087137587782266\t0.2561686242795911\n+0.23041243413517845\t0.1737538702448312\n+0.35509813601276213\t0.26205605431582735\n+0.7426853375518987\t0.2956353381388062\n+0.7561597843951798\t0.2645122651788845\n+0.5993300493864279\t0.29344389950041466\n+0.7209296373935035\t0.2975788393706798\n+0.3040784591461396\t0.22896184573957126\n+0.9019949688753013\t0.21201337818936575\n+0.7497504838038797\t0.26523052818617104\n+0.710694840844321\t0.2982217338539133\n+0.8234299988118662\t0.25309661721167287\n+0.8183020094462835\t0.2542011399249339\n+1.0802218174190048\t0.14816641128480504\n+0.49406731234618884\t0.3233579267785503\n+1.0228659051024958\t0.1608150211802086\n+0.8909255793132065\t0.21507874488738749\n+0.3781659546390094\t0.27513528271865223\n+0.8119428225571325\t0.2555204774341058\n+1.0384280905426386\t0.14658329545339127\n+0.9572485158822553\t0.17070075554997147\n+0.6462242746087918\t0.29814464871707297\n+1.030462438405125\t0.16973803908010238\n+0.26739840865070974\t0.2021007325766286\n+0.3897018856986143\t0.28122832962736855\n+0.42532787417496676\t0.2981910549773167\n+0.3744649788998871\t0.27311739811856206\n+0.6894738445618068\t0.2989906652637997\n+0.3659834156812915\t0.2683769533340012\n+0.7265155325938699\t0.29715425217688735\n+0.9327795708602975\t0.18816774129642894\n+0.6248920739918503\t0.29650380613775323\n+0.5670318367778188\t0.2878300879958361\n+0.658522008163567\t0.26780341413307435\n+0.3618966388373674\t0.26603505895295065\n+0.35184517609770877\t0.2601154783533435\n+0.3962837823814707\t0.28457196118824823\n+0.3897835337133333\t0.2812703964212488\n+0.6196628725592699\t0.29597563891237055\n+0.611410856474712\t0.2950402148232131\n+0.44770711534893665\t0.30745036536048626\n+0.7534434251783654\t0.2648249209014339\n+0.4046972039745766\t0.2887065732828538\n+0.4293608071023479\t0.2999379305335418\n+1.0127981827744166\t0.16174980915904885\n+0.9045059934319415\t0.1954873034806193\n+0.9931812506744123\t0.15594679950042586\n+1.0043804534355187\t0.16271831438605358\n+0.9145613496362927\t0.1928697738111152\n+0.9934384842583838\t0.1559155386839213\n+0.38834534281098715\t0.280527239033953\n+0.9408483843773181\t0.20111355223361968\n+0.3705412887264126\t0.2709445003855083\n+0.3744522843625436\t0.2731104237761506\n+0.9077963027732773\t0.210388958304018\n+1.0946505468638124\t0.13661215921825462\n+0.3206365112839805\t0.2403086404270786\n+0.5708531084816525\t0.28859779607877606\n+0.41466286888892395\t0.2934032686281066\n+0.8152181228350323\t0.2548479985888227\n+0.694165545417337\t0.26854320466918025\n+1.098409734079681\t0.14520569060611815\n+0.42532593630243704\t0.23867767405832788\n+0.5672347304315699\t0.28787155586507074\n+0.3030274678891571\t0.2282233679383151\n+0.6478776937696753\t0.2982376560599995\n+0.9673809515227353\t0.17982433392040795\n+0.8030475825817696\t0.2572698023661788\n+0.9267581574271272\t0.20504786950269568\n+0.7311307179257496\t0.2967644802967746\n+1.013879635587044\t0.1616374055661056\n+0.813749693392389\t0.2551513562286425\n+1.058374800847919\t0.1526378739867504\n+0.38940275729700624\t0.2810740856097636\n+0.6601987677597654\t0.29877753882979224\n+0.9849953007851809\t0.1761463179249889\n+0.7170682752239004\t0.29784198745163265\n+0.9196086192491736\t0.20706225923802618\n+0.6368028716882841\t0.29752119831131374\n+0.6361823235321716\t0'..b'79210878\n+0.9726041684149543\t0.15883031652496402\n+0.7485866376666537\t0.26535369140738396\n+0.5893541290359126\t0.29192055727883837\n+0.8201533945084601\t0.25380650218836487\n+0.931361763144829\t0.20375518371942564\n+0.6596714122065276\t0.2987599429839115\n+0.819805652715358\t0.253880984417774\n+0.3315009457861924\t0.24744493349972066\n+0.6649005476133472\t0.2989127406094553\n+0.4185737736552525\t0.295187444670044\n+0.5898641669409223\t0.3432146807033369\n+0.6200861347898955\t0.29602024893740936\n+0.6788261789430159\t0.2685061768981933\n+0.8498383232030872\t0.2468687115400429\n+0.35360725384982755\t0.26116960236481906\n+1.0688049777699111\t0.14730866477928733\n+0.34642805947772964\t0.256831314665136\n+0.7755714635747524\t0.26193037875438624\n+0.9143408658252499\t0.2085475610695636\n+1.1150304545476306\t0.14778552914932774\n+0.6141128648722997\t0.29536029373090916\n+0.7832871045705448\t0.26073813770769194\n+0.8139316105467876\t0.255113938444678\n+0.9959774183588527\t0.17413716342647517\n+0.8497173641243085\t0.24689917478091106\n+0.6473379021330388\t0.2982078288465926\n+0.8343839397221121\t0.2506198596247329\n+0.687350009670622\t0.2990252573802831\n+1.0148804579648527\t0.17131952746799442\n+0.9928573611061682\t0.17468284180309004\n+1.0215520164278806\t0.1533029609563672\n+0.9241864680294853\t0.2057718170709521\n+0.6424737247643977\t0.29791554650901664\n+0.3194397766728125\t0.2395072696107159\n+0.6559153438259062\t0.29862039324661727\n+0.37095279810320086\t0.2711740135796347\n+0.994306725665988\t0.1640831017894032\n+0.678055527875843\t0.2990850226201896\n+0.9262667692810892\t0.2051861180164325\n+0.8108909151336615\t0.25573324629750105\n+0.40187495309151433\t0.2873370055417824\n+0.6262213094287471\t0.2966301109153165\n+0.7539224920779654\t0.2427700524740976\n+0.6430979508076609\t0.29795542531250147\n+0.9356670744367184\t0.22231552218545475\n+0.4223557008483703\t0.29688138042729034\n+0.7280542724047185\t0.2671539274582725\n+0.3513480296616929\t0.2598168137622624\n+0.6050711870073908\t0.2942362819015585\n+0.6770074835360047\t0.2990823681964528\n+1.0506303768811633\t0.15248639253649765\n+1.0093846967431506\t0.16212281636916728\n+0.6564938824911337\t0.2986435148117375\n+0.4187112245922003\t0.2952495490789696\n+0.39341863157346674\t0.2831282633961746\n+0.32260144761239357\t0.24161790930976149\n+0.41606608293388453\t0.29404723002054295\n+0.7412917118214144\t0.2660744058510391\n+0.8606729330850884\t0.24407400682517524\n+0.822459279026873\t0.25330844205317876\n+0.657844761985329\t0.29869519524227045\n+0.4349016544080557\t0.30228153046448714\n+0.7482385659371712\t0.2950154298682092\n+0.5542594478805604\t0.2850594486968247\n+0.4102159205150578\t0.29133423674789766\n+0.2944406288751901\t0.2221132563281476\n+0.9515103477637135\t0.1835107202995593\n+0.8119708593952131\t0.2555147850127299\n+0.9519518688980295\t0.18340438917591312\n+0.6338414595696431\t0.2972922382927267\n+0.8627231868900813\t0.24353101791014117\n+1.0223911730033772\t0.15325100726394364\n+0.6471464222845764\t0.2981971232375375\n+0.4285845001708905\t0.29960436804831575\n+0.8709317221725978\t0.22043801554732004\n+0.4144221963890022\t0.29329239066870655\n+0.38804797962222776\t0.28037300784468716\n+0.7824267710485932\t0.26087566643587673\n+0.37357821014004583\t0.2726293397251401\n+0.41045277692182536\t0.2914455239048292\n+0.8664777634568681\t0.24252560859541486\n+0.9020000937382548\t0.15620148228250627\n+0.8026094954788531\t0.25735300182244997\n+0.6173640754457205\t0.29572762353172805\n+0.9648960699672491\t0.1803811969895106\n+0.6174822401387567\t0.2957406084504126\n+0.9049961910953458\t0.211174188292788\n+1.0016436668842221\t0.18643988232219938\n+0.9096241805950299\t0.23015377050622465\n+0.7915503769489123\t0.2593593404285475\n+0.38437681201141555\t0.27845267013370945\n+0.9705336742372228\t0.1791301648370778\n+0.5012969754506618\t0.32545426385142856\n+0.7676897473749488\t0.2237760415802802\n+0.408211340336399\t0.29038748252246527\n+0.39582247296661516\t0.2843407433929141\n+0.43216630601993405\t0.3011327056380514\n+0.2627276657824346\t0.1985523090655405\n+0.8115642537787877\t0.2555972308204999\n+0.6234118024240579\t0.29635936031905524\n+0.6350800609256917\t0.29738992630862804\n'
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/small_2gauss_with_posterior.png
b
Binary file test-data/small_2gauss_with_posterior.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/test1_res.png
b
Binary file test-data/test1_res.png has changed
b
diff -r 000000000000 -r 02b1fe7aed76 test-data/test2_res.pdf
b
Binary file test-data/test2_res.pdf has changed