data(wetland) ## Load species data
wetkm = kmeans(wetland, centers=3)## Create three clusters using
strassoc(wetland, wetkm$cluster, func="IndVal.g") ## Compute Dufrene and Legendre's IndVal
strassoc(wetland, wetkm$cluster, func="r", nboot =100) ## Compute point-biserial correlation, with bootstrap 95 percent confidence intervals
Run the code above in your browser using DataLab