-
Notifications
You must be signed in to change notification settings - Fork 180
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
Fix star index issues #350
Conversation
|
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.
You can do it like this, but alternatively you can use the getGenomeAttribute function and assign the igenomes star index directly to params.star_index
like it is done here. This way, you could leave line 70 as it is and params.star_index
will always have the correct content.
My proposed approach does not overwrite params.star_index
with an igenomes index if it has already been set manually, because nextflow does not allow overwriting parameters.
Will leave this up to you
That's great, wasn't aware of that tooling. I should probably do the same for gtf/fasta then. |
star_index = params.star ? file(params.star_index, checkIfExists: true) : ( params.genome ? file( getGenomeAttribute('star'), checkIfExists: true ) : [] ) | ||
ch_star_index = [[id: star_index.baseName], star_index] |
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 it like this:
params.star_index = getGenomeAttribute('star')
ch_star_index = [[id: star_index.baseName], file(params.star_index, checkIfExists: true)]
You are utilizing params.star
which I have not seen elsewhere in the pipeline
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.
dammit, good catch
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'll fix it in a separate PR
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.
params.star_index = getGenomeAttribute('star')
But that logic seems to overwrite params.star_index
with the one from iGenomes even when star_index
is explicitly set. I find that a bit strange... or am I overlooking something here?
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.
ah, just saw your edit above now
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.
That simplifies things
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.
💯
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.
Feel free to request my review on a new PR
Close #277
Close #327
PR checklist
nf-core 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).