Eu tenho este conjunto de coordenadas x e y:

x<-c(1.798805,2.402390,2.000000,3.000000,1.000000)
y<-c(0.3130147,0.4739707,0.2000000,0.8000000,0.1000000)
as.matrix(cbind(x,y))->d

e eu quero calcular o elipsóide que contém esse conjunto de pontos, eu uso a funçãoellipsoidhull() no pacote "cluster", e recebo:

> ellipsoidhull(d)
'ellipsoid' in 2 dimensions:`
 center = ( 2.00108 0.36696 ); squared ave.radius d^2 =  2`
and shape matrix =
x 0.66590 0.233106
y 0.23311 0.095482
  hence, area  =  0.60406

No entanto, não é óbvio para mim como posso obter esses resultados, os comprimentos dos eixos semi-maiores dessa elipse.

Qualquer ideia?

Muito obrigado antecipadamente.

Tina.