explain pi0 estimation in manual for tutorial + add input checks in R functions