def zinb_ll(data, P, R, Z):
lls = nb_ll(data, P, R)<EOL>clusters = P.shape[<NUM_LIT:1>]<EOL>for c in range(clusters):<EOL><INDENT>pass<EOL><DEDENT>return lls<EOL>
Returns the zero-inflated negative binomial log-likelihood of the data.