msmc2 icon indicating copy to clipboard operation
msmc2 copied to clipboard

How many samples should be used in the analysis

Open enluo211 opened this issue 1 year ago • 5 comments
trafficstars

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

enluo211 avatar Apr 22 '24 02:04 enluo211

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.

stschiff avatar Apr 23 '24 15:04 stschiff

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

enluo211 avatar Apr 23 '24 15:04 enluo211

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

enluo211 avatar Apr 23 '24 15:04 enluo211

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

stschiff avatar Apr 23 '24 15:04 stschiff

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

enluo211 avatar Apr 23 '24 15:04 enluo211