问题描述
这些天,我试图在纯素食主义者的包装上进行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笔记本,但是我不知道如何在这里提交,但是可能会更清楚。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)