эта ссылка

аю простую задачу оптимизации. Набор данных имеет 26 столбцов и более 3000 строк. Исходный код выглядит так

Means   <- colMeans(Returns)
Sigma   <- cov(Returns)
invSigma1 <- solve(Sigma)

И все работает отлично, но затем я хочу сделать то же самое для более короткого периода (только 261 строк) и функция решения пишет следующую ошибку:

solve(Sigma)
Error in solve.default(Sigma) : 
  Lapack routine dgesv: system is exactly singular 

Это странно, потому что, когда я делаю то же самое с некоторыми случайными числами:

Returns<-matrix(runif(6786,-1,1), nrow=261)
Means   <- colMeans(Returns)
Sigma   <- cov(Returns)
invSigma <- solve(Sigma)

ошибка не возникает вообще. Может ли кто-нибудь объяснить мне, где может быть проблема и как ее решить. Спасибо большое, Алекс

Ответы на вопрос(3)

Ваш ответ на вопрос