The purpose of this study was to estimate (co) variance components of three milk production traits for genetic evaluation using a multiple lactation model. Each of the first five lactations was treated as different traits. For the parameter estimation study, a data set was set up including lactations from cows calved from 2001 to 2009. The total number of raw lactation records in first to fifth parities reached 1,416,589. At least 10 cows were required for each contemporary group, herd-year-season effect. Sires with fewer than 10 daughters were discarded. Lactations with 305d milk yield exceeding 15,000 kg were removed. In total, 1,456 sires of cows were remained after all the selection steps. A complete pedigree consisting of 292,382 records was used for the study. A sire model containing herd-year-season, caving age, and sire additive genetic effects was applied to the selected lactation data and pedigree for estimating (co) variance components via VCE. Heritabilities and genetic or residual correlations were then derived from the (co) variance estimates using R package. Genetic correlations between lactations ranged from 0.76 to 0.98 for milk yield, 0.79~1.00 for fat yield, 0.75~1.00 for protein yield. On individual lactation basis, relatively low heritability values were obtained 0.14~0.23, 0.13~0.20 and 0.14~0.19 for milk, fat, and protein yields, respectively. For the combined lactation heritability values were 0.29, 0.28, and 0.26 for milk, fat, and protein yields. The estimated parameters will be used in national genetic evaluations for production traits.