subsets <- function(n, k, set = 1:n)

# POTENZMENGE
# n ÜBER k AUSZÄHLEN
# n Durchläufe, um alle Kombinationen auszuzählen ohne (!) Nullmenge
# aus
# http://www.stat.unibe.ch/~duembgen/EinfStat/R/splus2000.pdf

if(k <= 0) NULL else
if(k >= n) set else
rbind(cbind(set[1], Recall(n - 1, k - 1, set[-1])),
Recall(n - 1, k, set[-1]))