Add genclusters.R