- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试创建 dist
来自大距离矩阵的对象。我使用 stats::as.dist
内存不足.例如,我在当前机器上有大约 128 Gb 可用,但 as.dist
处理 73,000 x 73,000 矩阵(约 42Gb)时内存不足。鉴于最后dist
对象应该小于矩阵大小的一半(即它是下三角形,存储为向量)在我看来应该可以以更节省内存的方式进行此计算 - 只要我们小心创建大型中间对象,只需将输入的相关元素直接复制到输出。
查看 getS3method('as.dist', 'default')
的代码,我看到它使用 ans <- m[row(m) > col(m)]
进行计算这需要创建 row
和 col
与输入具有相同维度的矩阵。
我想我可以使用来自 here 的算法对此进行改进。生成下三角形的索引。这是我使用这种方法的尝试。
as.dist.new = function(dm, diag = FALSE, upper = FALSE) {
n = dim(dm)[1]
stopifnot(is.matrix(dm))
stopifnot(dim(dm)[2] == n)
k = 1:((n^2 - n)/2)
j <- floor(((2 * n + 1) - sqrt((2 * n - 1) ^ 2 - 8 * (k - 1))) / 2)
i <- j + k - (2 * n - j) * (j - 1) / 2
idx = cbind(i,j)
remove(i,j,k)
gc()
d = dm[idx]
class(d) <- "dist"
attr(d, "Size") <- n
attr(d, "call") <- match.call()
attr(d, "Diag") <- diag
attr(d, "Upper") <- upper
d
}
N = 10
dm = matrix(runif(N*N), N, N)
diag(dm) = 0
x = as.dist(dm)
y = as.dist.new(dm)
as.dist
相同的内存问题。 .
N = 73000
dm = matrix(runif(N*N), N, N)
gc()
diag(dm) = 0
gc()
as.dist.new(dm)
RcppArmadillo
来做到这一点。 ,但我没有使用它的经验。
最佳答案
嗯,遍历下三角的逻辑比较简单,
如果你用 C++ 来做,那么它也可以很快:
library(Rcpp)
sourceCpp(code='
// [[Rcpp::plugins(cpp11)]]
#include <cstddef> // size_t
#include <Rcpp.h>
using namespace Rcpp;
// [[Rcpp::export]]
NumericVector as_dist(const NumericMatrix& mat) {
std::size_t nrow = mat.nrow();
std::size_t ncol = mat.ncol();
std::size_t size = nrow * (nrow - 1) / 2;
NumericVector ans(size);
if (nrow > 1) {
std::size_t k = 0;
for (std::size_t j = 0; j < ncol; j++) {
for (std::size_t i = j + 1; i < nrow; i++) {
ans[k++] = mat(i,j);
}
}
}
ans.attr("class") = "dist";
ans.attr("Size") = nrow;
ans.attr("Diag") = false;
ans.attr("Upper") = false;
return ans;
}
')
as_dist(matrix(1:100, 10, 10))
1 2 3 4 5 6 7 8 9
2 2
3 3 13
4 4 14 24
5 5 15 25 35
6 6 16 26 36 46
7 7 17 27 37 47 57
8 8 18 28 38 48 58 68
9 9 19 29 39 49 59 69 79
10 10 20 30 40 50 60 70 80 90
关于r - 从距离矩阵创建 dist 对象的内存高效方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/59907035/
前言: 有时候,一个数据库有多个帐号,包括数据库管理员,开发人员,运维支撑人员等,可能有很多帐号都有比较大的权限,例如DDL操作权限(创建,修改,删除存储过程,创建,修改,删除表等),账户多了,管理
所以我用 Create React App 创建并设置了一个大型 React 应用程序。最近我们开始使用 Storybook 来处理和创建组件。它很棒。但是,当我们尝试运行或构建应用程序时,我们不断遇
遵循我正在创建的控件的代码片段。这个控件用在不同的地方,变量也不同。 我正在尝试编写指令来清理代码,但在 {{}} 附近插入值时出现解析错误。 刚接触 Angular ,无法确定我错过了什么。请帮忙。
我正在尝试创建一个 image/jpeg jax-rs 提供程序类,它为我的基于 post rest 的 Web 服务创建一个图像。我无法制定请求来测试以下内容,最简单的测试方法是什么? @POST
我一直在 Windows 10 的模拟器中练习 c。后来我改用dev C++ IDE。当我在 C 中使用 FILE 时。创建的文件的名称为 test.txt ,而我给出了其他名称。请帮助解决它。 下面
当我们创建自定义 View 时,我们将 View 文件的所有者设置为自定义类,并使用 initWithFrame 或 initWithCode 对其进行实例化。 当我们创建 customUITable
我正在尝试为函数 * Producer 创建一个线程,但用于创建线程的行显示错误。我为这句话加了星标,但我无法弄清楚它出了什么问题...... #include #include #include
今天在做项目时,遇到了需要创建JavaScript对象的情况。所以Bing了一篇老外写的关于3种创建JavaScript对象的文章,看后跟着打了一遍代码。感觉方法挺好的,在这里与大家分享一下。 &
我正在阅读将查询字符串传递给 Amazon 的 S3 以进行身份验证的文档,但似乎无法理解 StringToSign 的创建和使用方式。我正在寻找一个具体示例来说明 (1) 如何构造 String
前言:我对 C# 中任务的底层实现不太了解,只了解它们的用法。为我在下面屠宰的任何东西道歉: 对于“我怎样才能开始一项任务但不等待它?”这个问题,我找不到一个好的答案。在 C# 中。更具体地说,即使任
我有一个由一些复杂的表达式生成的 ILookup。假设这是按姓氏查找人。 (在我们简单的世界模型中,姓氏在家庭中是唯一的) ILookup families; 现在我有两个对如何构建感兴趣的查询。 首
我试图创建一个 MSI,其中包含 和 exe。在 WIX 中使用了捆绑选项。这样做时出错。有人可以帮我解决这个问题。下面是代码: 错误 error LGH
在 Yii 中,Create 和 Update 通常使用相同的形式。因此,如果我在创建期间有电子邮件、密码、...other_fields...等字段,但我不想在更新期间专门显示电子邮件和密码字段,但
上周我一直在努力创建一个给定一行和一列的 QModelIndex。 或者,我会满足于在已经存在的 QModelIndex 中更改 row() 的值。 任何帮助,将不胜感激。 编辑: QModelInd
出于某种原因,这不起作用: const char * str_reset_command = "\r\nReset"; const char * str_config_command = "\r\nC
现在,我有以下由 original.df %.% group_by(Category) %.% tally() %.% arrange(desc(n)) 创建的 data.frame。 DF 5),
在今天之前,我使用/etc/vim/vimrc来配置我的vim设置。今天,我想到了创建.vimrc文件。所以,我用 touch .vimrc cat /etc/vim/vimrc > .vimrc 所
我可以创建一个 MKAnnotation,还是只读的?我有坐标,但我发现使用 setCooperative 手动创建 MKAnnotation 并不容易。 想法? 最佳答案 MKAnnotation
在以下代码中,第一个日志语句按预期显示小数,但第二个日志语句记录 NULL。我做错了什么? NSDictionary *entry = [[NSDictionary alloc] initWithOb
我正在使用与此类似的代码动态添加到数组; $arrayF[$f+1][$y][$x+1] = $value+1; 但是我在错误报告中收到了这个: undefined offset :1 问题:尝试创
我是一名优秀的程序员,十分优秀!