seq <- readRNAStringSet(system.file("examples/RNA_example.fasta",
package="rRDP"))
### decode the actual classification
actual <- decode_Greengenes(names(seq))
### use RDP to predict the classification
pred <- predict(rdp(), seq)
### calculate accuracy
confusionTable(actual, pred, "genus")
accuracy(actual, pred, "genus")
Run the code above in your browser using DataLab