通过Rprof监测，发现closepairs函数是限制瓶颈，因此需要再写个c++函数来加快速度，还未完成