>> p.ov <- overlay[onePLOT["Id"],p] > Error in overlay[onePLOT["Id"], p] : > object of type 'closure' is not subsettable > > thanks for help or suggestions 'overlay' is a function, you are using square brackets! Try overlay(a,b) instead of overlay[a,b] Barry