I have a question regarding coordinate transformations. When we have a vector, say $X^\mu\partial_\mu$, in some coordinates $\xi^\alpha$ and perform a change of coordinates $\xi^\alpha\to\eta^\alpha$, the components of the original vector written in the new coordinates can be retrieved using the vector transformation law:
$$\tilde{X}^\alpha=\dfrac{\partial\eta^\alpha}{\partial\xi^\beta}X^\beta$$
In certain cases we perform an infinitesimal change of coordinates $\xi^\mu \to \xi^\mu + \epsilon v^\mu(\xi)$. Applying the above formula, one would expect to get
$$\tilde{X}^\alpha=\dfrac{\partial(\xi^\alpha+\epsilon v^\alpha(\xi))}{\partial\xi^\beta}X^\beta=\delta^\alpha_\beta X^\beta+\epsilon\dfrac{\partial v^\alpha}{\partial \xi^\beta}X^\beta=X^\alpha+\epsilon\dfrac{\partial v^\alpha}{\partial \xi^\beta}X^\beta$$
However, in How does a vector field transform under an infinitesimal coordinate transformation? it is stated that the transformation requires the use of the Lie derivative. The result is similar if we perform a Taylor expansion:
$$\tilde{X}^\alpha=X^\alpha(\xi^\beta+\epsilon v^\beta(\xi))\simeq X^\alpha(\xi^\beta)+\epsilon v^\nu\dfrac{\partial X^\alpha}{\partial\xi^\nu}=X^\alpha+\epsilon v^\nu\dfrac{\partial X^\alpha}{\partial\xi^\nu}$$
However, $\epsilon\dfrac{\partial v^\alpha}{\partial \xi^\beta}X^\beta\neq \epsilon v^\nu\dfrac{\partial X^\alpha}{\partial\xi^\nu}$.
Why does this happen? Which is the correct approach to this change of coordinates?
If the correct approach is the Lie derivative, why is the formula for the change of coordinates used in Metric components transformation under change of coordinates instead of the Lie derivative?