swh:1:snp:16c54c84bc54885e783d4424d714e5cc82f479a1
Tip revision: 948ad5967298b4dc174f4d96511af60f38e9279e authored by Roger Koenker on 27 July 2019, 09:38:23 UTC
version 5.42.1
version 5.42.1
Tip revision: 948ad59
combos.R
"combos" <- function(n,p){
if(length(n) != 1){
n <- n[1]
warning("Using first element as n")
}
if(length(p) != 1){
p <- p[1]
warning("Using first element as p")
}
if(n != as.integer(n)){
warning("Coercing n to integer")
n <- as.integer(n)
}
if(p != as.integer(p)){
warning("Coercing p to integer")
p <- as.integer(p)
}
if(p > n) stop("p is greater than n")
m <- choose(n,p)
z <- .Fortran("combin",
as.integer(n),
as.integer(p),
as.integer(m),
a = integer(p*m),
integer(n),
integer(n),
integer(n))
matrix(z$a,p)
}