tuning_zipath function

find optimal path for penalized zero-inflated model