I am trying to plot species abundance for intestinal bacteria analysed using a 16S RNA sequencing method. I am using RStudio (1.1.456) and R 3.5.1 My input is a file called "Speciespercentagestest.csv", here. I was following the tutorial on "Importing phyloseq Data", more specifically, loading the table into phyloseq using the sample_data command. Here is my code:

BiocManager::install('phyloseq')
library(readr)
library(phyloseq)
library(ggplot2)
species_analysis <- read_csv("Species_percentages_test.csv")
theme_set(theme_bw())
species_p <- sample_data(species_analysis, errorIfNULL=TRUE)
taxa_are_rows(species_p)
plot_bar(species_p)

My issue is with the "taxaarerows" function. I believe I need to use it to specify the orientation of operational taxonomic units (OTUs), in my case "species" but I am none the wiser on how to apply this function in my specific example. I am sure this is a rookie question but this is the first time I deal with microbiome data. Your help is truly appreciated!