-
Notifications
You must be signed in to change notification settings - Fork 28
Description
Dear authors
I use rfmix recently, while the result is always 1, that cannot be true.
~/software/rfmix/rfmix -f na.vcf.gz -r hg.vcf.gz -m pop.list -g allchrs.b37.gmap --chromosome=22 -o test
and resuts are below:
#reference_panel_population: HG
chromosome physical_position genetic_position genetic_marker_index NA19152:::hap1:::HG NA19152:::hap2:::HG NA18542:::hap1:::HG NA18542:::hap2:::H
22 16854565 1.58688 0 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
22 16854792 1.58702 5 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
22 16854988 1.58713 10 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
22 16855262 1.58730 15 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000
#Subpopulation order/codes: HG=0
#chm spos epos sgpos egpos n snps NA19152.0 NA19152.1 NA18542.0 NA18542.1 NA18740.0 NA18740.1 NA18989.0 NA18989.1
22 16854565 51244237 1.59 74.12 317408 0 0 0 0 0 0 0 0
#rfmix diploid global ancestry .Q format output
#sample HG
NA19152 1.00000
NA18542 1.00000
NA18740 1.00000
NA18989 1.00000
I donot know why, can I get your feed back?
Best
Xue