DongmeiDai

Results 15 comments of DongmeiDai

Thank you for your reply. Is the 2 of g=cNe2 the square of the effective population size? If the g=57667,Ne=30000, so the c is 0.000064, is it right? Thank you...

Hello, If I want to calculate β2 using an outgroup to call substitutions and ancestral alleles, can I just get epo files from outgroup vcf file, without common ancestor? Thank...

Thank you. I will try it ! If I use vcfm2acf to convert VCF to ACF , then I use usepopsrootanc to tell the root(outgroup), finally I use the replaceanc,...

Hello, I have another question. Is the file2.acf.gz in https://github.com/grenaud/glactools the result using the glactools with vcfm2acf --onlyGT --fai parameter? I am not sure if using these parameters is correct?...

Hello, the file1.acf.gz only have a outgroup individual: ![image](https://github.com/grenaud/glactools/assets/76728696/61930e93-0364-4bde-887a-4d8c2c64f8ce) when I use the "usepopsrootanc" parameter, I only tell the outgroup, then the error is as follows ![image](https://github.com/grenaud/glactools/assets/76728696/8b1340ad-87fa-403e-85d7-02f969f87dd2) I think it...

Hello, What does "ChimpHumanAncestor " refer to specifically? Thank you very much

Hello, Thank you for your reply! 1. May I ask if I have only root and ancestor, but no groups like other apes, is it OK to generate the follow...

I think I misunderstood your meaning. Is file4.acf.gz in your example (glactools replaceanc file2.acf.gz file3.acf.gz > file4.acf.gz) equal to the 1kg_chr6.AFRsubset.acf.gz in following code ( glactools vcfm2acf --onlyGT --epo all.epo.gz...

sorry, 1. glactools replaceanc file2.acf.gz file4.acf.gz 2. glactools vcfm2acf --onlyGT --epo all.epo.gz --fai human_g1k_v37.fasta.fai chr6.AFRsubset.recode.vcf > 1kg_chr6.AFRsubset.acf.gz What I understand is that if there is no epo file, file4 generated...

Yes, I used 3 commands to generate the file successfully. Since my species is not human, I want to know if the output file of 1 is the equivalent of...