I have been using the vegdist to run the nmds with no issues, but now I wanted to try for betapart to get the turnover and nestedness. So I tried this code below and got the following error when I tried to run the scores code to then plug into the ggplot: Error in x$species[, choices, drop = FALSE] : incorrect number of dimensions
I doubled checked my objects for the right numbers, and they looked clear. Not sure what I'm missing. Any help is greatly appreciated.
"# Shannon diversity and evenness"
"matrix_ten <- build_community_data(all_species, plots_tens, tens)"
"# Tens"
"dist_plots_10 <- beta.pair(matrix_ten, index.family = "Sorensen")"
"dist10_sim<- dist_plots_10$beta.sim"
"nmds_plots_10 <- metaMDS(dist10_sim, k=3, try=200, trace=TRUE)"
"nmds_plots_scores_10 <- plots_tens %>%"
"bind_cols(NMDS1 = scores(nmds_plots_10)[,1], NMDS2 = scores(nmds_plots_10)[,2])"
