这是我之前关于椭圆的一个问题的后续内容。我期待找到椭圆的半长轴和半短轴以及椭圆体（3维）的轴。我查看了其他帖子。最相关的是How to obtain the lengths of semi axes of an ellipse? in R 我收到一条混合信息。半轴的平方是形状矩阵的特征值乘以平均半径？例如继承人我的代码。战神a，b和c是半主轴吗？
library(MASS)
set.seed(1234)
m<-1
x<-NULL
k<-20
Sigma2 <- matrix(c(.58,.49,.37,.49,.58,.38,.37,.38,.34),3,3)
eta<-mvrnorm(k, mu=c(-0.85,-2.09,-1.90), Sigma2)
p1<-exp(eta[,1])/(1+exp(eta[,1]))
p2<-exp(eta[,2])/(1+exp(eta[,2]))
p3<-exp(eta[,3])/(1+exp(eta[,3]))
p<-cbind(p1,p2,p3)
n<-60
x1<-replicate(m,rbinom(k,n,p1))
x2<-replicate(m,rbinom(k,n,p2))
x3<-replicate(m,rbinom(k,n,p3))
x<-cbind(x1,x2,x3)
rate<-x/60
sigma<-var(rate)
library(cluster)
r <- ellipsoidhull(rate)
e <- sqrt(eigen(r$cov)$values)
a <- sqrt(r$d2) * e[1] # semi-principal axes
b <- sqrt(r$d2) * e[2]
c <- sqrt(r$d2) * e[3]