ps <- c(0.2,0.8) N <- 30 library(combinat) calc.survival.2 <- function(X){ ret <- 0 for(i in 1:length(X)){ tmp <- expand.grid(0:(length(X[[i]][,1])-1),0:(length(X[[i]][1,])-1)) tmp.2 <- apply(tmp,1,sum) ret <- ret + sum(tmp.2 * X[[i]]) } r…