Ylaine/2009-8-17

From ZhangLabWiki
Revision as of 18:47, 17 August 2009 by >Ylaine (New page: ==New 1K Genomes Reference== * Get 5 more matches and 5 fewer misses than when trio file used for comparison ** Should be more. Too many mismatches? Debugging code. ** Found bug: using re...)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

New 1K Genomes Reference

  • Get 5 more matches and 5 fewer misses than when trio file used for comparison
    • Should be more. Too many mismatches? Debugging code.
    • Found bug: using reference instead of 1K call for comparison
../Scripts/hetOnly.pl no_dbSNP.new1K >no_dbSNP.new1K.het
  • Quality and coverage alone are the best combination of parameters
  • For n=50, p=0.3, f = 0.75:
' Mean Std
FP 25.26 7.32
MD 15.20 6.17
Overall 18.32 4.02

Thoughts about classification

  • We need a universal classification scheme that doesn't rely on knowing the classes of the training set.
  • A more analogous problem is one of clustering