1### R code from vignette source 'partitionspaper.Rnw'
2
3###################################################
4### code chunk number 1: partitionspaper.Rnw:122-122
5###################################################
6
7
8
9###################################################
10### code chunk number 2: partitionspaper.Rnw:123-124
11###################################################
12require(partitions)
13
14
15###################################################
16### code chunk number 3: partitionspaper.Rnw:132-133
17###################################################
18P(100)
19
20
21###################################################
22### code chunk number 4: partitionspaper.Rnw:137-138
23###################################################
24diffparts(10)
25
26
27###################################################
28### code chunk number 5: partitionspaper.Rnw:170-176
29###################################################
30f <- function(x){prod(factorial(x),factorial(tabulate(x)))}
31prob <- function(a,n){
32  jj <- restrictedparts(n,a,include.zero=FALSE)
33  N <- factorial(a)*factorial(n)*sum(1/apply(jj,2,f))
34  return(N/a^n)
35}
36
37
38