In Rcpp, how to get a user-defined structure from C into R

后端 未结 1 1886
萌比男神i
萌比男神i 2020-12-04 02:29

Am using Rcpp packages and can get my C function to compile and run in R, but now I want to return a large, user-defined data structure to R. The fields in the structure are

相关标签:
1条回答
  • 2020-12-04 02:57

    The right structure in R depends on what your struct looks like exactly. A named list is the most general one. Here a simple sample implementation for a wrap function as referred to in the comments:

    #include <RcppCommon.h>
    
    typedef struct {
      char*   firstname[128];
      char*   lastname[128];
      int      nbrOfSamples;
    } HEADER_INFO;
    
    namespace Rcpp {
      template <>
      SEXP wrap(const HEADER_INFO& x);
    }
    
    #include <Rcpp.h>
    
    namespace Rcpp {
      template <>
      SEXP wrap(const HEADER_INFO& x) {
        Rcpp::CharacterVector firstname(x.firstname, x.firstname + x.nbrOfSamples);
        Rcpp::CharacterVector lastname(x.lastname, x.lastname + x.nbrOfSamples);
        return Rcpp::wrap(Rcpp::List::create(Rcpp::Named("firstname") = firstname,
                          Rcpp::Named("lastname") = lastname,
                          Rcpp::Named("nbrOfSamples") = Rcpp::wrap(x.nbrOfSamples)));
      };
    }
    
    //  [[Rcpp::export]]
    HEADER_INFO getHeaderInfo() {
      HEADER_INFO header;
      header.firstname[0] = (char*)"Albert";
      header.lastname[0] = (char*)"Einstein";
      header.firstname[1] = (char*)"Niels";
      header.lastname[1] = (char*)"Bohr";
      header.firstname[2] = (char*)"Werner";
      header.lastname[2] = (char*)"Heisenberg";
      header.nbrOfSamples = 3;
      return header;
    }
    
    /*** R
    getHeaderInfo()
     */
    

    Output:

    > getHeaderInfo()
    $firstname
    [1] "Albert" "Niels"   "Werner"
    
    $lastname
    [1] "Einstein"   "Bohr"       "Heisenberg"
    
    $nbrOfSamples
    [1] 3
    

    However, for this particular case a data.frame would be more natural to use, which can be achieved by replacing above wrap with:

      template <>
      SEXP wrap(const HEADER_INFO& x) {
        Rcpp::CharacterVector firstname(x.firstname, x.firstname + x.nbrOfSamples);
        Rcpp::CharacterVector lastname(x.lastname, x.lastname + x.nbrOfSamples);
        return Rcpp::wrap(Rcpp::DataFrame::create(Rcpp::Named("firstname") = firstname,
                                                  Rcpp::Named("lastname") = lastname));
      };
    

    Output:

    > getHeaderInfo()
      firstname   lastname
    1    Albert   Einstein
    2     Niels       Bohr
    3    Werner Heisenberg
    
    0 讨论(0)
提交回复
热议问题