ошибка при записи файла после анализа адониса - PullRequest
1 голос
/ 12 июля 2020

Я пытаюсь записать результат анализа перестановки, используя приведенный ниже код, но он показывает ошибку, не могли бы вы предложить, как я могу решить эту ошибку?

R code

###adonis_calculation
BC <- phyloseq::distance(phyloseq_prop, "bray")
pmv =adonis(BC ~Season*Region, as(sample_data(physeqN3), "data.frame"),permutations = 5000)
pmv

write.table(pmv, "permutation.16s.Root.txt")

Error:

              Df SumsOfSqs MeanSqs F.Model      R2 Pr(>F)
Season          1    1.2435 1.24348  17.229 0.07405  2e-04 ***
Region          2    5.3491 2.67457  37.058 0.31855  2e-04 ***
Season:Region   2    1.3222 0.66110   9.160 0.07874  2e-04 ***
Residuals     123    8.8773 0.07217         0.52866
Total         128   16.7921                 1.00000
---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) :
  cannot coerce class ‘"adonis"’ to a data.frame
Calls: write.table ... data.frame -> as.data.frame -> as.data.frame.default
Execution halted

Большое спасибо

1 Ответ

0 голосов
/ 12 июля 2020

Вывод не в формате data.frame, можно записать вывод с помощью capture.output

capture.output(pmv, file = "permutation.16s.Root.txt")

Используя небольшой воспроизводимый пример

library(vegan)
data(dune)
data(dune.env)
pmv <- adonis(dune ~ Management*A1, data = dune.env)
capture.output(pmv, file  = "/Users/akrun/tmp/permutation.16s.Root.txt")
 

-выход

введите описание изображения здесь

...