0

我已经用 Rcpp 和 C++ 为 R 编写了一些代码，以尝试更加熟悉它：

#include <Rcpp.h>
#include <vector>

using namespace Rcpp;

// [[Rcpp::export]]
CharacterMatrix reduce_sequences(CharacterMatrix completeDNA)
{
  std::vector<int> informativeSites; 
  for(int i = 0; i < completeDNA.ncol(); i++)
  {
    CharacterVector bpsite(completeDNA.nrow());
    for(int n = 0; n < completeDNA.nrow(); n++)
    {
      bpsite[n] = completeDNA(n,i); 
    }
    if(any(bpsite != bpsite[0]).is_true()) informativeSites.push_back(i);
  }
  CharacterMatrix cutDNA(3, informativeSites.size());
  for(int i = 0; i < informativeSites.size(); i++)
  {
    for(int n = 0; n < cutDNA.nrow(); n++)
    {
      cutDNA(n,i) = completeDNA(n,informativeSites[i]);
    }
  }
  return cutDNA;
}

但我得到一个完整的错误，但不是来自我的源文件，而是来自 Comparator_With_One_Value.h：

我不会假装完全理解这些错误，因为我还处于 C++ 的初级阶段，但是通过适当地注释掉我的代码并找出导致它的原因，这是我的第 17 行：

if(any(bpsite != bpsite[0]).is_true()) informativeSites.push_back(i);

我认为这与我使用任何（）有关。我做错了什么？

编辑：更改行以反映上面解决的所有问题，除了两个：控制台输出：

Error in Rcpp::sourceCpp("reduceseq.cpp") : 
  Error 1 occurred building shared library.

从 Comparator_With_One_Value.h operands to ?: have different types 'SEXPREC*' and 'int' 和 返回的问题invalid conversion from 'SEXPREC* const' to 'int'

谢谢，本。

4

1 回答 1

2

我们故意阻止这种情况，因为R: TRUE, FALSE,中有 3 个逻辑值NA。你应该可以这样使用is_true：

if(any(bpsite != bpsite[0]).is_true()) informativeSites.push_back(i);
于 2013-05-15T08:24:51.473 回答