dirPath <- system.file("extdata/Sogin_2006", package="OTUbase")
## read in data into OTUset object
soginOTU <- readOTUset(dirPath=dirPath, level="0.03", samplefile="sogin.groups", fastafile="sogin.fasta", otufile="sogin.unique.filter.fn.list", sampleADF="sample_metadata.txt")
## get the number of OTUs
notus(soginOTU)
## get the number of samples
nsamples(soginOTU)
Run the code above in your browser using DataLab