optional data.frame output for pi/theta/tau

Changes for Rjournal resubmission

fixing half-closed notation and examples

adding mean.equals.sd arg and updating docs

Need to use long longs to avoid overrun problem on 32 bit machines.

Updated version.

Deleted blank lines.