R predict.ellipsoid 椭球体预测方法

说明

计算椭球边界上的点，主要用于绘图。

用法

predict.ellipsoid(object, n.out=201, ...)
## S3 method for class 'ellipsoid'
predict(object, n.out=201, ...)
ellipsoidPoints(A, d2, loc, n.half = 201)

参数

object	类 ellipsoid 的对象，通常来自 ellipsoidhull() ；或者任何具有适当组件的 list-like 对象，请参阅下面的详细信息。
n.out , n.half	要创建的点数的一半。
A , d2 , loc	辅助 ellipsoidPoints 的参数，见下文。
...	传入和传出方法。

细节

注意 ellipsoidPoints 是 predict.ellipsoid 的主力函数，是 ellipsoid 对象的独立函数和方法，请参阅 ellipsoidhull 。 object的类没有被检查；它必须仅具有有效组件loc(长度p)、p \times p矩阵cov(对应于A)和d2作为中心、形状(“covariance”)矩阵和平方平均半径(或距离)或qchisq(*, p) 分位数。

不幸的是，目前仅针对 p = 2 实现；非常欢迎为 p \ge 3 做出贡献。

值

维度为 2*n.out 乘以 p 的数字矩阵。

例子

 ## see also  example(ellipsoidhull)

## Robust vs. L.S. covariance matrix
set.seed(143)
x <- rt(200, df=3)
y <- 3*x + rt(200, df=2)
plot(x,y, main="non-normal data (N=200)")
mtext("with classical and robust cov.matrix ellipsoids")
X <- cbind(x,y)
C.ls <- cov(X) ; m.ls <- colMeans(X)
d2.99 <- qchisq(0.99, df = 2)
lines(ellipsoidPoints(C.ls, d2.99, loc=m.ls), col="green")
if(require(MASS)) {
  Cxy <- cov.rob(cbind(x,y))
  lines(ellipsoidPoints(Cxy$cov, d2 = d2.99, loc=Cxy$center), col="red")
}# MASS

也可以看看

ellipsoidhull、volume.ellipsoid。