1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
"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)
}
|
