> x
V1 V2 V3
1 aa 2001 55
2 aa 2002 66
3 aa 2003 77
4 bb 2002 88
5 bb 2003 99
> y <- unlist(tapply(x$V2, x$V1, order))
> y
aa1 aa2 aa3 bb1 bb2
1 2 3 1 2
> x <- cbind(x, y)
> x
V1 V2 V3 y
aa1 aa 2001 55 1
aa2 aa 2002 66 2
aa3 aa 2003 77 3
bb1 bb 2002 88 1
bb2 bb 2003 99 2