merge the overlaping bases with sort -k1,1 -k2,2n *.bed | bedtools merge
This problem was solved by changing na.rm value from na.rm=T to na.rm=TRUE in the R scripts.
Problem in preparing analysis