This directory contains the Mathematica script used to compare the joint deconvolution algorithm to the sequential one.
The script compareScript.wls is compatible with Mathematica v11. Due to a complete exploration of possible hyper-parameter values this computation is quite long (about 4 hours). Hence the saved result is also available in the comparisonResult.wl file.
Using this saved file, the plots can be regenerated by:
Get["./ComparisonSeqJoin/comparisonResult.wl"];
FormatGlobalResult[comparisonResult["resultSeqcNeg"],comparisonResult["resultJoincNeg"]] (* convex baseline *)
FormatGlobalResult[comparisonResult["resultSeqcPos"],comparisonResult["resultJoincPos"]] (* concave baseline *)The true peaks (deconvolved and convolved with the Gaussian peak shape) are:
