Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

SweeD does not report correct CLR for folded spectra #12

Open
oneoverx opened this issue Aug 27, 2020 · 0 comments
Open

SweeD does not report correct CLR for folded spectra #12

oneoverx opened this issue Aug 27, 2020 · 0 comments

Comments

@oneoverx
Copy link

Hello,
SweeD does not output a correct CLR (and probably alpha) for folded spectra as in the following example of a simulated strong sweep:

Output SweeD:
24954877.3877 0.000000e+00 2.660754e-02 24953940.0000 24955327.0000
24964879.3125 0.000000e+00 1.318682e-01 24964668.0000 24964971.0000
24974881.2373 0.000000e+00 1.008404e-01 24974129.0000 24974999.0000
24984883.1621 0.000000e+00 1.188119e-01 24984783.0000 24985169.0000
24994885.0869 0.000000e+00 6.417114e-02 24994602.0000 24995071.0000
25004887.0117 0.000000e+00 2.400000e+00 25004883.0000 25005130.0000
25014888.9365 0.000000e+00 1.008404e-01 25014769.0000 25015216.0000
25024890.8613 0.000000e+00 3.000000e+00 25024888.0000 25025193.0000
25034892.7861 0.000000e+00 8.695655e-02 25034568.0000 25035030.0000
25044894.7109 0.000000e+00 1.176471e-01 25044794.0000 25045055.0000

Output Sweepfinder:
24954877.725339 1.182454 3.270155e-04
24964879.650005 0.093961 7.456625e-03
24974881.574671 0.829999 8.739903e-04
24984883.499337 46.348489 3.772824e-05
24994885.424003 57.058797 2.168375e-05
25004887.348669 0.259692 2.800138e-02
25014889.273335 7.959813 2.066296e-04
25024891.198001 0.343199 3.594658e-02
25034893.122667 0.624442 8.493904e-04
25044895.047334 0.187886 9.855205e-03

The input file was in the Sweepfinder format, but 'ms' format yields the same problem. I tried version 3.2.1 and 4.0.0 of SweeD.

If the same input data is characterized as unfolded, (a zero instead of a 1 in the last column) both tools yield very similar results. Applying command line option "-folded" instead of marking the status in the file makes no difference.

I guuess this is a bug and not a feature...

Alex

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant