信頼楕円の描画に役立ちます

車のパッケージの楕円関数は、radiusパラメーターを要求します。ヘルプでそれは「楕円を生成する円の半径」であると言います。これはどのサークルか教えていただけませんか？

どうもありがとうございました

r confidence-interval multivariate-analysis

私が考えることが楕円形である場合

は、G半径として与えられるべきです。

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

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

— George Dontas

ジョージは正しいです。

— 暗がり

楕円は、任意の円のアフィン画像としてパラメーター化できます。単位円を考えると：

x = a \cos (t)

$x=a \cos (t)$

y = b \sin (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, ...)

ellipse関数は、円の中心と半径、およびアフィン変換のパラメーターを指定することと同じである共分散行列を要求することに注意してください。

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

carパッケージ関数を見てみましょう：

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

radiusあなたが直接、共分散行列を使用したい場合は、パラメータを1に設定することができますshapeパラメータ。これは、人々が望むなら正規化された行列を代わりに使用するのを助けるために導入されたと思います。

編集：whuberのコメントで述べたように、以下の2つの楕円は同じです。

> 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)

— 中華鍋
ソース

実際には、単位円から楕円へのアフィン変換は、「形状」パラメーターの平方根によって与えられます。たとえば、行列{{4,0}、{0,4}}は、円のサイズを2倍にするだけです。

— whuber