path <- system.file("extdata", "virushunter.tsv", package = "Virusparies")
vh_file <- ImportVirusTable(path)
vh_file_filtered <- VhgPreprocessTaxa(vh_file,"Family")
processed_taxa <- VhgAddPhylum(vh_file_filtered,"ViralRefSeq_taxonomy")
print(unique(processed_taxa$Phylum))
Run the code above in your browser using DataLab