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)