Description
Hi,
I am trying to use this package to plot quartet values on a species tree. I have done everything step-by-step exactly as it is in the demo, but I run into an error that I can't solve.
After setting my working directory, here's the code I've run successfully so far:
"setwd(work.dir)
tree.dir = "trees"
outgroups <- c("Ganc-1_merged_v2","Golym-1_merged_v2","Gyro-1_merged_v2","Gyro-3_merged_v2","Hydnot-2_S31_L003_v2","Hydnot-3_S19_L003_v2","Hydnot-4_merged_v2")
save.name = "morch_v2_BUSCO_pruned"
setupAstral(genetree.folder = tree.dir,
output.name = save.name,
min.n.samples = 18,
min.sample.prop = 0.5,
make.polytomy = TRUE,
polytomy.limit = 10)
runAstral(input.genetrees = "morch_v2_BUSCO_pruned_genetrees.tre",
output.name = save.name,
astral.path = astral.path,
astral.t = 4,
quiet = FALSE,
load.tree = FALSE,
multi.thread = TRUE,
memory = "1g")
astral.data = createAstralPlane(astral.tree = save.name,
outgroups = outgroups,
tip.length = 1)"
I get this error:
Error in names(x) <- value :
'names' attribute [11] must be the same length as the vector [3]
Any guidance or suggestions? Thanks