用其他软件有什么好方法也可以说说呀各位，毕竟都是相通的

dat1 = data.frame(x = c("a","a","a","b","c","c"), y = c(1,2,3,2,2,2))
  
  
  
  
  ff = function(x0){
    return(tabulate(dat1$y[dat1$x == x0], nbins = 3))
  }
  t(sapply(levels(dat1$x), ff))

>   t(sapply(levels(dat1$x), ff))
  [,1] [,2] [,3]
a    1    1    1
b    0    1    0
c    0    2    0

> dat1 = data.frame(x = c("a","a","a","b","c","c"), y = c(1,2,3,2,2,2))
>   
>   ff = function(x0){
+     return(tabulate(dat1$y[dat1$x == x0], nbins = 3))
+   }
>   t(sapply(levels(dat1$x), ff))
  [,1] [,2] [,3]
a    1    1    1
b    0    1    0
c    0    2    0