以前にも同様の質問を投稿しました。ポイントが楕円内にあるかどうかを判断しようとしていました。基本的に、二変量の正規データを生成し、楕円を作成します。私が使用するコードはこちら
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]) が上の楕円の領域内にあるかどうかを確認する必要があります。