Non-Gaussian Chi-squared method with the multivariate Edgeworth expansion