An R Implementation of the Polya-Aeppli Distribution