来自贾克卡指数的P值

Question

我的表有40个原始列和4列，其中4列第一列属于一个组，其余列构成另一个组。 使用以下命令计算jaccard的索引

x <- read.csv(file name,header=T, sep= )
jac <- vegdist(x,method="jaccard")

1. 从这个文件（jac）中我怎样才能找到两组的p值？
2. 我如何绘制这两组的缺口框图？

当我使用boxes(as.matrix(jac)~x$first column,notch=TRUE)时 它显示了40个箱形图。为什么会这样？