Closed petrelharp closed 9 years ago
Currently, the writeup and code try to make H[x,y] fit D[x,y]; but D is symmetric and H is not; instead, make (H[x,y]+H[y,x])/2 fit D[x,y]. Also, should add a vertical shift to D to account for sequencing error and within-pop divergence.
vertical shift added in 72e41b616c6d9d55fc9a597d78f9610b0f5476ed to trust optimization code
symmetrization added in f994bf8e2a7dcecfbf937b418fc27a06bb815550
Currently, the writeup and code try to make H[x,y] fit D[x,y]; but D is symmetric and H is not; instead, make (H[x,y]+H[y,x])/2 fit D[x,y]. Also, should add a vertical shift to D to account for sequencing error and within-pop divergence.