Reordering matrix in bLP

[stouffer]

Is it necessary to reorder the adjacency matrix like in the third line of the function bLP (i.e., x <- x[sample(c(1:NROW(x))),sample(c(1:NCOL(x)))]) or at all?

I already get classic cannot allocate vector of size errors at that stage with large-ish datasets.

[tpoisot]

That's a good question. I think it is because Liu & Murata have a paragraph on the paper about random order for LP to avoid resonance/priority effect, which is most likely why I added this line (but I wrote that in early 2011, so I can't guarantee anything about my motivations).

Let me push some code that would (perhaps) allow to keep this in place.