testlist <- list(rate = 9.16734271117862e-18, tab = structure(c(2.07040287343983e-172, 3.68210697042526e-69, 4.12614494748254e-57, 3394740484192042, 5.45427831416536e+168, 7.9643700123634e+78, 1.6797026028097e-79, 1.65570446607059e-290, 7.20894415361623e-167, 3.72047100402077e-255, 3.21710596525932e-160, 9.2848922257689e+294, 2.55348323297481e-228, 1.81305341023638e-301, 9.33317725725485e-82, 4.18219016778294e+110, 3.91066232231884e+152, 2.92539749011086e-56, 1.63099155258338e+21, 1.3450241930386e-95), .Dim = 4:5)) result <- do.call(TransPhylo:::probSubtree,testlist) str(result)