To use pdiff, you'll need to load it with load("pdiff"). The load statement got cut off my example. Barton