ooops... dashed too fast on the "Enter" key. The last line of the example should be: seq.pm <- theProbeSet$RESIDUES[(xy.pm[,1] == theProbeSet$POSITIONX + 1) & (xy.pm[,2] == theProbeSet$POSITIONY + 1)] L.