The peaks in the cosmic fields have been attracting a lot of interest in the analysis of cosmological fields. In this talk, analytic formulas for the abundance and correlation of peaks in the presence of weak non-Gaussianity are newly derived and presented in N dimensions with N=1,2,3. These formulas provide theoretical methods for the analysis of peak statistics in cosmological fields, such as the galaxy density field or the weak lensing field, etc.