代码之家 › 专栏 › 技术社区 › user5054

有没有办法从矩阵中矢量化列的选择(重复)?

vectorization matrix r

user5054 · 技术社区 · 7 年前

我有一个大小为L的矩阵 n x k 和一个大小为Z的向量 p 。Z由表示L的列索引的整数组成。我想创建一个大小为的矩阵X n x p 它是基于Z中的值选择的L的相应列的聚合。

Z = c(1, 3, 1, 2)
L = matrix(c(73,50,4,14,87,5,34,51,17,57,47,65),nrow=4)
> L
     [,1] [,2] [,3]
[1,]   73   87   17
[2,]   50    5   57
[3,]    4   34   47
[4,]   14   51   65

我希望X是

> X
     [,1] [,2] [,3] [,4]
[1,]   73   17   73   87
[2,]   50   57   50    5
[3,]    4   47    4   34
[4,]   14   65   14   51

在我的原始数据中, p ,则, k 和 n 非常大(分别为30K、500和2K),在所有Z值上循环以选择和组合L中的列需要很长时间。是否有一种矢量化的方法(无循环)来执行此任务?

1 回复 | 直到 7 年前

thelatemail 7 年前

很确定这只是:

L[,Z]
#     [,1] [,2] [,3] [,4]
#[1,]   73   17   73   87
#[2,]   50   57   50    5
#[3,]    4   47    4   34
#[4,]   14   65   14   51

当您从大多数对象中进行选择时,R并不关心是否有重复的列索引。

推荐文章

Amp · 使用R ggplot2删除geom_radial中axis.line和panel.border之间的空格

3 月前

Hard_Course · 用另一列中的值替换行的最后一个非NA条目

3 月前

Mark R · 使用geom_sf()删除地球仪上不需要的网格线

3 月前

Joe · 根据对工作日和本周早些时候的日期的了解,找到一个日期

3 月前

Ben · 统计向量中的单词在字符串中出现的频率

3 月前

TheCodeNovice · R中符号格式的尾随零和其他问题[重复]

3 月前

katefull06 · 在R中使用terra修改范围时,会为单独的SpatRaster重写范围

3 月前

dez93_2000 · 在R管道子功能中引用管道对象的当前状态

4 月前

accibio · 在ggplot2中为同一变量创建两个连续的颜色渐变比例

4 月前

Mankka · 如何在Ggplot2中绘制均匀的径向图

4 月前