将概率应用于多变量矩阵的样本函数

Question

目前，我正在尝试引导一个包含 114 个 obs 和 16 个变量的数据集。

我使用的示例函数如下：

x[sample(nrow(x),size=114,replace=TRUE),] where x is my dataset. 

但是，我想使用分配给特定列的概率进行采样，因为示例函数包含这种可能性。例如，我想以 1-5 0.1 和 5-200 0.9 之间的数字概率对第 5 列进行采样。

我该怎么做？

Answer 1

如果我了解您在寻找什么，这可能是 mapply 的工作。

# fake data
x <- as.data.frame(matrix(1:10,nrow=2))
x
  V1 V2
1  1  6
2  2  7
3  3  8
4  4  9
5  5 10

# fake probabilities of each row, for each column
probs <- as.data.frame(matrix(c(.1,.1,.1,.2,.5,.5,.2,.1,.1,.1),ncol=2))
probs
   V1  V2
1 0.1 0.5
2 0.1 0.2
3 0.1 0.1
4 0.2 0.1
5 0.5 0.1

# then a little mapply magic - change size as needed
mapply(sample, x=x, prob=probs, replace=T, size=10)
      V1 V2
 [1,]  5  7
 [2,]  5  7
 [3,]  2  6
 [4,]  5  6
 [5,]  1  6
 [6,]  5  9
 [7,]  1  9
 [8,]  5  9
 [9,]  5  6
[10,]  4  7