这些天，我试图在纯素食主义者的包装上进行db-RDA分析。我有24个解释变量和46个样本。问题是，在开发db-RDA时，我只会获得前7个环境变量的结果，而与它们的排列顺序无关。可能是什么原因？ 我用代码解释了问题
＃1我向素食主义者收取费用
library(vegan)
＃2我上传了OTU数据帧，并采用了正确的格式
featuretable16 <-read.table(file = '16Sfeaturetable.txt',sep = '\t',header = TRUE)
featuretable162 <- featuretable16[,2:57]
row.names(featuretable162) <- featuretable16$Rep
ft16stransponed <- t (featuretable162)
＃3我上传了环境条件，并采用了正确的格式（解释性变量）
ambiental16 <-read.table(file = 'environmental16.txt',header = TRUE)
rownames (ambiental16) = ambiental16$SampleID
ambiental162 <- ambiental16[,2:24]
＃4进行dbrda分析和绘图，仅出现前7个解释变量
dbrda <- dbrda(ft16stransponed ~ .,ambiental162,dist="bray")
plot(dbrda)
＃5我执行了方差分析。尽管分析评估了所有环境因素，但仅显示了前七个解释变量的结果。
anova(dbrda,by="terms",permu=999)
Permutation test for dbrda under reduced model
Terms added sequentially (first to last)
Permutation: free
Number of permutations: 999
Model: dbrda(formula = ft16stransponed ~ time + Location + pH + DO + Conductivity + Temperature + Water.Flow + Water.depth + Nitrogen + Nitrate + Ammonium + Phosphorus + TOC + Amoxicillin + Azithromycin + Ciprofloxacin + Clarithromycin + Erythromycin + Lincomicina + Metronidazol + Sulfametoxazol + Ofloxacino + Trimetoprima,data = ambiental162,distance = "bray")
Df SumOfSqs F Pr(>F)
time 1 0.6731 3.5504 0.004 **
Location 1 4.7222 24.9074 0.001 ***
pH 1 0.9193 4.8488 0.002 **
DO 1 0.6811 3.5926 0.004 **
Conductivity 1 0.5441 2.8700 0.004 **
Temperature 1 0.6336 3.3418 0.004 **
Water.Flow 1 0.9761 5.1486 0.001 ***
Residual 48 9.1003
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
我也做了R笔记本，但是我不知道如何在这里提交，但是可能会更清楚。