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

Insufficient number of observations #16

Open
srbehera11 opened this issue Oct 17, 2017 · 1 comment
Open

Insufficient number of observations #16

srbehera11 opened this issue Oct 17, 2017 · 1 comment

Comments

@srbehera11
Copy link

Can you please suggest what values should be manually assigned for resolving the following error?

[17/10/2017 11:26:33] Sorting splice-graphs by read count
[17/10/2017 11:26:33] Finished sorting splice-graphs by read count

[17/10/2017 11:26:33] Spawning 1 thread(s) for fetching alignments and 1 i/o thread
[17/10/2017 11:26:33] Estimating fragment length distribution from 4 transcripts longer than 2500 nucleotides
ERROR: Insufficient number of observations for estimating the fragment length distribution parameters. Please specify these manually using the --frag-mean and --frag-sd options.

@lassemaretty
Copy link
Contributor

Hi,

Thank you for reporting the problem. Bayesembler requires at least 1000 paired-end reads that align to single-isoform genes for estimating the length distribution of the sequencing fragments. You can try any of the following:

  • If you are only running on a subset of your data, you should try running on the full set
  • If you already have an idea about the mean and standard deviation of the fragment length distribution you can provide these as input using the options suggested in the error message
  • If you have multiple samples that you think should have the same fragment length distribution, you could try to pool the data across samples and do a prerun of Bayesembler on the pooled data to get the estimates from the Bayesembler output and then provide these as input to runs on individual samples as described above.

I hope this helps! Best regards,

Lasse

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

2 participants