3

以前にも同様の質問を投稿しました。ポイントが楕円内にあるかどうかを判断しようとしていました。基本的に、二変量の正規データを生成し、楕円を作成します。私が使用するコードはこちら

 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]) が上の楕円の領域内にあるかどうかを確認する必要があります。

4

2 に答える 2