Without a reprex (see the FAQ), all I can offer is the example from the help page
library("indicspecies")
#> Loading required package: permute
## Loads species data
data(wetland)
## Creates three clusters using kmeans
wetkm <- kmeans(wetland, centers=3)
## Runs the combination analysis using IndVal.g as statistic
wetpt <- multipatt(wetland, wetkm$cluster)
## Lists those species with significant association to one combination
summary(wetpt)
#>
#> Multilevel pattern analysis
#> ---------------------------
#>
#> Association function: IndVal.g
#> Significance level (alpha): 0.05
#>
#> Total number of species: 33
#> Selected number of species: 5
#> Number of species associated to 1 group: 3
#> Number of species associated to 2 groups: 2
#>
#> List of species associated to each combination:
#>
#> Group 1 #sps. 2
#> stat p.value
#> Psespi 0.951 0.005 **
#> Elesp. 0.716 0.020 *
#>
#> Group 3 #sps. 1
#> stat p.value
#> Phynod 0.771 0.005 **
#>
#> Group 1+2 #sps. 2
#> stat p.value
#> Orysp. 0.960 0.005 **
#> Ludads 0.882 0.005 **
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Created on 2023-07-29 with reprex v2.0.2