如何解决如何从 Julia 中的大矩阵中随机子集一个矩阵
我有一个大矩阵(28960×45807 Array{Float64,2}),其中行代表个体,列代表snpID。现在我想得到一个子集矩阵(28960x4580),它是从大矩阵中随机选择的。 我怎样才能在 Julia 中做到这一点?
解决方法
假设您的矩阵是 x
:
using StatsBase
@view x[:,sample(1:size(x,2),4580,replace=false,ordered=true)]
说明:
- 使用
@view
可避免数据复制。可以在此命令中跳过它,但会降低性能 - 冒号
:
选择所有行的数据,数组切片的第二个参数用于选择一堆列 -
size(x,2)
返回列数 - 我们使用
sample
中的StatsBase
对列号进行采样而不进行替换。此外,我假设您不想更改列的顺序
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。