小编Dat*_*'oh的帖子

我在 R 中使用 qgraph。它会自动缩短节点标签，这非常有用。但我如何要求它使用完整标签呢？

这是一个例子。

require(qgraph)
m <- matrix(1,3,3)
colnames(m) <- c("AAAAAAAA", "BBBBBBBB", "CCCCCCCC")
qgraph(m)

我想从矩阵列表中创建交错矩阵.

输入示例:

> l <- list(a=matrix(1:4,2),b=matrix(5:8,2))
> l
$a
     [,1] [,2]
[1,]    1    3
[2,]    2    4

$b
     [,1] [,2]
[1,]    5    7
[2,]    6    8

预期产量:

1    3
5    7
2    4
6    8

我已经检查了gdata中的interleave函数,但它没有显示列表的这种行为.任何帮助赞赏.

我想从表示为2d数组的对称矩阵的下三角形创建一个字典.例如,如果numpy数组是;

array([[0, 2, 3],
       [2, 0, 4],
       [3, 4, 0]])

然后我希望字典看起来像;

{('1', '0'): 2, ('2', '0'): 3, ('2', '1'): 4}

矢量有类似的帖子;

将Numpy数组转换为稀疏字典的最快方法？

我对python比较新,所以任何帮助/建议都赞赏.

目的是在 R 中实现正交投影非负矩阵分解 (opnmf) 的快速版本。我正在翻译此处提供的 matlab 代码。

我实现了一个普通的 R 版本，但对于 20 因子解决方案，它比我的数据 (~ 225000 x 150) 上的 matlab 实现慢得多（大约慢 5.5 倍）。

所以我认为使用 c++ 可能会加快速度，但它的速度与 R 类似。我认为这可以优化，但不确定如何优化，因为我是 c++ 的新手。这是讨论类似问题的线程。

这是我的 RcppArmadillo 实现。

// [[Rcpp::export]]
Rcpp::List arma_opnmf(const arma::mat & X, const arma::mat & W0, double tol=0.00001, int maxiter=10000, double eps=1e-16) {
  arma::mat W = W0;
  arma::mat Wold = W;
  arma::mat XXW = X * (X.t()*W);
  double diffW = 9999999999.9;
  
  Rcout << "The value of maxiter : " << …