Problem
As part of a proof for the commutation relation of a vector operator and the angular momentum operator, I need the evaluate the expression
$$(R_\omega)_{ij}V_j = \big([1-\cos(\omega)]\hat{\omega}_i\hat{\omega}_j + \cos(\omega)\delta_{ij} + \sin(\omega)\varepsilon_{ijk}\hat{\omega}_k\big) V_j$$
Here, $(R_\omega)_{ij}$ is the representation for the rotation generator I was given and advised to use, and $V_j$ is the $j$-component of said vector operator. (Eventually, I'm trying to show that $[L_i,V_j]=i\hbar\varepsilon_{ijk}V_k$.)
Attempt
Now, on the LHS, $j$-index is repeated, so I should sum over $j=i,j,k$. I'm still clunky with this convention, so I computed the sum term-wise, (with $\hat{\omega}_i\hat{\omega}_j=\delta_{ij}$ assuming orthonormal unit vectors)
$$ (R_\omega)_{ij}V_j = \big(1-\sin(\omega)\hat{\omega}_k\big)V_j $$
Is this correct? I have reservations, because (a) why was I given $\hat{\omega}_i\hat{\omega}_j$ instead of another $\delta$-function? and (b) this result leaves me with a scalar minus a vector, which seems erroneous to me.
I feel like my execution of the Einstein Summation notation is flawed. I'm familiar with the mechanics of it, but again, I'm clunky with implementation. For instance, am I suppose to sum over $j$ on the LHS and then $i,j$ on the RHS?