Commits

Davide Cittaro committed d314609

dump pca data into npz file

  • Participants
  • Parent commits 07f2c41

Comments (0)

Files changed (1)

         file_out.write('\t'.join(fields))
         x += 1  
   
+  if cli_options.dump:
+    sys.stderr.write("Dumping PCA data...\n")
+    ext_ind = cli_options.vcf.rfind('.vcf')
+    if ext_ind > 1:
+      dump_file = "PCAdata_%s" % cli_options.vcf[:ext_ind]
+    else:
+      dump_file = "PCAdata"
+    np.savez(dump_file, results.Wt, results.Y)