f1 <- as.data.frame(table(aaa))
f1 <- f1[f1$Freq != 0, ]

aaa<-data.frame(x=c(1,1,1,2,2,2,3,3),y=c(1,1,2,3,2,4,5,5))
f <- c(2,1,1,1,1,2)
e <- cbind(aaa[-c(1,8),],f)

   x y f
2 1 1 2
3 1 2 1
4 2 3 1
5 2 2 1
6 2 4 1
7 3 5 2

但是rownames,变为 2:7，为什么不是1:6？？？？？？

rownames  2:7有它的好处

如果要1:6
h <-data.frame(e,row.names=1:length(e[,1]))
   x y f
1 1 1 2
2 1 2 1
3 2 3 1
4 2 2 1
5 2 4 1
6 3 5 2

dxystata :

能写的再清楚点撒？ f是否跟 x,y 关联，x, y是否要重新排序。。。还要删掉 头尾？

楼主要的是把相同x，y的组合的频率求出来