I am using
godon to find genes under positive selection with the branch-site test. However, after running for a little while I get the error "panic: division of zero by zero or infinity by infinity".
Below is the full output with the error report. I can supply the phylogenetic tree and fasta file I am using if that would be helpful. Thank you.
panic: division of zero by zero or infinity by infinity
goroutine 1 [running]:
math/big.(*Float).Quo(0xc4273a8f40, 0xc4273a8f68, 0xc4273a8f18, 0xc42e88c118)
bitbucket.org/Davydov/godon/cmodel.(*BranchSite).BEBPosterior(0xc422164120, 0xc4259bb900, 0x91, 0x91)
main.computeFinal(0x1f263a0, 0xc422164120, 0xc4201522a0, 0x101, 0x6a80, 0x0)
main.performSingleTest(0xc420162320, 0xc428266000, 0x9399, 0x0, 0x0, 0x0, 0x0, 0x0, 0x0, 0x0, …)
main.hypTest(0x1ee03a7, 0x1ee03a7, 0x4, 0xc4201aa501, 0x1, 0x1)
Ah wait you mentioned a branch sites test, right that makes sense. You've got two or more identical sequences in the analysis. Ensure there are no identical sequences in your alignments and ... PRESTO! It will work :-)
Below are my original thoughts before I re-read your question. I've kept the notes because it helps explain the theory a bit. There is no need to supply the tree, it should be resolved.
The tree likelihood looks to be decreasing nicely, so it can see both the tree and the codons in the alignment.
At a guess the error relates to two or more indentical nucleotide sequences in your alignment. Thus zero non-synonymous sites (amino acid change) and zero synonymous (silent) sites. This error would only occur if you were calculating selection for each 'node' in the tree (internal branch). Alternatively, possibly if the number of synonymous sites was zero, regardless of the non-synonymous mutations this error could appear, but this doesn't seem to fit the error output.
I don't know "Godon", but recommend
PAML for this analysis. Admittedly it is complicated to deploy, but it is a very robust algorithm and hugely respected.
Answered by M__ on April 2, 2021
2 Asked on August 21, 2021
1 Asked on August 21, 2021 by najeha-mohamed
1 Asked on August 21, 2021 by mariam
2 Asked on August 21, 2021 by felipe-almeida
1 Asked on August 21, 2021 by ace8910
1 Asked on August 19, 2021
4 Asked on August 19, 2021 by jameshfisher
2 Asked on August 19, 2021 by clove444
1 Asked on August 19, 2021 by ophelia-vallee
2 Asked on August 19, 2021
0 Asked on August 18, 2021 by sasam
0 Asked on August 17, 2021 by unicorn
1 Asked on August 15, 2021
1 Asked on August 13, 2021
Get help from others!