The variants included
sex and neuter status, weight, age, and vaccine group.
The model was assessed for significance by use of the le
Cessie–van Houwelingen test. A value of P < 0.05 was considered
significant. Statistical analysis was performed with
R version 2.11.1