朱莉娅的扫描相当于

Question

来自R文档:

sweep:通过扫描摘要统计信息返回从输入数组中获取的数组.

例如,以下是我如何将每一行除以行总和:

> rs = rowSums(attitude)
> ratios = sweep(attitude, 1, rs, FUN="/")
> head(ratios)
     rating complaints privileges  learning    raises  critical    advance
1 0.1191136  0.1412742 0.08310249 0.1080332 0.1689751 0.2548476 0.12465374
2 0.1518072  0.1542169 0.12289157 0.1301205 0.1518072 0.1759036 0.11325301
3 0.1454918  0.1434426 0.13934426 0.1413934 0.1557377 0.1762295 0.09836066
4 0.1568123  0.1619537 0.11568123 0.1208226 0.1388175 0.2159383 0.08997429
5 0.1680498  0.1618257 0.11618257 0.1369295 0.1473029 0.1721992 0.09751037
6 0.1310976  0.1676829 0.14939024 0.1341463 0.1646341 0.1493902 0.10365854
> rowSums(ratios) # check that ratios sum up to 1
 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

我在朱莉娅的尝试:

x = rand(3, 4)
x[1, 1] = 10
x[2, 1] = 20
x[3, 1] = 30
rowsum = sum(x, 2)
rowsum_mat = repmat(rowsum, 1, size(x, 2))
x = x ./ rowsum_mat

这很有效,但很笨重.有没有更优雅和有效的方式这样做？

Answer 1

无需使用repmat- 所有Julia的.操作员都默认进行"广播".这意味着它匹配两个参数的尺寸,然后展开任何长度为1(单例尺寸)的尺寸以匹配另一个数组.由于缩减保持源阵列的相同维度,因此它们可以直接与任何点运算符一起使用.

在您的情况下,您可以使用:

x ./ sum(x, 2)

以来:

julia> x ./ rowsum_mat == x ./ rowsum
true