バイモーダルデータに対して対称の特定のカットを設定する

Question

私は両端でカットしたいバイナリアシメティック分布を持っています。その特定の部分は、私が各ベルの適切な側で対称境界を計算したいということです。この図は、単純化のためにベル間の分離の極端な場合を示しています。この場合

赤カットは、眼及び各サイドで使用1550本の青い線によって選択された潜在的にトリムするための機能を通過することができ、任意の値を表します。私の目標は、青い線の間のすべての部分集合です。

hist(p3_cut$x,50) 
abline(v=c(6200,7600),col='red') 
abline(v=c(6200-1500,7600+1500),col='blue') 

私の推測では、ここでの問題は基本的に各カーブの「エッジ」を見つけることです。私は手段の間に半分の距離を使用することはできません、私は0（または非常に低い値）から比較的高い何かに周波数変化を認識するものが必要です。

Answer 1

多少一般的な回答です。問題によっては、密度関数にビン幅を調整する必要があるかもしれません：

# get density of x and normalize so max is one 
dens <- density(x,adjust=0.1) 
dens$y <- dens$y/max(dens$y) 

# keep all x where density is higher than some fraction of max (here 1%) 
min_frac <- 0.01 
x_keep <- dens$x[dens$y > 0.01] 

# find position of gap in x, and get x just before and after gap 
gap_pos <- which.max(diff(x_keep)) 
left_cut <- x_keep[gap_pos] 
right_cut <- x_keep[gap_pos + 1] 

このコードを使用すると、私は、少なくともこの場合のほぼ完璧なカットを計算することができたdensity機能で調整パラメータを変更します。私は、このアプローチがこれに似ているほとんどの状況で十分柔軟であるということは肯定的です。私は提案されたカットの結果を示します。