The aim of this study was to estimate genetic parameters for lactation yields of milk (MY), fat (FY), protein (PY), and somatic cell score (SCS) of New Zealand dairy goats. The analysis used 64,604 lactation records from 23,583 does, kidding between 2004 and 2017, distributed in 21 flocks and representing 915 bucks. Estimates of genetic and residual (co) variances, heritabilities, and repeatabilities were obtained using a multiple-trait repeatability animal model. The model included the fixed effects of contemporary group (does kidding in the same flock and year), age of the doe (in years), and as covariates, kidding day, proportion of Alpine, Nubian, Toggenburg, and "unknown" breeds (Saanen was used as the base breed), and heterosis. Random effects included additive animal genetic and doe permanent environmental effects. Estimates of heritabilities were 0.25 for MY, 0.24 for FY, 0.24 for PY, and 0.21 for SCS. The phenotypic correlations between MY, FY, and PY ranged from 0.90 to 0.96, and the genetic correlations ranged from 0.81 to 0.93. These results indicate lactation yield traits exhibit useful heritable variation and that multiple trait selection for these traits could improve milk revenue produced from successive generations of New Zealand dairy goats.