Hi Wacek, Wacek Kusnierczyk schrieb: > ... or gain a bit on performance by doing the threshold comparison on > the whole matrix just once at once: > > dd = d <= 1 > d[apply(d, 1, any), apply(d, 2, any)] d[apply(dd, 1, any), apply(dd, 2, any)] Or not? Cheers, Stephan