1
の1%分位数のための分位数推定量として
私は最近、(非統計的、適用された)論文の連続確率変数に対する次の変位値推定量を見つけました:100の長いベクトルxxx場合、1%変位値はmin(x)min(x)\min(x)推定されます。以下に、その実行方法を示します。以下は、N(0,1)分布からの100個の長いサンプルの100,000回のシミュレーション実行からのmin(x)min(x)\min(x)推定量の実現のカーネル密度プロットです。垂直線は真の値、つまりN(0,1)分布の理論上の1%分位です。シミュレーションのコードも示します。N(0,1)N(0,1)N(0,1)N(0,1)N(0,1)N(0,1) M=10e5; n=100 quantiles=rep(NA,M) for(i in 1:M){ set.seed(i); quantiles[i]=min(rnorm(n)) } plot(density(quantiles),main="Kernel density estimate of quantiles from M=100,000 simulation runs"); abline(v=qnorm(1/n)) グラフは、分布(例)については質的に似ています。どちらの場合も、推定量は下方に偏っています。ただし、他の推定量と比較しないと、それ以外の点でどれほど優れているかを言うのは困難です。したがって、私の質問:予想される絶対誤差または予想される二乗誤差の意味でより優れている代替推定器はありますか?t(3)t(3)t(3)