num_permutations_hack

DeclareDesign · Sep 12, 2017 · 0510b0f · 0510b0f
1 parent 3067917
commit 0510b0f
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 0 deletions.
diff --git a/R/obtain_num_permutations.R b/R/obtain_num_permutations.R
@@ -159,6 +159,9 @@ multinomial_coefficient <-
     # https://www.statlect.com/mathematical-tools/partitions
     if (N != sum(m_each)) {
       stop("the sum of m_each must equal N.")
+    } 
+    if (lfactorial(N) > 300){
+      return(Inf) # this is a hack!
     }
     factorial(N) / prod(factorial(m_each))
   }

diff --git a/tests/testthat/test-permutations.R b/tests/testthat/test-permutations.R
@@ -171,3 +171,9 @@ rowMeans(perms)
 perms %*% perm_probs
 
 
+declaration <- declare_ra(N = 958, m = 479)
+
+obtain_num_permutations(declaration)
+obtain_permutation_matrix(declaration)
+
+