Dear Laura
If you use multiSpeciesTP() or multiModelTP() to calculate trophic position, then the object you use to store the output has a "TPs" object which is a list. You just have to combine the output of your different TPs calculated with different TDFs and compare them as usual with a logical test (< or <=, > or >=) or the Bhattacharyya coefficient (bhat).
For example
object1 <- multiSpeciesTP()
object2 <- multiSpeciesTP()
combined_list <- do.call(c, list(object1$TPs, object2$TPs))
pairwiseComparisons(combined_list, "bhat")
If you have any questions, please let me know.
Best wishes,
Claudio