p_zipath function

compute p-values from penalized zero-inflated model with multi-split data