예제 #1
0
 def compute_covariance(self):
     breakpoints = BedTool(paths.reference + self.params.breakpointsfile)
     blocks = SnpPartition(self.refpanel, breakpoints, remove_mhc=True)
     myranges = self.ranges_in_chunk(blocks.ranges())
     print('working on', len(myranges), 'ld blocks')
     return BlockDiag.ld_matrix_blocks(self.refpanel, myranges)
예제 #2
0
 def compute_covariance(self):
     breakpoints = BedTool(paths.reference + self.params.breakpointsfile)
     blocks = SnpPartition(self.refpanel, breakpoints, remove_mhc=True)
     myranges = self.ranges_in_chunk(blocks.ranges())
     print('working on', len(myranges), 'ld blocks')
     return BlockDiag.ld_matrix_blocks(self.refpanel, myranges)