-
Notifications
You must be signed in to change notification settings - Fork 51
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
Error in seq_gtf(gtf, seqpath, ...) : unused argument (readlen = 70) #24
Comments
Hey @irenerodriguez - sorry for the delay here! Looking in to this and plan to push out a fix soon. Thanks for your patience. |
Hi @alyssafrazee , This is a great package, I'm very keen to try out those bias model you have included, unfortunately I've run into the same issue using version 1.7.1 installed via bioconductor, it looks like none of those additional arguments is working in simulate_experiment()? Best, |
I've also hit this bug after installing the package from bioconductor even though its almost 5 months from the bugfix. I guess this hasn't made it into the bioconductor release yet? |
Yikes, I'm so sorry about that -- I thought we had pushed to BioC, but it doesn't look like it yet. I'll look into it this week. In the meantime, installing from github directly should work. Thanks for your patience. |
@nickschurch Hey Nick, try installing again from the devel branch using BioC. Please let me know if it has since been resolved. Will also update the release branch if confirmed this issue is resolved
|
Thanks. JMF47: I tried this but it doesn't seem to have solved the problem. The version of polyester I have according to sessionInfo() is 1.6.0 - is that the version in the dev branch? |
Hi @alyssafrazee
First of all, thank you for writing this package.
I am trying to simulate some reads from GTF format + FASTA sequence files. I adapted the example provided in https://github.com/alyssafrazee/polyester --> simulating reads.
Everything works fine when I don't use the readlen parameter in the simulate_experiment function, and I get my simulations with read length equals 100 (default).
However, when I use this parameter in the simulate_experiment function (for example, readlen=70), I obtain the following error:
Error in seq_gtf(gtf, seqpath, ...) : unused argument (readlen = 70)
It works:
simulate_experiment(gtf='mygenes.gtf', seqpath=seqpath,
reads_per_transcript=readspertx, num_reps=c(2,2,2),
fold_changes=fold_changes, outdir='output_reads',idfield='name')
It doesn't work:
simulate_experiment(gtf='mygenes.gtf', seqpath=seqpath,
reads_per_transcript=readspertx, num_reps=c(2,2,2),
fold_changes=fold_changes, outdir='output_reads',idfield='name',readlen=70)
Thank you for your help.
Best
Irene.
The text was updated successfully, but these errors were encountered: