Я провожу анализ PGLS (phylogeneti c обобщенные наименьшие квадраты) 5 независимых переменных (AE) по сравнению с зависимой переменной (X) с использованием функции phylolm в R.
вот строка кода, который я использую:
X <- seq(1,300,1)
model_fit <- phylolm(X~factor(data$A)+factor(data$B)+factor(data$C)+factor(data$D)+factor(data$E),...
data=data, phy=pruned.tree, model="lambda")'
pruned.tree - это дерево филогенетических данных c, которое было сокращено, чтобы соответствовать названиям видов, представляющим каждое значение X.
множителями AE являются категориальные присвоения различных экологических переменных (символов / ячеек), которые варьируются от 2 до 4 категорий на фактор.
Когда я запускаю все факторы вместе, я получаю следующую ошибку:
Error in solve.default(comp$XX) :
Lapack routine dgesv: system is exactly singular: U[186,186] = 0
This ошибка не имеет смысла для меня (или указать, какая часть моей матрицы данных вызывает проблемы). Я могу успешно запустить филолм, используя различные комбинации меньшего количества факторов, но не все вместе. Что может быть не так? Спасибо за помощь!