DNA methylation patterns have profound impacts on genome stability, gene expression, and development. The molecular base of DNA methylation patterns has long been focused at single CpG sites level. Here, we construct a kinetic model of DNA methylation with collaborations between CpG sites, from which a correlation function was established based on experimental data. The function consists of three parts that suggest three possible sources of the correlation: movement of enzymes along DNA, collaboration between DNA methylation and nucleosome modification, and global enzyme concentrations within a cell. Moreover, the collaboration strength between DNA methylation and nucleosome modification is universal for mouse early embryo cells. The obtained correlation function provide insightful understanding for the mechanisms of inheritance of DNA methylation patterns.