Giúp vẽ hình elip tự tin

Chức năng hình elip của gói xe yêu cầu một radiustham số. Trong trợ giúp nói rằng đó là "bán kính của vòng tròn tạo ra hình elip". Bạn có thể vui lòng cho tôi biết đây là vòng tròn nào?

Cảm ơn rât nhiều

r confidence-interval multivariate-analysis

— George Dontas
nguồn

Tôi nghĩ rằng nếu ellipsoid có dạng

, thì g nên được đưa ra dưới dạng bán kính.

(x - \bar{x})^{T} S^{- 1} (x - \bar{x}) = g^{2}

$(x-\bar{x})^TS^{-1}(x-\bar{x})=g^2$

— George Dontas

George đúng.

— u ám

Một hình elip có thể được tham số hóa như hình ảnh affine của bất kỳ hình tròn nào. Nếu chúng ta xem xét vòng tròn đơn vị:

x = = một \cos (t)

$x=a \cos (t)$

y = = b tội (t)

$y=b \sin (t)$

ellipse(center, shape, radius, log="", center.pch=19, center.cex=1.5, 
  segments=51, add=TRUE, xlab="", ylab="", 
   col=palette()[2], lwd=2, fill=FALSE, fill.alpha=0.3, grid=TRUE, ...)

Bạn có thể nhận thấy ellipsehàm yêu cầu tâm và bán kính của vòng tròn, cũng như ma trận hiệp phương sai, tương đương với việc đưa ra các tham số của phép biến đổi affine.

center  2-element vector with coordinates of center of ellipse.
shape   2 * 2 shape (or covariance) matrix.
radius  radius of circle generating the ellipse.

Chúng ta hãy xem carchức năng gói:

ellipse <- t(center + radius * t(unit.circle %*% chol(shape)))

Các radiustham số có thể được thiết lập để 1 nếu bạn muốn sử dụng ma trận hiệp phương sai trực tiếp cho các shapetham số. Tôi tin rằng nó đã được giới thiệu để giúp mọi người sử dụng ma trận chuẩn hóa thay vào đó nếu họ thích như vậy.

Chỉnh sửa: Như đã đề cập trong nhận xét của người đăng ký, hai hình elip bên dưới giống nhau.

> library(car)
> s=matrix(c(1,0,0,1), nrow=2, ncol=2)
> plot(0, 0, xlim=c(-5,5), ylim=c(-5,5))
> ellipse(c(0,0), 4*s, 1)
> ellipse(c(0,0), s, 2)

— Chảo
nguồn

Trên thực tế, phép biến đổi affine từ vòng tròn đơn vị sang hình elip được cho bởi căn bậc hai của tham số "hình dạng". Ví dụ: ma trận {{4,0}, {0,4}} sẽ chỉ nhân đôi kích thước của vòng tròn.

— whuber