-
Notifications
You must be signed in to change notification settings - Fork 12
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
Minor release: 1.1.0 with CNV benchmark subworkflow #164
Conversation
typo in README.md
add zenado id
…ols-vcfeval implement rtgtools for somatic small bench
161 cnv benchmark subworkflow
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
LGTM!
I left a couple of comments/questions - but they are not blockers, so I leave it approved :)
bin/merge_reports.py
Outdated
@@ -191,6 +191,29 @@ def get_happy_resuls(file_paths): | |||
|
|||
return merged_df | |||
|
|||
def get_intersect_resuls(file_paths): |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
def get_intersect_resuls(file_paths): | |
def get_intersect_results(file_paths): |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Or is this meant to be called resuls ? I see you have it throughout the script so maybe I am missing something
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this is a sytex mistake :D I realize it but leave it as it is now
truth_id = "SEQC2" | ||
truth_vcf = "${params.test_data_base}/testdata/hg38/truth/somatic/cnv_benchmark_calls.vcf" | ||
regions_bed = "${params.test_data_base}/testdata/hg38/truth/somatic/ngs_benchmark_cnv_gain_loss_loh.bed" | ||
rename_chr = "${projectDir}/assets/rename_contigs/grch37_grch38.txt" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Should this data be in nf-core/testdatasets? could it be useful for other modules?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
It is in test-database, I am linking it through params.test_data_base
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'm confused by this comment
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I meant the last row - it is in projectDir not test_data_base, so I was wondering whether there was a reason for that
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
For some reason I highlighted everything :D
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
test_data_base = 'https://raw.githubusercontent.com/nf-core/test-datasets/variantbenchmarking'
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
it is already in test-datasets, that is what I meant :D
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this one too? "${projectDir}/assets/rename_contigs/grch37_grch38.txt"
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this is almost default for the pipeline itself, It is not test purposes.
@@ -0,0 +1,54 @@ | |||
process BEDTOOLS_INTERSECT { |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Is there no way to use https://nf-co.re/modules/bedtools_intersect ?
Or a patch to it? Or is it doing something different?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The python script I have, is with more functionalities like checking the proper file formats, and summarizing statistics.
workflows/variantbenchmarking.nf
Outdated
}else{ | ||
log.error "Please specify params.truth_id and params.truth_vcf to perform benchmarking analysis" | ||
log.error "Please specify params.truth_id and params.truth_vcf or regions_bed to perform benchmarking analysis" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Could this not be catched with nf-schema?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Or is it to better error messages? :)
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Might be possible, but I couldnt figure out how..
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
this is not being validated with nf-schema by the way. This is coming from nextflow_schema.json.
CHANGELOG.md
Outdated
@@ -3,14 +3,37 @@ | |||
The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/) | |||
and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html). | |||
|
|||
## v1.0dev - [date] | |||
## 1.1.0 - [06.03.2025] |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The date is of yesterday, but I guess it is not a big deal :) In case you decide to do a new PR t could be updated
apply luisas review
fix download_pipeline.yml
PR checklist
nf-core pipelines lint
).nextflow run . -profile test,docker --outdir <OUTDIR>
).nextflow run . -profile debug,test,docker --outdir <OUTDIR>
).docs/usage.md
is updated.docs/output.md
is updated.CHANGELOG.md
is updated.README.md
is updated (including new tool citations and authors/contributors).