search

stschiff / msmc2

GNU General Public License v3.0
53 stars 9 forks source link

How many samples should be used in the analysis #61

Open enluo211 opened 6 months ago

enluo211 commented 6 months ago

When I used 16 haplotypes for analysis(two poputation,four samples per population), the number of bits in my MSMC input file was very small, only a few hundred, when I used 8 haplotypes for analysis, the number of bits in the MSMC input file was several thousand (about 10 times the number of 16 haplotypes).Will too few sites in the input file affect the result? If so, how can I avoid this from happening? chr1.8hapmap.multihetsep.txt chr1.16hapmap.multihetsep.txt

stschiff commented 6 months ago

Well, your files have extremely different lengths. One has a few dozen sites, the other many more. You seem to be missing a hole bunch of sites from your 16-haplotype input file.

enluo211 commented 6 months ago

Yeah,but I don't know what causes this phenomenon

enluo211 commented 6 months ago

Is it because most loci are masked, and if so, why does this happen?

stschiff commented 6 months ago

Well, but how can I possibly know what went wrong with the creation of your file?

enluo211 commented 6 months ago

You are right, but all my operations follow the process provided by you on github (refer to msmc2 and msmc-tools).