S?rensen 指數(shù)和 Jaccard 指數(shù)作為配對樣方間β 多樣性的度量，使用 Baselga 在 2010 和 2012 年提出的 β 多樣性分解方法結(jié)合 Legendre 在 2014 年對該方法的總結(jié)和改進，將群落間物種組成差異（即 β 多樣性）分解成周轉(zhuǎn)組分（用 Repl 表示）和嵌套組分（用 Nest 表示）（余等，2022）。

Jaccard指數(shù)計算：

library(adespatial)
#數(shù)據(jù)為豐度數(shù)據(jù)（行：樣本，列：物種）
data=read.csv("asv_table.csv",row.names = 1)
#利用 beta.div.comp函數(shù)計算
out1 <- beta.div.comp(t(data), coef ='J', quant = FALSE)
#樣本對之間的相異指數(shù)計算
D <- as.matrix(out1$D)
diag(D) <- "NA" 
D[upper.tri(D)] <- "NA"
D1=melt(D)
D1=subset(D1,!D1$value=="NA")
#樣本對之間的物種替換（repl）
repl <- as.matrix(out1$repl)
diag(repl) <- "NA" 
repl[upper.tri(repl)] <- "NA"
repl1=melt(repl)
repl1=subset(repl1,!repl1$value=="NA")
#樣本對之間的豐富度差異（nest）
nest <- as.matrix(out1$nest)
diag(nest) <- "NA" 
nest[upper.tri(nest)] <- "NA"
nest1=melt(nest)
nest1=subset(nest1,!nest1$value=="NA")
#合并數(shù)據(jù)nest-repl至data.frame：
dff=merge(D1,repl1,by=c("Var1","Var2"))
dff=merge(dff,nest1,by=c("Var1","Var2"))
colnames(dff)=c("site1","site2","D","repl","nest")
dff$site=paste(dff$site1,dff$site2,sep="_")

library(ggtern)
dff$D=as.numeric(dff$D)
dff$repl=as.numeric(dff$repl)
dff$nest=as.numeric(dff$nest)
#計算平均值
dm=as.data.frame(t(data.frame(mean=colMeans(dff[,3:5]))))
#繪圖
p=ggtern(data=dff,aes(1-D,repl,nest))+
  geom_point(alpha=0.2,size=3)+
  stat_density_tern(
    geom='polygon',
    base = "ilr",
    aes(fill=..level..,
        alpha = ..level..))+
  geom_point(data=dm,col="green4",size=5)+
  scale_fill_gradient(low = "blue",
                      high = "red")+
  labs(x ="1-β",y = "Repl",z = "Nest",
       title = "",fill = "level")+
  guides(alpha = "none",position="bottom")+
  theme_rgbw()+
  theme(legend.position = "bottom",
        legend.key = element_rect(fill = NA),
        text=element_text(size=16,  family="serif"))
p
ggsave(filename="三元相圖.pdf",width = 10,height = 10)
#輸出表格
write.csv(dff,"repl_nset.csv",row.names = F)

如在余等（2022）的文章中，Jaccard 指數(shù)法顯示在灌木植物中，周轉(zhuǎn)和嵌套對 β 多樣性的貢獻分別為 82%和 18%；在草本植物中，周轉(zhuǎn)和嵌套對 β 多樣性的貢獻分別為 95%和 5%。
有需要也可@TB：西北工作室