r - 二変量楕円を使用したカバレッジの検索

Question

見積もりの​​問題で悩んでいます。いくつかの例で、二変量正規分布によって生成された点のベクトルを使用して二変量楕円を計算する方法を示しました。コードは正常に動作しますが、カバレッジ (生成されたまたは真の ps(p1,p2) が推定された楕円に含まれる回数) が非常に低いようです。また、古いバージョンの R では、新しいバージョンと比較して大幅に異なる結果が得られたことも述べておく必要があります。現在、R 3.0.1 を使用しています。問題を再現できるコードは次のとおりです。

  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(-1.01,-2.39), Sigma2)) 
  p1<-exp(eta)/(1+exp(eta)) # true p's
  n<-60
  x1<-replicate(300,rbinom(k,n,p1[1,])) 
  x2<-replicate(300,rbinom(k,n,p1[2,]))

  rate1<-x1/60  # Estimated p's
  rate2<-x2/60
  library(car)
  ell <- dataEllipse(rate1, rate2, levels=c(0.05, 0.95))
  library(sp)
  within<-point.in.polygon(p1[1,], p1[2,], ell$`0.95`[,1], ell$`0.95`[,2])
  mean(within)    # coverage