From Genome Analysis Wiki
Jump to navigationJump to search
6 bytes added
, 07:01, 18 May 2010
Line 99: |
Line 99: |
| mach1 -d sample.dat -p subset.ped -s chr20.snps -h chr20.hap --compact --greedy --autoFlip -r 100 -o par_infer > mach.infer.log | | mach1 -d sample.dat -p subset.ped -s chr20.snps -h chr20.hap --compact --greedy --autoFlip -r 100 -o par_infer > mach.infer.log |
| <br> | | <br> |
− | mach1 -d sample.dat -p sample.ped -s chr20.snps -h chr20.hap --compact --greedy --autoFlip --errorMap par_infer.erate --crossoverMap par_infer.rec --mle --mldetails > mach.imp.log | + | mach1 -d sample.dat -p sample.ped -s chr20.snps -h chr20.hap --compact --greedy --autoFlip <br> |
| + | --errorMap par_infer.erate --crossoverMap par_infer.rec --mle --mldetails > mach.imp.log |
| | | |
| '''Q: Where can I find combined HapMap reference files? '''<br> A: http://www.sph.umich.edu/csg/yli/mach/download/HapMap-r21.html <br><br> | | '''Q: Where can I find combined HapMap reference files? '''<br> A: http://www.sph.umich.edu/csg/yli/mach/download/HapMap-r21.html <br><br> |
Line 144: |
Line 145: |
| The combined files should look like:<br> '''comb.ped'''<br> r1 r1 0 0 1 A/C C/C G/G G/A A/A<br> 1 1 0 0 1 A/A ./. ./. G/G ./. | | The combined files should look like:<br> '''comb.ped'''<br> r1 r1 0 0 1 A/C C/C G/G G/A A/A<br> 1 1 0 0 1 A/A ./. ./. G/G ./. |
| | | |
− | '''comb.dat'''<br> M SNP1<br> M SNP2<br> M SNP3<br> M SNP4<br> M SNP5<br> | + | '''comb.dat'''<br> M SNP1<br> M SNP2<br> M SNP3<br> M SNP4<br> M SNP5<br> |
| | | |
| == Examples == | | == Examples == |