ca4516b7f5442f06d7eca5bffd7a82d18ef18a3b
kate
  Wed Sep 16 11:18:35 2020 -0700
Interim check-in for exploratory work with COVID GWAS V2 track (moving to branch). refs #26129

diff --git src/hg/makeDb/doc/covid/covidHgiGwas.txt src/hg/makeDb/doc/covid/covidHgiGwas.txt
index c9ee1c0..a9da78a 100644
--- src/hg/makeDb/doc/covid/covidHgiGwas.txt
+++ src/hg/makeDb/doc/covid/covidHgiGwas.txt
@@ -120,15 +120,54 @@
     echo $b.$a.$db
     if ($a == "B2") then
         set studies = 8
     else
         set studies = 18
     endif
     perl $bin/makeCovidHgiGwas.pl $db $studies $f.txt > $f.bed
     bedSort $f.bed $f.sorted.bed
     bedToBigBed -type=bed9+12 -as=$bin/covidHgiGwas.as -tab $f.sorted.bed $sizes $bb.bb
     ln -s `pwd`/$bb.bb /gbdb/hg19/covidHgiGwas
 end
 'EOF'
 
 make makeHgi.csh >&! makeHgi.out &
 
+##################
+# Get effect size distribution to inform color gradient selection
+
+% textHistogram -real -minVal=-25  -binSize=1 -maxBinCount=50  effectSizes.txt 
+-13.000000  1
+-12.000000  1
+-11.000000  0
+-10.000000  2
+-9.000000  3
+-8.000000  6
+-7.000000  6
+-6.000000 * 19
+-5.000000 *** 53
+-4.000000 ********************************************************* 944
+-3.000000 ************************************************************ 1000
+-2.000000 ************************************************************ 1000
+-1.000000 ************************************************************ 1000
+0.000000 ************************************************************ 1001
+1.000000 ************************************************************ 1000
+2.000000 ************************************************************ 1000
+3.000000 ************************************************************ 1000
+4.000000 ************************************************************ 994
+5.000000 ********************************************************** 971
+6.000000 **************************************************** 866
+7.000000 ******************************************* 717
+8.000000 ****************************** 504
+9.000000 ********************** 371
+10.000000 *********** 183
+11.000000 ***** 78
+12.000000 *** 43
+13.000000 ** 37
+14.000000 * 13
+15.000000  4
+16.000000  6
+17.000000  2
+18.000000  0
+19.000000  0
+20.000000  1
+