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。