All.joiner problem: Difference between revisions

From genomewiki
Jump to navigationJump to search
(all.joiner problem to be aware of)
 
(Replacing page with 'This page is no longer maintained.')
 
Line 1: Line 1:
There is a problem with all.joiner rules and the -database parameter.  This is taken from an email from Ann in 2007:
This page is no longer maintained.
 
Bob and I just ran into something sneaky with running joinerCheck or runJoiner.csh.  There are some rules that can be tricky and we are not exactly sure why.  What we do know is the kind of output you can be on the lookout for to know that you're dealing with a "tricky" rule.  The output sort of makes it look like the rule has run (with no problems), but it really hasn't run at all.
 
Example of rules that aren't really running from joinerCheck:
[hgwdev:~/joiner> joinerCheck -keys -database=fr2 -identifier=blastHg18KGPep04Id all.joiner
Checking keys on database fr2
 
from runJoiner.csh:
[hgwdev:~/joiner> runJoiner.csh fr2 blastHg18KG
found identifiers:
blastHg18KGRef04Id
blastHg18KGPep04Id
Checking keys on database fr2
Checking keys on database fr2
 
 
If the rule had really run there would be some output saying something like X hits out of Y ok.
 
The way to force these rules to run is to leave off the "database" parameter (you can only do this when using the joinerCheck (runJoiner.csh can't handle that)).  This will run the rule for all databases that are listed -- just pick yours out of the output.
 
[hgwdev:~/joiner> joinerCheck -keys -identifier=blastHg18KGPep04Id all.joiner Checking keys on database (null)
  anoCar1.blastHg18KG.qName - hits 45332 of 45332 ok
  equCab1.blastHg18KG.qName - hits 61628 of 61628 ok
  ornAna1.blastHg18KG.qName - hits 35761 of 35761 ok
  oryLat1.blastHg18KG.qName - hits 38111 of 38111 ok
  bosTau2.blastHg18KG.qName - hits 56832 of 56832 ok
  strPur1.blastHg18KG.qName - hits 30907 of 30907 ok
  mm8.blastHg18KG.qName - hits 70754 of 70754 ok
  fr2.blastHg18KG.qName - hits 40510 of 40510 ok
  dm2.blastHg18KG.qName - hits 22239 of 22239 ok
  felCat3.blastHg18KG.qName - hits 37715 of 37715 ok
  canFam2.blastHg18KG.qName - hits 62741 of 62741 ok
  danRer4.blastHg18KG.qName - hits 44692 of 44692 ok
  xenTro2.blastHg18KG.qName - hits 40291 of 40291 ok
  galGal3.blastHg18KG.qName - hits 35851 of 35851 ok
  monDom4.blastHg18KG.qName - hits 63264 of 63264 ok
  hg18.blastKGRef04.acc - hits 36727 of 36727 ok
  hg18.knownGene.name - hits 0 of 56722
Error: 56722 of 56722 elements of hg18.knownGene.name are not in key blastKGPep04.name line 1995 of all.joiner
Example miss: UC001AAA.1
 
[[Category:Browser QA tracks]]

Latest revision as of 19:31, 10 March 2011

This page is no longer maintained.