count.data.in.bins<-function(x,bins){
        n=length(x)
        m=length(bins)
        y=rep(0,m+1)
        for(i in 1:n){
                if(x[i]<=bins[1])y[1]=y[1]+1
                else if(x[i]>bins[m])y[m+1]=y[m+1]+1
        else for(j in 1:(m-1)){if(x[i]>bins[j]&x[i]<=bins[j+1])y[j+1]=y[j+1]+1}
        }
return(y)
}