bias correction for help

[Sithions Brand]

Hi everyone,

I encountered a problem while performing bias correction for SNP sequences generated by GBS-SNP-CROP. When I ran RAxML, I received the following error message: 'For partition No Name Provided, you specified that the likelihood score should be corrected for invariant sites via an ascertainment bias correction. However, some sites in this partition are already invariant. This is not allowed, please remove all invariant sites and try again, exiting...' I cannot understand why this error is occurring. My command was: ./raxmlHPC-SSE3 -s 1.phy -f a -x 1234 -p 1234 -# 1000 -m ASC_GTRCAT -V --K80 --asc-corr=lewis -n 3.phy -o S19,S20,S21,S18.

Thank you!

Xu Jun

[Alexandros Stamatakis]

Dear Xu,

Which part of the error message which I believe is rather verbose don't
you understand?

It tells you that your MSA already contains some invariant columns which
should be removed prior to doing an analysis under the ascertainment
bias correction.

The reason for this is that the maths do not work and break down if you
are trying to correct for asc. bias in a MSA that contains some
invariant sites already.

I believe that somebody wrote a script to remove such sites for exactly
this purpose (using asc. bias correction) a couple of years ago.

Please also consider switching to the new version of RAxML:

https://github.com/amkozlov/raxml-ng

Alexis

> --
> You received this message because you are subscribed to the Google
> Groups "raxml" group.
> To unsubscribe from this group and stop receiving emails from it, send
> an email to raxml+un...@googlegroups.com
> <mailto:raxml+un...@googlegroups.com>.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/raxml/2aa8d2aa-ec44-4cd6-ab7a-fc82c6485253n%40googlegroups.com <https://groups.google.com/d/msgid/raxml/2aa8d2aa-ec44-4cd6-ab7a-fc82c6485253n%40googlegroups.com?utm_medium=email&utm_source=footer>.

--
Alexandros (Alexis) Stamatakis

ERA Chair, Institute of Computer Science, Foundation for Research and
Technology - Hellas
Research Group Leader, Heidelberg Institute for Theoretical Studies
Full Professor, Dept. of Informatics, Karlsruhe Institute of Technology

www.biocomp.gr (Crete lab)
www.exelixis-lab.org (Heidelberg lab)

[Sithions Brand]

Dear Alexandros Stamatakis，

I found a script for deleting invariant sites (https://github.com/btmartin721/raxml_ascbias), but it doesn't seem to achieve my goal because when I run it again, I encounter the same error, and it appears that the invariant sites are still present. 

I also need to remove invariant sites as input for raxml-ng, and my input file consists of concatenated SNPs.

Thank you！

Xu Jun

[Alexis]

Maybe the script is not removing sites containing a single character but also some gaps and undetermined characters, such sites are also considered as being invariable.