Ошибка в G * t(hat): неконформируемые массивы
Пытаясь выполнить PERMANOVA, используя adonis в R. Добавленный в мою таблицу данных и факторы, тем не менее, я продолжаю получать ошибку выборки - Ошибка в G * t(hat): несоответствующие массивы
Нет данных N/A, и он работает с готовой матрицей данных и подмножеством моей собственной таблицы данных.
Пример данных, которые не работают;
subset<-matrix(c(0,0.000666667,0.001333333,0.000333333,0.000333333,0.019833333,0.007666667,0.014666667,0.0005,0.022833333,0.016833333,0.018166667,0.000666667,0.009666667,0.008833333,0.009166667,0.001333333,0.018666667,0.0295,0.031833333),ncol=4)
test_groups<-c(1,1,2,2)
adonis(subset~test_groups)
Ошибка в G * t(hat): неконформируемые массивы
Тестовые данные, которые действительно работают;
test<-matrix(data=c(0.1,0.1,0.3,0.1,0.2,0.3,0.1,0,0.3,0.1,0.2,0.3,0.1,0.2,0.3,0.1),ncol=4)
Почему это не работает?
1 ответ
Длина test_groups
вектор равен 4, а число строк subset
это 5.
Ниже я установил на 5 количество рядов subset
и adonis
Команда работает хорошо.
library(vegan)
subset <- matrix(c(0,0.000666667,0.001333333,0.000333333,0.000333333,0.019833333,0.007666667,0.014666667,0.0005,0.022833333,0.016833333,0.018166667,0.000666667,0.009666667,0.008833333,0.009166667,0.001333333,0.018666667,0.0295,0.031833333),ncol=5)
test_groups <- c(1,1,2,2)
adonis(subset~test_groups)
Вот результат.
'nperm' >= set of all permutations: complete enumeration.
Set of permutations < 'minperm'. Generating entire set.
Call:
adonis(formula = subset ~ test_groups)
Permutation: free
Number of permutations: 23
Terms added sequentially (first to last)
Df SumsOfSqs MeanSqs F.Model R2 Pr(>F)
test_groups 1 0.22804 0.22804 1.0592 0.34624 0.3333
Residuals 2 0.43059 0.21530 0.65376
Total 3 0.65864 1.00000