Question

我正在尝试确定哪些列是从每行中随机采样的矩阵中采样的。函数示例似乎无法告诉您实际采样的位置。现在，如果所有值都是唯一的，那么简单的匹配例程就可以解决问题。但是，它们不在我的情况下，所以这不起作用。

x <- c(2,3,5,1,6,7,2,3,5,6,3,5)
y <- matrix(x,ncol=4,nrow=3)
random <- t(apply(y,1,sample,2,replace=FALSE))

ý

   [,1] [,2] [,3] [,4]
[1,]    2    1    2    6
[2,]    3    6    3    3
[3,]    5    7    5    5

随机

     [,1] [,2]
[1,]    2    6
[2,]    3    3
[3,]    5    5

对于原始矩阵中的重复值，我无法判断是否从第1列或第3列中采样了随机[1,1]，因为它们都具有值2.因此，匹配在此处不起作用。

随着矩阵“随机”，我还想要一个矩阵，它以相同大小的矩阵给出每个值被采样的列。例如，如：

     [,1] [,2]
[1,]    1    4
[2,]    1    3
[3,]    3    4

谢谢！

Answer 1

您需要单独保存sample的随机选择，以便日后不必担心匹配。例如，再次使用y：

y
#     [,1] [,2] [,3] [,4]
#[1,]    2    1    2    6
#[2,]    3    6    3    3
#[3,]    5    7    5    5

set.seed(42)
randkey <- t(replicate(nrow(y),sample(1:ncol(y),2)))
#     [,1] [,2]
#[1,]    4    3
#[2,]    2    3
#[3,]    3    2

random <- matrix(y[cbind(c(row(randkey)), c(randkey))], nrow(y))
#     [,1] [,2]
#[1,]    6    2
#[2,]    6    3
#[3,]    5    7

R哪些样品是从基质中随机抽取的？

1 个答案: