I will discuss the holomorphic anomaly equation for GW invariants of several Calabi-Yau geometries. Holomorphic anomaly equation was originated from Physics B-model which recursively determine higher genus GW invariants in terms of lower genus invariants. I will explain how one can study holomorphic anomaly equation in mathematics for several examples. This talk is based on joint work with Rahul Pandharipande.