[R] generate combination multiset (set with repetition)
Reuben Cummings
reubano at gmail.com
Thu Dec 11 22:59:31 CET 2008
Thanks for the suggestions!! I upgraded R which solved a few of my
problems and got the following functions to work:
# works for m > 2
# requires {prob]
multisets1 <- function(n, m) {
as.matrix(urnsamples(1:n, size = m, ordered = FALSE, replace = TRUE))
}
# works for m > 2
# requires {multic]
multisets2 <- function(n, m) {
subsets(n, m, allow.repeat=T)
}
# doesn't work for m > 2
multisets3 <- function(n, m) {
brks <- cumsum(n:1)
k <- 1:choose(n+1, m)
j <- findInterval(k, brks+1) + 1
i <- k - (brks-brks[1])[j]
cbind(j, i)
}
# doesn't work for m > 2
multisets4 <- function(n, m) {
which(lower.tri(diag(n), diag=TRUE), arr.ind=TRUE )[,m:1]
}
my questions now are... how would I generalize functions 3 and 4 for m
> 2? And, which of the 4 functions would be best for varying ranges of
n and m? I am expecting values for n to range between 1 and 1e3, while
m will range between 1 and 1e6.
Reuben
On Thu, Dec 11, 2008 at 1:26 PM, Charles C. Berry <cberry at tajo.ucsd.edu> wrote:
> On Thu, 11 Dec 2008, Reuben Cummings wrote:
>
>> Hi,
>>
>> This has been asked before but not sufficiently answered from what I
>> could find. How do you create combinations
>> with repetitions (multisets) in R?
>>
>> If I have
>>>
>>> set <- array(1:3)
>
> Why wrap 1:3 in array() ??
>
>>
>> And I want to choose all combinations of picking 2 numbers, I want to
>> get a print out like
>>
>> [,1] [,2]
>> [1,] 1 1
>> [2,] 1 2
>> [3,] 1 3
>> [4,] 2 2
>> [5,] 2 3
>> [6,] 3 3
>>
>
> For small problems (n < 100, say) :
>
> which( lower.tri( diag( n ), diag=TRUE), arr.ind=TRUE )[,2:1]
>
> For larger problems, something like :
>
> foo <- function(n) {
> brks <- cumsum( n:1 )
> k <- 1:choose( n+1, 2 )
> j <- findInterval( k, brks+1 ) + 1
> i <- k - ( brks-brks[1] )[ j ]
> cbind( j, i ) }
>
> If the number in 'set' are not 1:n, you can do a lookup using the results
> from above.
>
>
> HTH,
>
> Chuck
>
>
>
>> subsets(set, 2, allow.repeat=T) should work, but I can't get the
>> multic package to install, t(combn(set,2)) was suggested but it
>> doesn't produce repetitions; expand.grid(rep(list(1:3), 2)) was also
>> suggested, but it produces permuations, not combinations.
>> Additionally, I would like to iterate through each resultant set for
>> large n (similar to the description for getNextSet {pcalg}). Any
>> suggestions?
>>
>> Reuben Cummings
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
> Charles C. Berry (858) 534-2098
> Dept of Family/Preventive
> Medicine
> E mailto:cberry at tajo.ucsd.edu UC San Diego
> http://famprevmed.ucsd.edu/faculty/cberry/ La Jolla, San Diego 92093-0901
>
>
>
More information about the R-help
mailing list