Changes

From Genome Analysis Wiki
Jump to navigationJump to search
1,118 bytes added ,  10:16, 26 May 2010
Line 511: Line 511:  
| DEFAULT_GENE_FLANK || If only a --refgene is given with no --flank, this variable sets the default flank that is used.  
 
| DEFAULT_GENE_FLANK || If only a --refgene is given with no --flank, this variable sets the default flank that is used.  
 
|}
 
|}
 +
 +
=== LD cache ===
 +
 +
LocusZoom attempts to remember LD calculations that were made on previous runs of the program to avoid having to re-calculate the same regional LD for subsequent runs. The process works as follows:
 +
 +
* For a given reference SNP and chr/start/stop:
 +
** If LD has not been previously computed, use new_fugue to compute LD with the reference SNP and all other SNPs in the region, and store this result to the LD cache
 +
** Else, retrieve the previously stored LD results
 +
 +
The cache intelligently stores LD from separate sources (hapmap, 1000G), populations, builds, and even different versions of genotype files separately.
 +
 +
Upon running LocusZoom, a file called "ld_cache.db" will automatically be created in the current directory, and LD computations will be stored there. If you wish to change the location of the LD cache, use <code>--cache <file></code>. If you wish to disable LD caching, you can use <code>--cache None</code>.
 +
 +
''Python note:'' The ld_cache.db is actually a shelve, and you can explore its contents using:
 +
 +
<code><pre>
 +
import shelve
 +
d = shelve.open("ld_cache.db")
 +
</pre></code>
    
== License  ==
 
== License  ==
    
[[Category:Software]]
 
[[Category:Software]]
239

edits

Navigation menu