Я использую R (и пакет CCA) и пытаюсь выполнить регуляризованный канонический корреляционный анализ с двумя переменными наборами (обилие видов и количество продуктов питания, хранящихся в виде двух матриц Y и X соответственно), в которых количество единиц (N = 15) меньше числа переменных в матрицах, что составляет > 400 (большинство из них являются потенциальными "пояснительными" переменными с 12-13 "ответными" переменными). Gonzalez et al. (2008, http://www.jstatsoft.org/v23/i12/paper) обратите внимание, что пакет "включает в себя упорядоченную версию CCA для обработки наборов данных с большим количеством переменных, чем единиц", что конечно, то, что у меня есть только с 15 "единицами". Таким образом, я пытаюсь выполнить упорядоченный канонический корреляционный анализ с использованием пакета CCA, чтобы посмотреть на отношения в моих наборах переменных. Я слежу за процессом Gonzalez и др. (2008) в своей статье. Однако я получаю сообщение об ошибке Error in chol.default(Bmat) : the leading minor of order 12 is not positive definite
, и я не знаю, что это значит или что с этим делать. Вот код, и любые идеи или знания по этому вопросу будут оценены.
library(CCA)
correl <- matcor(X, Y)
img.matcor(correl, type = 2)
res.regul <- estim.regul(X, Y, plt = TRUE,
grid1 = seq(0.0001, 0.2, l=51),
grid2 = seq(0, 0.2, l=51))
Error in chol.default(Bmat) : the leading minor of order 12 is not positive definite
(Примечание: estim.regul()
занимает много времени (~ 30-40 минут), чтобы завершить, когда вы используете образец данных из CCA).
Любые советы? Кто-нибудь знает, что делать с этой ошибкой? Это потому, что у некоторых моих столбцов есть NA? Может ли это быть из-за столбцов со слишком большим количеством 0? Заранее благодарим за любую помощь, которую вы можете предложить этой комбинированной статистике и R novice.