I just tried out the rarecurve function in R for one of my samples, and it worked perfectly- however, I have 72 samples to include in the curve, but I can't find a way to run them all in the code and make it work. Here is the script I am using- I presume the samples would all go in the first square brackets, but just typing them in like I did for the single one I used as a test doesn't work.

rarecurve(t(obj1.1$data$OTU_counts_bulk_barcoding[ "34b02"]), step = 20, sample = min(colSums(obj1.1$data$OTU_counts_bulk_barcoding[, Bulk_Barcoding_Sample_data$`Sample ID`])), col = "blue", cex = 1.5)