I'm trying to get a handle on how LD patterns define blocks which are often interpreted as haplotype boundaries. I have LD data (D', r^2) and would like to come up with or use an existing function which takes this matrix and can generate row IDs for the boundaries of each "block".
Is there a BioConductor package or other tool which does this, or a paper which describes a simple approach I could try? The matrix looks "block-diagonal" so I imagine even a technique from linear algebra might have a simple procedure, e.g. this Matlab function (which I don't have access to)?