To build on user304539's answer: what you want to know is whether starting in the same initial state $\mathbf{x}_0$ again will move you to the same next state $\mathbf{x}_1=\mathbf{f}(\mathbf{x}_0)$ (the system is deterministic) or to some other state $\mathbf{x}_1+\mathbf{y}$ (the system is indeterministic/noisy - $\mathbf{y}$ doesn't depend on $\mathbf{x}_0$). In practice, since you cannot get a perfect re-run, the question is whether the difference between realizations starting close tends to be small or big. Note that we are ideally talking about a short time between state 0 and 1: at late times chaos and randomness are indistinguishable.
One way of checking this is to make a return map of the system: plot the points $(\mathbf{x}_0,\mathbf{x}_1)$. This is obviously simpler for scalar states rather than vector states since you get a 2D plot, but you can often just select some simple projection function $p(\mathbf{x})$ that projects the state onto a single variable. If the system is a deterministic function of state (and later unpredictability is due to chaos) you will get a plot of $\mathbf{f}(\mathbf{x})$. If it is random, you will just get a point cloud. And if it is somewhere in between, you will get a fuzzy shape.
Here is a return plot of iterates of the logistic map $x_{n+1}=\lambda x_n(1-x_n)$ for $\lambda=3.96$ where I plot $(x_n,x_{n+1})$ (leftmost pane). It is a simple parabolic function, although the points $x_n$ are jumping around chaotically. In the right three panes I add normally distributed noise of S.D. 0.01, 0.1 and 1 to the $x_n$ (clamping them to stay in [0,1]). For the strongest noise there is no real pattern, the system is indeterministic. In the intermediate cases there is a fuzzy shape: the system has a mix of deterministic dynamics (which would be chaotic if noise-free) plus random noise.
One thing worth noting in this example is that even adding a tiny bit of noise makes the attractor shape different - the parabola does not touch 0 or 1 in the noise-free case, but does when subjected to noise. Deducing the "true" deterministic behaviour that would happen without noise can be very hard. In this case you can just fit a simple function and when the fit is used for iteration you will get the right dynamics. But noise can make the dynamics drift away and explore parts of state space that would normally be inaccessible - indeed, without the clamping to [0,1] the iterates would run off to infinity in this case, making the return plot useless.