I'm using the benford.analysis package to explore the benford's law. I'm being asked to plot distribution of the first digit of any given data set and overlay it on a benford graph. Right now, i'm looking at the plot() of the benford.analysis package, i guess this will fit my requirement, its just that the plot() displays so many "plots"...
bfd.cp <- benford(corporate.payment$Amount,1)
plot(cp, except=c("mantissa", "chi square", "abs diff", "second order"))