Ах, я пропустил это, поскольку в данных вашего образца нет двух свиноматок в один и тот же день. Вы можете сделать как .numeric (factor (paste (piglets $ fdate, piglets $ sow, sep = ""))), чтобы посмотреть комбинацию этих двух.
ел бы создать новую переменную, помет, чтобы указать каждую свиноматку или помет в разные даты опороса (fdate). Каждый помет должен быть пронумерован от 1 до N с шагом 1, как показано в последнем столбце.
sow season piglet fdate litter
1M521 1 5702 14/09/2009 1
1M521 1 5703 14/09/2009 1
1M521 2 22920 17/02/2010 2
1M521 2 22920 17/02/2010 2
1M521 2 22920 17/02/2010 2
1M584 1 8516 28/09/2009 3
1M584 1 8516 28/09/2009 3
1M584 1 8516 28/09/2009 3
1N312 1 6192 16/09/2009 4
1N312 1 6193 16/09/2009 4
1N312 1 6194 16/09/2009 4
1N312 2 21818 11/02/2010 5
1N312 2 21819 11/02/2010 5
1N312 2 21820 11/02/2010 5
2A1526 1 9293 3/10/2009 6
2A1526 1 9294 3/10/2009 6
2A1526 1 9295 3/10/2009 6
2A1526 1 9296 3/10/2009 6