r - 如何选择组合之和服从特定分布的样本组合？

Question

我有两个 30 个样本数据集。

A1 <- c(368.4777, 387.9693, 403.7576, 417.4397, 429.7728, 441.1835, 451.9409, 462.2286, 472.1809, 481.9012, 491.4738, 500.9705, 510.4558, 519.9902, 529.6330, 539.4454, 549.4926, 559.8473, 570.5931, 581.8296, 593.6794, 606.2988, 619.8938, 634.7472, 651.2666, 670.0744, 692.1994, 719.5452, 756.3235, 815.7298)

B1 <- c(414.2523,  471.8087,  519.1119,  560.4502,  597.9043,  632.6668,  665.5000,  696.9308,  727.3477,  757.0533,  786.2953,  815.2867,  844.2194,  873.2745,  902.6306,  932.4719,  962.9959,  994.4228, 1027.0068, 1061.0519, 1096.9331, 1135.1297, 1176.2765, 1221.2488, 1271.3134, 1328.4178, 1395.8053, 1479.5340, 1593.1899, 1780.2717)

如果我想对两个分布求和，

1. A1+B1会给我一个绿色虚线密度图。

2. rowSums(expand.grid(A1, B1)) 会给我蓝色实线密度图。 (30*30=900种可能的组合)

但是，我想使总和的分布类似于红线 (Target)。换句话说，我想从 A1 和 B1 中选择组合，使组合之和遵循红线密度分布(每个值只能选择一次)。我该怎么做？

Target <- c(1691.60, 1041.01, 1581.92, 1392.43, 1806.88, 1788.24, 1838.67, 1953.35, 1253.21, 1791.36, 1644.61, 1425.05, 1218.66, 1596.36, 1305.67, 1031.98,  828.52, 1375.77, 1475.28, 1718.25)

Answer 1

我不太理解你问题的最后一部分，所以答案可能需要调整。

如果你想对 A1 和 B1 的组合进行采样，那么你可以只使用 sample 函数。所以

s1 = sample(A1) + sample(B1)
plot(density(s1))

这是没有放回的抽样。如果要替换，只需在 sample 函数中添加 replace=TRUE 即可。