【发布时间】:2015-03-02 20:10:25
【问题描述】:
下面的代码 sn-p 创建二元正态数据,绘制它,然后沿着 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轴之间的角度是多少?然后调整窗口大小,就像你在这个问题中所做的那样。现在是什么角度? -
相关:stackoverflow.com/questions/4587794/…,stackoverflow.com/questions/6322603/…;
MASS:eqscplot可能有用。
标签: r