R自组织子集

Question

我有一个有100行的矩阵。 有没有办法获得十行最接近第一行的子集。

res2 <- matrix(rexp(200, rate=.1), ncol=10, nrow=100) 

set1 <- subset(res2, res2 >condition1) 
set1[with(set1, order(condition)), ] 
set2 <- head(set1,10) 

Answer 1

也许：

生成数据：

set.seed(101) 
res2 <- matrix(rexp(200, rate=.1), ncol=10, nrow=100) 

计算距离矩阵。这是非常低效的，因为我们计算所有成对距离，但它是有效的编码和易于使用，你有很多选择的距离度量（见?dist，寻找method）。对于这个尺寸问题，它非常快。

dd <- dist(res2) 
rr <- rank(as.matrix(dd)[1,]) 

你会发现，在第一行的第一个元素（其是行1和它本身之间的距离）的秩为1，其值（as.matrix(dd)[1,1]）为零。因此，我们现在需要的是具有接下来十个最小距离的行...

res2[rr>1 & rr<=11,]