Here is a deviant of a question I'veI feel like I have seen several times on truncated exponentials and similar distributions for finding sufficient statistics:
Let $$\mathbb{P}(Y=y)=\theta^y(1-\theta)^{\mathbb{I}\{y\leq z\}}$$ for some $z\in\mathbb{N}$ and $y\in\mathbb{N}$. Is the MLE for $\hat{\theta}$, $\frac{n\bar{Y}-n}{n-\sum\limits_{i=1}^n \mathbb{I}\{y=z+1\}}$, sufficient? If not, what else aside from the entire data $Y=(y_1,\dots,y_n)$, would constitute a sufficient statistic?
My Attempt: I believe that $T(Y)=\left(\sum y_i,\sum \mathbb{I}\{y=z+1\})\right)$ would easily work as a sufficient statistic. The MLE is a function (not 1-1 though) of this statistic which checks the box. As for if the MLE itself is sufficient (which intuitively I don't think is), we have
$$\mathcal{L}(\theta)=\theta^{\sum y_i}(1-\theta)^{n-\sum \mathbb{I}\{y_i=z+1\}},$$
but by inspection, I really see no way of (successfuly) using the factorization theorem here nor anything dealing with the exponential family as this does not reside in that family of distributions.
Question: Is there a more rigorous or correct way of going about this? I am being very hand-wavy here and was wondering if there's an alternative?