如何在具有2个嵌套因子水平的R中执行PERMANOVA？

Question

我有一个实验，比较两种鱼的饮食性能。每种饮食被随机分配到海水中的三重网箱中，每个网箱包含90条鱼。收集肠道粘膜和鱼的内容物并进行测序，以分析肠道菌群。目的是评估饮食和样本来源（粘膜对内容物）的影响及其相互作用。鱼和网围是随机效应，鱼嵌套在网围内，网围嵌套在饮食内。以下是我尝试使用 vegan 包中的df = pd.DataFrame.from_records([list1, list2]).T.fillna('') # df.columns = ['list1', 'list2'] 0 1 0 new book revealing potentially embarrassing 1 library for topic modelling several international 2 daily newspaper 函数运行受限PERMANOVA的尝试，但是我不确定是否正确。我是否正确定义了排列方案？运行PERMANOVA的公式正确吗？

adonis()

^{由reprex package（v0.3.0）于2019-09-05创建}