The stochastic part of t-SNE is presumably from the randomness of the initial placement of the points in the low-dimensional space.
Does this mean that t-SNE should be re-run on re-seeded randomness to ascertain whether the final embedding of high-dimensional data onto the low-dimensional space is or is not stable? That is at least a rudimentary sensitivity analysis of comparing the results of different runs should be performed?
The minimization of the Kullback–Leibler divergence with respect to the output points is performed using gradient descent. Is stochastic gradient descent usually used here, and if so does this introduce randomness of the output? Would this be of a lesser amount to that introduced by 1., or would the two not be comparable?