Извлечь p-значения и тесты `mvabund::anova.manyglm` в виде одной таблицы

Я хотел бы получить раздел «Одномерные тесты» из результата anova.manyglm, но атрибуты по умолчанию дают две отдельные таблицы, поэтому я пытаюсь объединить их, чтобы тест и столбцы p-значения одной и той же переменной оставались рядом друг с другом. .

library(mvabund)
data("Tasmania")
attach(Tasmania)
tasmvabund=mvabund(copepods)
tas_mod = manyglm(tasmvabund~ block*treatment)
tas_sum <- summary.manyglm(tas_mod)
tas_aov <- anova.manyglm(tas_mod, resamp = "perm.resid",  p.uni= "adjusted", nBoot = 50)



  > str(tas_aov)
List of 11
 $ family      : chr "negative.binomial"
 $ p.uni       : chr "adjusted"
 $ test        : chr "Dev"
 $ cor.type    : chr "I"
 $ resamp      : chr "perm.resid"
 $ nBoot       : num 50
 $ shrink.param: num [1:4] 0 0 0 0
 $ n.bootsdone : num 50
 $ table       :'data.frame':   4 obs. of  4 variables:
  ..$ Res.Df  : int [1:4] 15 12 11 8
  ..$ Df.diff : num [1:4] NA 3 1 3
  ..$ Dev     : num [1:4] NA 117.5 66.9 37.4
  ..$ Pr(>Dev): num [1:4] NA 0.0196 0.0196 0.0588
  ..- attr(*, "heading")= chr [1:2] "Analysis of Deviance Table\n" "Model: tasmvabund ~ block * treatment"
  ..- attr(*, "title")= chr "\nMultivariate test:\n"
 $ uni.p       : num [1:4, 1:12] NA 0.4706 0.0196 0.451 NA ...
  ..- attr(*, "title")= chr "Univariate Tests:"
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:4] "(Intercept)" "block" "treatment" "block:treatment"
  .. ..$ : chr [1:12] "Ameira" "Adopsyllus" "Ectinosoma" "Ectinosomat" ...
 $ uni.test    : num [1:4, 1:12] NA 4.93 13.94 6.35 NA ...
  ..- attr(*, "title")= chr "Univariate Tests:"
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : chr [1:4] "(Intercept)" "block" "treatment" "block:treatment"
  .. ..$ : chr [1:12] "Ameira" "Adopsyllus" "Ectinosoma" "Ectinosomat" ...
 - attr(*, "class")= chr "anova.manyglm"

При этом печатается только многомерный раздел.

tas_aov$table
'data.frame':   4 obs. of  4 variables:
 $ Res.Df  : int  15 12 11 8
 $ Df.diff : num  NA 3 1 3
 $ Dev     : num  NA 117.5 66.9 37.4
 $ Pr(>Dev): num  NA 0.0196 0.0196 0.0588
 - attr(*, "heading")= chr [1:2] "Analysis of Deviance Table\n" "Model: tasmvabund ~ block * treatment"
 - attr(*, "title")= chr "\nMultivariate test:\n"

Это печатает только p-значения одномерных тестов.

tas_aov$uni.p
 str(tas_aov$uni.p)
 num [1:4, 1:12] NA 0.4706 0.0196 0.451 NA ...
 - attr(*, "title")= chr "Univariate Tests:"
 - attr(*, "dimnames")=List of 2
  ..$ : chr [1:4] "(Intercept)" "block" "treatment" "block:treatment"
  ..$ : chr [1:12] "Ameira" "Adopsyllus" "Ectinosoma" "Ectinosomat" ...

И это печатает только тестовую таблицу.

tas_aov$uni.test
 str(tas_aov$uni.test)
 num [1:4, 1:12] NA 4.93 13.94 6.35 NA ...
 - attr(*, "title")= chr "Univariate Tests:"
 - attr(*, "dimnames")=List of 2
  ..$ : chr [1:4] "(Intercept)" "block" "treatment" "block:treatment"
  ..$ : chr [1:12] "Ameira" "Adopsyllus" "Ectinosoma" "Ectinosomat" ...

Итак, я пытаюсь объединить их, а затем упорядочить столбцы по именам, но пока мне это не удалось.

tab <- cbind(tas_aov$uni.p, tas_aov$uni.test)

tab[ , order(names(tab))]
Error in order(names(tab)) : argument 1 is not a vector

В идеале мне нужен только исходный раздел Univariate, потому что в нем аккуратно напечатаны имена столбцов, т. Е. Тест рядом с p-значениями, но я доволен просто объединением таблиц. Большое спасибо за любые подсказки.