下面的代码段创建二元正态的数据,绘制它然后与线Y = X当图形窗口不是方形时,为什么椭圆会改变方向?
require(ellipse); require(MASS)
mu <- c(30,30)
Sigma <- matrix(c(900,630,630,900),2,2,byrow=TRUE)
dt <- data.frame(mvrnorm(n=1000,mu,Sigma))
names(dt) <- c("x","y")
plot(dt$x,dt$y)
df_ell <- data.frame(ellipse(cor(dt$x, dt$y), scale=c(sd(dt$x),sd(dt$y)), centre=c(mean(dt$x),mean(dt$y))))
lines(df_ell)
abline(a=0,b=1)
的直线y = x应通过适当的协方差结构的椭圆的长轴沿着绘制椭圆和平等的手段。
在一个正方形图形窗口,一切似乎很动听:
但是,如果调整窗口的大小,使其不再是方形,椭圆形,似乎自己从线定向远:
是什么导致了这种情况的发生?
绘制'y = x'线。线和'x'轴之间的角度是多少?然后调整窗口大小,就像你在这个问题中所做的那样。现在是什么角度? – nicola 2015-03-02 20:25:19
相关:http://stackoverflow.com/questions/4587794/drawing-non-intersecting-circles,http://stackoverflow.com/questions/6322603/trying-to-draw-circles-based-on-distance-between点; 'MASS:eqscplot'可能会有用。 – 2015-03-02 21:19:13