diff --git a/R/vg.R b/R/vg.R
index c1a32e8..4911c34 100644
--- a/R/vg.R
+++ b/R/vg.R
@@ -118,7 +118,7 @@ find_vg <- function(dat,
 #' 
 compute_vg_stats <- function(dat,
                              verbose = FALSE) {
-  
+  scaled_data = dat
   if(verbose) { cat("Computing scaled means\n") }
   means <- gene_means(scaled_data)