Times (in days) for non-infected cows: timenoninfected <- c(60, 74, 37, 45, 75, 40, 50, 50, 146, 70, 50, 84, 60, 149, 50, 90, 259, 40, 90, 101, 70, 90, 254, 130, 80, 40) Time (in days) for infected cows: timeinfected <- c(50, 130, 100, 130, 50, 140, 129, 76, 138, 69, 70, 144, 70, 130, 70, 150, 251, 110, 120)