Mam ten zestaw współrzędnych xiy:

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

i chcę obliczyć elipsoidę zawierającą ten zestaw punktów, używam funkcjiellipsoidhull() w pakiecie „klaster” i otrzymuję:

> 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

Jednak nie jest dla mnie oczywiste, jak mogę uzyskać te wyniki, długości pół-głównych osi tej elipsy.

Dowolny pomysł?

Z góry bardzo dziękuję.

Tina.