From 89d5fd0bf412bd9635c4c54327c2e522ab49e500 Mon Sep 17 00:00:00 2001
From: David Gerard <gerard.1787@gmail.com>
Date: Thu, 3 Aug 2023 20:20:39 -0400
Subject: [PATCH] export_vcf, but still not exported

---
 R/multidog.R | 8 +++++---
 1 file changed, 5 insertions(+), 3 deletions(-)

diff --git a/R/multidog.R b/R/multidog.R
index b2497e7..e520098 100644
--- a/R/multidog.R
+++ b/R/multidog.R
@@ -690,6 +690,9 @@ export_vcf <- function(obj, filename) {
       requireNamespace("S4Vectors", quietly = TRUE) &&
       requireNamespace("IRanges", quietly = TRUE)) {
 
+    ## Get sample neames
+    indvec <- colnames(format_multidog(x = obj, varname = "postmean"))
+
     ploidy <- unique(obj$snpdf$ploidy)
     stopifnot(length(ploidy) == 1)
     obj$inddf$alt <- obj$inddf$size - obj$inddf$ref
@@ -738,11 +741,11 @@ export_vcf <- function(obj, filename) {
         strand=NULL,
         seqinfo = NULL,
         names = obj$snpdf$snp),
-      colData = as(matrix(nrow = nind, ncol = 0), "DataFrame"),
+      colData = DataFrame(row.names = indvec),
       exptData = list(
         header = VariantAnnotation::VCFHeader(
           reference = character(),
-          samples = character(),
+          samples = indvec,
           header = IRanges::DataFrameList(
             fileformat = S4Vectors::DataFrame(row.names = "fileformat", Value = "VCFv4.3"),
             fileDate = S4Vectors::DataFrame(row.names = "fileDate", Value = gsub("-", "", Sys.Date())),
@@ -777,4 +780,3 @@ export_vcf <- function(obj, filename) {
                    " \"GenomicRanges\", \"S4Vectors\", \"IRanges\"))\n"))
   }
 }
-