DiaNN icon indicating copy to clipboard operation
DiaNN copied to clipboard

Convert speclib to tsv format

Open Niveda-S opened this issue 2 years ago • 3 comments

Hi,

I am interested in doing a comparison between two pre-existing libraries. Is there a way to convert library in the speclib format to a readable format such as tsv?

Niveda-S avatar Sep 23 '22 16:09 Niveda-S

Yes, just specify an output library in the GUI

vdemichev avatar Sep 29 '22 12:09 vdemichev

Okay thanks! Is there a way to directly read the speclib format also (with the help of some sort of a database)?

Niveda-S avatar Sep 30 '22 17:09 Niveda-S

No, the only way is to convert with DIA-NN

vdemichev avatar Sep 30 '22 17:09 vdemichev

Dear @vdemichev

Although I specify the "\report-lib.tsv" in the GUI, and the options seem to be passed to the diann.exe

diann.exe --lib "" --threads 16 --verbose 1 --out "D:\WEW_o29132_\BENCHMARKING\report.tsv" --qvalue 0.01 --matrices  --out-lib "D:\WEW_o29132_\BENCHMARKING\report-lib.tsv" --gen-spec-lib --predictor --prosit --fasta "D:\WEW_o29132_\BENCHMARKING\UP000000625_83333.fasta" --fasta-search --min-fr-mz 200 --max-fr-mz 1800 --met-excision --cut K*,R* --missed-cleavages 1 --min-pep-len 7 --max-pep-len 30 --min-pr-mz 300 --max-pr-mz 1800 --min-pr-charge 1 --max-pr-charge 4 --unimod4 --use-quant --reanalyse --relaxed-prot-inf --smart-profiling --peak-center --no-ifs-removal 

neither the report.tsv file nor the report-lib.tsv is being written by DIA-NN (see log below).

0 files will be processed
[0:00] Loading FASTA D:\WEW_o29132_\BENCHMARKING\UP000000625_83333.fasta
[0:00] Processing FASTA
[0:02] Assembling elution groups
[0:03] 498603 precursors generated
[0:04] Prosit input saved to D:\WEW_o29132_\BENCHMARKING\report-lib.prosit.csv
[0:04] Gene names missing for some isoforms
[0:04] Library contains 4392 proteins, and 4392 genes
[0:04] Encoding peptides for spectra and RTs prediction
[0:05] Predicting spectra and IMs
[5:27] Predicting RTs
[5:56] Decoding predicted spectra and IMs
[5:58] Decoding RTs
[5:58] Saving the library to D:\WEW_o29132_\BENCHMARKING\report-lib.predicted.speclib
[6:00] Initialising library

[6:00] Log saved to D:\WEW_o29132_\BENCHMARKING\report-lib.log.txt
Finished

wolski avatar Jun 16 '23 20:06 wolski

Same problem here, the report-lib.tsv is not being written: diann.exe --lib "MYPATHWITHOUTSPACES\report-lib.predicted.speclib" --threads 28 --verbose 1 --out "MYPATHWITHOUTSPACES\report.tsv" --qvalue 0.01 --matrices --out-lib "MYPATHWITHOUTSPACES\report-lib.tsv" --gen-spec-lib --predictor --reanalyse --relaxed-prot-inf --smart-profiling --peak-center --no-ifs-removal I have tried 1.8.1 and 1.8.2 (beta) but got the same result.

Arthfael avatar Sep 20 '23 07:09 Arthfael

@Arthfael

Running the command

diann.exe --lib "D:\WEW_o29132_\BENCHMARKING\report-lib.predicted.speclib" --threads 16 --verbose 1 --out "D:\WEW_o29132_\BENCHMARKING\report.tsv" --qvalue 0.01 --matrices  --out-lib "D:\WEW_o29132_\BENCHMARKING\report-lib.tsv" --gen-spec-lib --reanalyse --relaxed-prot-inf --smart-profiling --peak-center --no-ifs-removal

On the generated speclib file generates then the report-lib.tsv file. (Some of the options I guess are not needed).

wolski avatar Sep 20 '23 07:09 wolski

It worked, thanks for your swift reply! So the "--predictor" flag was the issue.

On Wed, Sep 20, 2023 at 9:22 AM Witold Wolski @.***> wrote:

@Arthfael https://github.com/Arthfael

Running the command

diann.exe --lib "D:\WEW_o29132_\BENCHMARKING\report-lib.predicted.speclib" --threads 16 --verbose 1 --out "D:\WEW_o29132_\BENCHMARKING\report.tsv" --qvalue 0.01 --matrices --out-lib "D:\WEW_o29132_\BENCHMARKING\report-lib.tsv" --gen-spec-lib --reanalyse --relaxed-prot-inf --smart-profiling --peak-center --no-ifs-removal

On the generated speclib file generates then the report-lib.tsv file. (Some of the options I guess are not needed).

— Reply to this email directly, view it on GitHub https://github.com/vdemichev/DiaNN/issues/510#issuecomment-1727116950, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABSHV4PMQHBCWRR6CYAHB4LX3KKUFANCNFSM6AAAAAAQUEXQPI . You are receiving this because you were mentioned.Message ID: @.***>

Arthfael avatar Sep 20 '23 07:09 Arthfael