add vcfreader@samples()

NEWS.md

@@ -1,3 +1,7 @@
+# vcfppR 0.3.7
+
+* add vcfreader@samples()
+
 # vcfppR 0.3.6
 
 * add `vcfreader::line`

R/RcppExports.R

@@ -61,6 +61,7 @@ heterozygosity <- function(vcffile, region = "", samples = "-", pass = FALSE, qu
 #' @field hasOTHER Test if current variant has a OTHER or not
 #' @field hasOVERLAP Test if current variant has a OVERLAP or not
 #' @field nsamples Return the number of samples
+#' @field samples Return a vector of samples id
 #' @field header Return the raw string of the vcf header
 #' @field string Return the raw string of current variant including newline
 #' @field line Return the raw string of current variant without newline

src/vcf-reader.cpp

@@ -51,6 +51,7 @@ using namespace std;
 //' @field hasOTHER Test if current variant has a OTHER or not
 //' @field hasOVERLAP Test if current variant has a OVERLAP or not
 //' @field nsamples Return the number of samples
+//' @field samples Return a vector of samples id
 //' @field header Return the raw string of the vcf header
 //' @field string Return the raw string of current variant including newline
 //' @field line Return the raw string of current variant without newline
@@ -245,6 +246,7 @@ class vcfreader {
         return v.size() / nsamples();
     }
     inline std::string header() const { return br.header.asString(); }
+    inline std::vector<std::string> samples() const { return br.header.getSamples(); }
     inline std::string string() const { return var.asString(); }
     inline std::string line() {
         std::string s = var.asString();
@@ -363,6 +365,7 @@ RCPP_MODULE(vcfreader) {
         .method("hasOTHER", &vcfreader::hasOTHER)
         .method("hasOVERLAP", &vcfreader::hasOVERLAP)
         .method("nsamples", &vcfreader::nsamples)
+        .method("samples", &vcfreader::samples)
         .method("header", &vcfreader::header)
         .method("string", &vcfreader::string)
         .method("line", &vcfreader::line)