Я пытаюсь выполнить анализ RDA для данных моего сообщества. Я хочу сравнить численность улиток 9 видов с 6 абиотическими и 9 биотическими параметрами на 6 участках. Количество переменных превышает количество сайтов, и я чувствую, что это проблема для моего анализа. Когда я суммирую свой RDA, я вижу только 5 осей RDA и первые 5 переменных image . Эта проблема не присутствует для численности улитки. Когда я проверяю R² и корректирую значения R², я получаю 1 и NA в качестве вывода. Когда я запускаю тот же код в других наборах данных, содержащих больше сайтов, чем переменных, я получаю ожидаемый вывод всех перечисленных переменных в файле среды и «нормальные» значения R².
Я использую пакет 'vegan' и команду rda () в R.
Итак, выполнение моего кода на других наборах данных работает, и я думаю, что это как-то связано с тем, что в моем наборе данных больше переменных, чем сайтов? Как я мог решить эту проблему?
physico = read.csv("site16.csv",row.names=1, sep = ';')
snail = read.csv("snails1_6.csv",row.names=1, sep = ';')
phys.rda <- rda(snail~.,physico)
summary(phys.rda)
RsquareAdj(phys.rda)$r.squared # R^2
RsquareAdj(phys.rda)$adj.r.squared # adjusted R^2
Я ожидаю, что все мои переменные будут перечислены в сводке вместе с большим количеством осей RDA вместо вывода, показанного на рисунке.