Question

有关获取子CharacterVector的Rcpp问题。我制作了一个滚动窗口样式函数，该函数可以在多个NumericVectors上使用，但是当包含CharacterVectors时无法使用。每当我尝试使用CharacterVector时，都会收到此错误：“获取临时对象类型代理的地址”。

这是导致问题的代码行：

  CharacterVector subchrom (&chrom[i], &chrom[i+4]);

我是C ++新手，因此目前解决方案超出了我的范围。根据我的阅读，问题在于将指针“＆”用于两种不同类型的向量，但是我不知道如何解决此问题。有人介意帮我解决这个错误吗？

我做了一些阅读，现在我有第二个问题。 Rcpp中的“ CharacterVectors”是否已编程为代理？这是我遇到此问题的唯一向量类型。最好的解决方法是什么？

这是一个完整的例子：

#example data for R   
start <- c(0, 2, 4, 6, 8, 10)
chrom <- c("chr1", "chr1", "chr1", "chr1", "chr1", "chr1")

df <- data.frame(start, chrom)

这是.cpp Rcpp程序：

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

using namespace Rcpp;

// [[Rcpp::export]]
List modifyDataFrame2(DataFrame df) {

  // access the columns
  CharacterVector chrom = df["chrom"];
  IntegerVector start = df["start"];


  List Out = List::create();

  // write the for loop     
  for(int i=0; i < df.nrow(); ++i){

  NumericVector substart (&start[i], &start[i+4]);
  CharacterVector subchrom (&chrom[i], &chrom[i+4]);

  Rcout << "substart is:" <<substart << "\n";
  Rcout << "subchrom is:" <<subchrom << "\n";

  List temp = List::create(subchrom, substart);
  Out.push_back(temp);

  }
  return Out;
}

此刻我在拍摄window函数时遇到了麻烦，因此为什么除了打印子矢量之外它什么都不做。在将计算复杂化之前，我想通过c ++进行输入和输出。

感谢您的帮助。非常感谢您的宝贵时间。

Answer 1

您正在对Rcpp 应该支持的内容做出一些英勇的假设。可悲的是，它并不总是与代码和库一起使用。

这些“灵活”子集不是Rcpp的最强子集，但支持 some 。但是您需要查看实际支持的内容-我只是通过查看现有测试的子集来提醒自己。最古老的方法之一是通过另一个向量建立索引。正如您一直想要的四个元素一样，我只是创建一个索引向量...，然后向其中添加i。

修改后的代码（也在IntegerVector上更正为NumericVector混乱，并更改了索引）：

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
List modifyDataFrame2(DataFrame df) {

  // access the columns
  CharacterVector chrom = df["chrom"];
  IntegerVector start = df["start"];

  List Out;
  IntegerVector ind = IntegerVector::create(0,1,2,3);

  // write the for loop
  for(int i=0; i < df.nrow()-3; ++i){
    IntegerVector currind = ind + i;
    IntegerVector substart = start[currind];
    CharacterVector subchrom = chrom[currind];

    Rcout << "substart is:" << substart << "\n";
    Rcout << "subchrom is:" << subchrom << "\n";
    List temp = List::create(subchrom, substart);
    Out.push_back(temp);
  }
  return Out;
}

/*** R
#example data for R
start <- c(0, 2, 4, 6, 8, 10)
chrom <- c("chr1", "chr1", "chr1", "chr1", "chr1", "chr1")
df <- data.frame(start, chrom)
modifyDataFrame2(df)
*/

当我sourceCpp()时，我们得到以下信息：

R> Rcpp::sourceCpp("~/git/stackoverflow/54190760/answer.cpp")

R> #example data for R
R> start <- c(0, 2, 4, 6, 8, 10)

R> chrom <- c("chr1", "chr1", "chr1", "chr1", "chr1", "chr1")

R> df <- data.frame(start, chrom)

R> modifyDataFrame2(df)
substart is:0 2 4 6
subchrom is:"chr1" "chr1" "chr1" "chr1"
substart is:2 4 6 8
subchrom is:"chr1" "chr1" "chr1" "chr1"
substart is:4 6 8 10
subchrom is:"chr1" "chr1" "chr1" "chr1"
[[1]]
[[1]][[1]]
[1] "chr1" "chr1" "chr1" "chr1"

[[1]][[2]]
[1] 0 2 4 6


[[2]]
[[2]][[1]]
[1] "chr1" "chr1" "chr1" "chr1"

[[2]][[2]]
[1] 2 4 6 8


[[3]]
[[3]][[1]]
[1] "chr1" "chr1" "chr1" "chr1"

[[3]][[2]]
[1]  4  6  8 10


R>

Rcpp：获取临时对象类型代理的地址，从CharacterVector

1 个答案: