Here is one way to extend it using a simple union bound to control the deviation of every coordinate at once. You could imagine other ways to extend it, but ensuring that every coordinate is near the true value is already very strong.
For $X_1, ..., X_n \in [0,1]$ be iid from some distribution with mean $\mu$, Hoefdding
says : $P ( | \bar{X} - \mu| > \epsilon) \leq 2\exp( - 2 n \epsilon^2) $.
Suppose we have iid vectors $Y_1, Y_2, \ldots, Y_n \in [0,1]^m$. For a vector $Z \in [0,1]^m$ let $Z(j)$ denote the $j$th coordinate. In particular, $(Y_k)(j)$ denotes the $jth$ coordinate of the $k$th vector in the sample.
Suppose $\mathbb{E}[Y_1(j)] = \mu_j$. We have $\mathbb{E}[Y](j) = \mathbb{E}[Y_1(j)]$ by linearity of expectation, and $\overline{Y}(j) = \overline {Y(j)}$ similarly, where the overline denotes the empirical mean.
We define the event that the $j$th coordinate of $\overline{Y}$ is far from the true mean of that coordinate: $A_j = \{ | \bar{Y}(j) - \mu_j| > \epsilon \}$. Hoeffding's tells us that $P(A_j) \leq 2\exp( - 2 n \epsilon^2) $.
We can bound the probability that any coordinate is $\epsilon$ away from the true mean of that coordinate by using the union bound. That is, $P ( \bigcup_{j = 1}^m A_j ) \leq \sum_{j = 1}^m P(A_j) \leq 2m \exp( -2n \epsilon^2)$.
This is useful as we only need to increase $n$ by adding $\log(m)/(2\epsilon^2)$ in order to eliminate the multiplier.