我想做一个更简单的测试,用于两个有两种不同处理的社区之间的比较.Ddata矩阵包含植物物种作为列和行是Year_Replication_Fertilizer(因为这是一个多年的复制研究和两种肥料处理)。我需要知道如何在simper函数中指定处理。
我使用了这段代码:
test<-read.csv("Fertilizer.csv") #Remove site names
m<-as.matrix(test[,-1]) #Use site names as row names
row.names(m)<- test[,1] #install.packages("reshape2")
#install.packages("ggplot2")
library(vegan)
library(reshape2)
library(ggplot2)
simpcom<- simper(comm=m, group=m$Group).
在原始数据框中,列包含处理,年份称为组。但我想我用我给出的代码删除了该组。在siper测试中,需要在治疗之间和之间进行比较,以查看社区的最大贡献者以及他们彼此之间的差异。请帮助我在我发送的代码的最后一行中如何做到这一点。