¿Cómo se podría demostrar que la función de masa de probabilidad de X + Y, donde X e Y son variables aleatorias independientes cada una distribuida geométricamente con parámetro p; es decir,
$p_X(n)=p_Y(n)=\left\{\begin{matrix} p(1-p)^{n-1} & n=1,2,...\\ 0 & de lo contrario \end{matrix}\right.$
es igual a $\mathbf{P_{X+Y}(n)= \color{Red}{(n-1)}\ p^2(1-p)^{n-2}}$
Utilizando la convolución, obtengo
$\mathbf{P(X+Y=n)=\sum_{n}^{k=0} Pr(X=k)*Pr(Y=n-k) =\sum_{k=1}^n p_X (1-p_x)^{k-1} p_Y(1-p_Y)^{n-k-1}$
ya que $p=p_X=p_Y$ se reduce a
$\mathbf{P(X+Y=n)=\sum_{k=1}^n p^2(1-p)^{n-2}$
¿Es este un método correcto? Estoy atascado aquí, no sé cómo llegar a la fórmula final. Me falta alguna transición para obtener el (n-1).