Sum of ((observed - expected)^2)/expected)
dividing by df (number of degrees of freedom