I have an equalized constellation for 16 QAM. The constellation is equalized by LMS algorithm.

I want to calculate the EVM for the equalized constellation. How can I calculate this so that the EVM lies within the threshold limit?