gpt4 book ai didi

r - 如何获得椭圆的半轴长度?在 R

转载 作者:行者123 更新时间:2023-12-04 11:42:28 25 4
gpt4 key购买 nike

我有这组 x 和 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

我想计算包含这组点的椭球,我使用函数 ellipsoidhull()在包“集群”中,我得到:
> 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

然而,我不清楚如何从这些结果中得到这个椭圆的半长轴的长度。

任何的想法?

非常感谢您提前。

蒂娜。

最佳答案

半轴的平方是
形状矩阵的特征值乘以平均平方半径。

x <- c(1.798805,2.402390,2.000000,3.000000,1.000000)
y <- c(0.3130147,0.4739707,0.2000000,0.8000000,0.1000000)
d <- cbind( x, y )
library(cluster)
r <- ellipsoidhull(d)
plot( x, y, asp=1, xlim=c(0,4) )
lines( predict(r) )
e <- sqrt(eigen(r$cov)$values)
a <- sqrt(r$d2) * e[1] # semi-major axis
b <- sqrt(r$d2) * e[2] # semi-minor axis
theta <- seq(0, 2*pi, length=200)
lines( r$loc[1] + a * cos(theta), r$loc[2] + a * sin(theta) )
lines( r$loc[1] + b * cos(theta), r$loc[2] + b * sin(theta) )

关于r - 如何获得椭圆的半轴长度?在 R,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/18278382/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com