小编JMF*_*MFA的帖子

我目前正在使用矩阵,我想找到每行中最低的正值.

使用apply(my.matrix,1,min)将不起作用,因为输出将始终为0 ...

有没有办法找到除0以外的最低值？

我试图计算100行x 2500列表中每行之间的成对差异数.

我有一个小的RScript,但运行时间(显然)非常高!我试图在C中编写一个循环,但我不断收到错误(compileCode).

您是否知道如何将以下循环"转换"为C？

pw.dist <- function (vec1, vec2) {

return( length(which(vec1!=vec2)) )

}

N.row <- dim(table)[1]
pw.dist.table <- array( dim = c(dim(table)[1], dim(table)[1]))

for (i in 1:N.row) {
    for (j in 1:N.row) {
        pw.dist.table[i,j] <- pw.dist(table[i,-c(1)], table[j,-c(1)])
    }
}

我正在尝试这样的事情:

sig <- signature(N.row="integer", table="integer", pw.dist.table="integer")
code <- "
  for( int i = 0; i < (*N.row) - 1; i++ ) {
    for( int j = i + 1; j < *N.row; j++ ) {
      int pw.dist.table = table[j] - …