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

数据表中带整数或双精度的矩阵乘法

data.table r

Z117 · 技术社区 · 7 年前

有人能给我解释一下为什么会这样:

a <- data.table(t(c(1, 2, 3, 4, 5)))
b <- matrix(data=1, nrow=10, ncol=5)
a[,lapply(.SD,function(x)(x*b[,x]))]

但事实并非如此:

a <- data.table(t(c(0.1, 0.2, 0.3, 0.4, 0.5)))
b <- matrix(data=1, nrow=10, ncol=5)
a[,lapply(.SD,function(x)(x*b[,x]))]

1 回复 | 直到 7 年前

akrun 7 年前

在第一种情况下,“a”的列值作为列索引来子集“b”,而在第二种情况下则不会,因为值是0.1到0.5。因此,不要循环遍历列,而是循环遍历列的序列,然后子集“b”和列值,这在这两种情况下都适用

a[, lapply(seq_along(.SD),function(x)(.SD[[x]]*b[,x]))]

或者只是复制“a”的行,然后与“b”相乘

a[rep(seq_len(.N), nrow(b))] * b

推荐文章

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 月前