Ошибка построения 3D PCA: младший младший номер 3 не является положительно определенным - PullRequest
Новая
       56

Ошибка построения 3D PCA: младший младший номер 3 не является положительно определенным

0 голосов
/ 27 ноября 2018

Я пытался решить проблему в течение последних нескольких часов, но безуспешно.

Я обнаружил одно сообщение об ошибке в другом посте, похожем здесь: Печать Ellipse3d в R Plotly withэллипс поверхности

Единственный ответ подчеркивает, что это из-за нехватки образцов.Но мой размер выборки = 3;что является минимумом для создания эллипса.

Другой ответ на другой пост подчеркивает, что это может быть потому, что матрица, вводимая в ellipse3d (), не является симметричной квадратной матрицей, но на самом деле это потому, что я беруполе ввода (data_frame_of_interest) при его вводе.

Подробнее о проблеме:

Я пытаюсь построить 3D PCA, используя пакет rgl в R и prcomp ();где top_noSR3 является объектом prcomp ().Меня интересуют первые 3 ПК для следующего примера. 

#Here is the exact PCA object obtained from prcomp(). It is 
#top_noSR3$x[,c(1,2,3)].    

top_noSR3 <- list(x = structure(c(-1.51867921165966, 4.7156146538954, -0.812795773332441, 
-2.38413966890339, -2.57305046487183, 8.4906222835565, 1.66680581870767, 
-7.5843776373923, 3.35199279752431, 2.4452307290032, -7.33666903970592, 
1.53944551317846, -1.72127982564219, 2.68563509110228, -2.75720515666432, 
1.79284989120424, -0.962414945057766, -1.44498205066367, -0.270435252418625, 
2.67783224814005, -0.780985238414212, 7.3028581587737, -1.3641698814598, 
-5.15770303889969, 0.989625278145346, 0.101301456051966, 1.96037141973899, 
-3.05129815393625, -2.29377625653868, -4.3601221831666, 5.9554046825257, 
0.698493757179566, 2.3724852157055, 0.0505080025883657, -1.97013906903266, 
-0.452854149261224), .Dim = c(12L, 3L), .Dimnames = list(c("SR1.4", 
"SR1.5", "SR1.6", "SR1.7", "SR2.2", "SR2.3", "SR2.4", "SR2.5", 
"CI1.4", "CI1.5", "CI1.6", "CI1.7"), c("PC1", "PC2", "PC3"))))

#col_object is an array of colours used for colouring the points in the PCA:
col_object <- c("blue","red","green","purple","blue","red",   
"green","purple","blue","red","green","purple")

plot3d(top_noSR3$x[,c(1,2,3)], type = "s", size = 1, col = col_object)

#Divide object into 3 components for 3d ellipse plotting function

for_plot3d <- function(casetime){

  return(top_noSR3$x[casetime, c(1,2,3)])
}

BASE <- c("SR1.4", "SR2.2", "CI1.4")
EA <- c("SR1.5", "SR2.3", "CI1.5")
LA <- c("SR1.6", "SR2.4", "CI1.6")
FU <- c("SR1.7", "SR2.5", "CI1.7")

BASE_for3d <- for_plot3d(BASE)
EA_for3d <- for_plot3d(EA)
LA_for3d <- for_plot3d(LA)
FU_for3d <- for_plot3d(FU)

#The following arguments used above (BASE,EA, etc.) ^ are 
#timepoints/rownames of interest. length(BASE) = ... = length(FU) = 3


center4plot <- function(x){
  k <- mean(x[,1])
  h <- mean(x[,2])
  z <- mean(x[,3])
  return(c(k,h,z))
}

center_base <- center4plot(BASE_for3d)
center_EA <- center4plot(EA_for3d)
center_LA <- center4plot(LA_for3d)
center_FU <- center4plot(FU_for3d)

ellipse_base <- ellipse3d(cov(BASE_for3d),level = 0.95, centre=center_base)
ellipse_EA <- ellipse3d(cov(EA_for3d), level = 0.95, centre=center_EA)

#ERROR appears here, in ellipse_LA
ellipse_LA <- ellipse3d(cov(LA_for3d), level = 0.95, centre=center_LA)
ellipse_FU <- ellipse3d(cov(FU_for3d), level = 0.95, centre=center_FU)

При запуске строки с ellipse_LA я получаю следующее сообщение: ellipse3d () Ошибка в chol.default (cov): ведущий младшийПорядка 3 не является положительно определенным

Спасибо всем за ваше время и усилия.

1 Ответ

0 голосов
/ 28 ноября 2018

Когда я запускаю текущую версию кода, я получаю ошибку, которую вы видели здесь: 

> ellipse_base <- ellipse3d(cov(BASE_for3d),level = 0.95, centre=center_base)
Error in chol.default(cov) : 
  the leading minor of order 3 is not positive definite

Если я смотрю на собственные значения этой матрицы, я вижу, что сообщение верное: 

> eigen(cov(BASE_for3d))
eigen() decomposition
$values
[1] 1.454977e+01 1.433775e+00 8.521084e-16

$vectors
           [,1]       [,2]       [,3]
[1,] 0.81230516  0.5213444  0.2614581
[2,] 0.04362183  0.3927276 -0.9186197
[3,] 0.58159906 -0.7576048 -0.2962727

Обратите внимание, что наименьшее собственное значение приблизительно равно нулю.

Объект BASE_for3d имеет 3 точки.Все они лежат в плоскости (как и любые 3 точки).Таким образом, ковариационная матрица для них - ранг 2, а не ранг 3. Вам нужно увеличить BASE, чтобы включить как минимум 4 строки не все в плоскости, а не только 3. Например, 

BASE <- c("SR1.4", "SR2.2", "CI1.4", "SR1.5")

Эторезультаты следующие: 

BASE_for3d <- for_plot3d(BASE)
center_base <- center4plot(BASE_for3d)
ellipse_base <- ellipse3d(cov(BASE_for3d),level = 0.95, centre=center_base)
shade3d(ellipse_base, alpha = 0.2)

screenshot

...