STITCH
STITCH copied to clipboard
how to produce input file
Hi, how to produce input file the input file gen.txt ? I am use tped format file generated by plink software to input file. However, I found that each SNP locus has two genotypes, which is not quite the same as the gen.txt format in the example. I hope you can give me some suggestions on how to generate gen.txt . Thank you very much !
Hi, can you post an example of what your file currently gives you? I suspect you are getting phased output from tped, whereas gen.txt is not phased but the unphased genotype, so you would add together phased output in the obvious way to get unphased genotype. I can offer suggestions on this if you post part of an example. Note that gen.txt is optional, you don’t need it to run STITCH, it is just helpful to show accuracy as you progress through the EM algorithm