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

astralProjection is not plotting the proportion of gene trees (f1/2/3) #2

Open
DoveOfGlass opened this issue Apr 30, 2024 · 0 comments

Comments

@DoveOfGlass
Copy link

Hi!!

Thank you for developing this package! :)

I am trying to plot my tree (from Astral III) using the following lines:

Using a premade tree:

outgroups = c("rr573")
astral.data = createAstralPlane(astral.tree = "path/mytree.tre",
outgroups = outgroups,
tip.length = 1)
quartetplots <- astralProjection(astral.plane = astral.data,
local.posterior = TRUE,
pie.plot = TRUE,
pie.data = "genetree",
pie.colors = c("blue", "purple", "green"),
node.color.text = c("black"),
node.color.bg = c("transparent"),
node.label.size = 1,
tip.label.size = 1,
pie.chart.size = 0.25)

It works, but If I change the argument "pie.data" to "pie.data = "qscore"", the piecharts that I obtain are the same (as the ones that I obtain with "pie.data = "genetree"") and correspond to the qscores.

I checked "astral.data" and it reads properly q1/2/3 and f1/2/3 as numbers.

The format of the labels of my tree is the following: "q1=0.3889582631081762;q2=0.257503503098498;q3=0.3535382337933258;f1=77.01373609541889;f2=50.985693613502605;f3=70.0005702910785;pp1=0.7825751429664773;pp2=0.03226467744511558;pp3=0.1851601795884253;QC=390;EN=198.0"

I don't see any warnings at any point of the process. I have tried erasing the q1/2/3 form the tree and I simply lose the piecharts, it does not use the f1/2/3. Any clues? Am I doing something incorrectly?

Thank you very much!!!
Paloma :)

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

1 participant