library(combinat) n <- 10 k1 <- k2 <- 2 k <- k1 + k2 knall <- nsimplex(k,0:n) kn1 <- nsimplex(k1,0:n) kn2 <- nsimplex(k2,n:0) kn1 * kn2 knall X <- list() n <- 4 for(i in 0:n){ X[[i+1]] <- as.matrix(xsimplex(k1,i)) } for(i in 1:length(X)){ …