search

poseidon-framework / poseidon-analysis-hs

A tool to analyse genotype data (optionally in the poseidon data format)
MIT License
3 stars 0 forks source link

Xerses option for jackknife output before mean and average #18

Closed giulioisac closed 1 year ago

giulioisac commented 1 year ago

Hello,

It would be great to have the option in the fstat command to output the value of the statistic for each block instead of mean and variance from jackknife. A possibility would be to also output the beginning and the end of the relative positions in the genomes for each block of SNPs, so that values of the F2 can be compared across different samples. Is this already possible in some way? I cannot find it in the documentation.

Thanks! Best, Giulio

stschiff commented 1 year ago

Yes, I agree. The other command ras can already do this. I can implement for Xerxes as well. Thanks for the feedback

stschiff commented 1 year ago

Done in the latest version 0.3.2.0. Note that the per-block output is not in every sense the "raw" data before the Jackknifing operation. I can explain in more detail in person, and I'm also working on a little Math writeup to add to this repo.

giulioisac commented 1 year ago

That will help, thanks!