У меня есть проблема в R. Я делаю несколько тестов NMDS с использованием пакета vegan, и пока он идет довольно хорошо, пока я не попытался использовать этот код для рисования полигонов в моем NMDS:
ordiplot(nmds,type = "n")
colors<-c(rep("blue",22),rep("red",24))
treat<-c(rep("ngs",22),rep("gs",24))
for(i in unique(treat)) {
ordihull(nmds$points[grep(i,treat),],groups = treat[treat==i],draw = "polygon",col = colors[grep(i,treat)],label = F)
}
При этом появляется сообщение об ошибке:
Error in chull(X) : finite coordinates are needed
In addition: Warning messages:
1: In complete.cases(pts) & !is.na(groups) :
longer object length is not a multiple of shorter object length
2: In groups == is & kk :
longer object length is not a multiple of shorter object length
Честно говоря, я не знаю, чем это вызвано.Я нашел кого-то, у кого была такая же проблема, и было высказано предположение, что в некоторых из них было мало наблюдений.Я сомневаюсь, что это так, но я не могу отклонить это просто потому, что не знаю, как это проверить.Мои данные содержат 46 наблюдений и две обработки.Вот данные, лечение меняется при наблюдении 23 введите здесь ссылку описания