An internal function
VP.CP.ZP.un(Z_mat, fixed_effects, control)
Z_mat
data frame
fixed_effects
formula specifying fixed effects to be included in model
control
a list