This paper analyzes the performance of cooperative hybrid automatic repeat request with incremental redundancy (HARQ-IR) and proposes a new approach of outage probability approximation for performance analysis. A general time-correlated Nakagami fading channel covering fast fading and Rayleigh fading as special cases is considered here. An efficient inverse moment matching method is proposed to approximate the outage probability in closed-form. The effect of approximation degree is theoretically analyzed to ease its selection. Moreover, diversity order of cooperative HARQ-IR is analyzed. It is proved that diversity order is irrelevant to the time correlation coefficient $ρ$ as long as $ρ<1$ and full diversity from both spatial and time domains can be achieved by cooperative HARQ-IR under time-correlated fading channels. The accuracy of the analytical results is verified by computer simulations and the results reveal that cooperative HARQ-IR scheme can benefit from high fading order and low channel time correlation. Optimal rate selection to maximize the long term average throughput given a maximum allowable outage probability is finally discussed as one application of the analytical results.