Could you please check whether this works properly? x<-matrix(1,nrow=100,ncol=1) stem(x, atom=1.e-7) The problem appears to be an integer overflow; this fixes it on my system. Duncan Murdoch