小编use*_*215的帖子

我正在尝试找到一个点位于椭圆内的概率？例如,如果我在95%椭圆体区域中绘制300个数据集的双变量数据(x,y),我如何计算300个点中椭圆内部的多少次？

继承人我正在使用的代码

   library(MASS)
   seed<-1234
   x<-NULL
   k<-1
   Sigma2 <- matrix(c(.72,.57,.57,.46),2,2)
   Sigma2
   rho <- Sigma2[1,2]/sqrt(Sigma2[1,1]*Sigma2[2,2])
   rho
   eta1<-replicate(300,mvrnorm(k, mu=c(-1.59,-2.44), Sigma2))

   library(car)
   dataEllipse(eta1[1,],eta1[2,], levels=c(0.05, 0.95))

谢谢你的帮助.

我之前发过一个类似的问题.我试图确定一个点是否位于椭圆内.基本上我生成一些双变量法线数据并创建一个椭圆.继承我使用的代码

 library(MASS)
 set.seed(1234)
 x1<-NULL
 x2<-NULL
 k<-1
 Sigma2 <- matrix(c(.72,.57,.57,.46),2,2)
 Sigma2
 rho <- Sigma2[1,2]/sqrt(Sigma2[1,1]*Sigma2[2,2])

 eta<-replicate(300,mvrnorm(k, mu=c(-2.503,-1.632), Sigma2)) 

 p1<-exp(eta)/(1+exp(eta))
 n<-60
 x1<-replicate(300,rbinom(k,n,p1[,1]))
 x2<-replicate(300,rbinom(k,n,p1[,2]))

 rate1<-x1/60
 rate2<-x2/60

 library(car)
 dataEllipse(rate1,rate2,levels=c(0.05, 0.95)) 

我需要找出对(p1 [,1],p1 [,2])是否位于上面椭圆的区域内.