Plot field vectors in a spherical shell

VectorPlot3D[{Sin[x], Sin[y] ^ 2, Cos[z] ^ 3 - Sin[x]}, {x, -5, 5}, {y, -5, 5}, {z, -5, 5}, VectorPoints -> 12, VectorColorFunction -> "TemperatureMap", RegionFunction -> Function[{x, y, z}, 4 ^ 2 < x ^ 2 + y ^ 2 + z ^ 2 < 5 ^ 2]]

-Graphics3D-