Hello,
Our group is trying to calculate delta statistics and p-values across ~16,000 gene trees. This is taking a long time to compute (~1 min/tree), which ends up taking a long time at the scale of our data, and we're looking for ways to make this computation faster.
It looks like the sim parameter of the delta() function is related to the number of MCMC simulations. If we reduce sim from 10000 to 100, the code runs faster, but we're not sure whether this affects the correctness of the delta statistic. Could you please let us know how this would affect our analysis?
We saw that your analysis is on the same scale (~14,000 gene trees), so we are wondering how you were able to accomplish this. Any advice would be greatly appreciated.
Thank you!
Heather
Hello,
Our group is trying to calculate delta statistics and p-values across ~16,000 gene trees. This is taking a long time to compute (~1 min/tree), which ends up taking a long time at the scale of our data, and we're looking for ways to make this computation faster.
It looks like the
simparameter of thedelta()function is related to the number of MCMC simulations. If we reducesimfrom10000to100, the code runs faster, but we're not sure whether this affects the correctness of the delta statistic. Could you please let us know how this would affect our analysis?We saw that your analysis is on the same scale (~14,000 gene trees), so we are wondering how you were able to accomplish this. Any advice would be greatly appreciated.
Thank you!
Heather