如何找出矩阵的多少行满足相当复杂的标准（在 R 中）？

Question

举个例子，这里有一种方法可以得到一个包含 4 个（公平）骰子的所有可能结果的矩阵。

z <- as.matrix(expand.grid(c(1:6),c(1:6),c(1:6),c(1:6)))

正如您可能已经理解的那样，我正在尝试解决一个已结束的问题，但在我看来，这是一个具有挑战性的问题。我使用计数技术来解决它（我的意思是手工），我最终得出了一些结果，子集之和为 5，等于 1296 中的 1083。该结果与提供给该问题的答案一致，在它关闭之前。我想知道如何使用 R 生成结果的子集（比如 z1，其中 dim(z1) = [1083,4] ）。你有什么想法吗？

谢谢你。

Answer 1

sum(apply(z, 1, function(x) 5 %in% unlist(sapply(1:4, function(i) combn(x, i, sum)))))