R-programs for the better than n-1 model