[論文レビュー] Developing methods for identifying the inflection point of a convex/concave curve

Demetris Christopoulos|arXiv (Cornell University)|Jun 24, 2012

Advanced Numerical Analysis Techniques被引用数 13

ひとこと要約

本稿では、回帰や曲率推定を用いずに、凸/凹曲線の拐点を特定する2つの新しい幾何的手法—極値面推定法（ESE）と極値距離推定法（EDE）—を提案する。代数的面積と左/右の弦からの距離を活用することで、曲線の拐点を含む区間を反復的に狭め、二分法に類似した収束を示す。ノイズのあるデータに対しても高い精度を発揮し、特にシグモイド曲線および3次多項式曲線に対して優れた性能を示す。

ABSTRACT

We are introducing two methods for revealing the true inflection point of data that contains or not error. The starting point is a set of geometrical properties that follow the existence of an inflection point p for a smooth function. These properties connect the concept of convexity/concavity before and after p respectively with three chords defined properly. Finally a set of experiments is presented for the class of sigmoid curves and for the third order polynomials.

研究の動機と目的

• 回帰やスプラインフィッティングに依存しない、曲線の拐点を特定するロバストで非パラメトリックな手法の開発。
• 左/右の弦からの距離と代数的面積に基づく幾何的基準を確立し、拐点を一意に特定する。
• 二分法に類似した収束行動を示す反復的アルゴリズムを構築し、真の拐点に収束させる。
• 3次多項式に対して、ESE推定値を解析的関係（拐点と弦の極値点の間）を用いて補正することで一貫した推定器を提供する。
• 測定誤差を伴う実世界のデータに対して実用的で計算効率の良いソリューションを提供し、シグモイド曲線および3次多項式モデルで妥当性を検証する。

提案手法

• 本手法は、曲線の端点と可変点xから定義される左弦および右弦を用い、これらの弦から曲線への代数的距離を計算する。
• 左および右の代数的面積積分、s_l(a,x)およびs_r(b,x)を定義し、その極値点x_lおよびx_rを用いて拐点の区間の境界とする。
• ESE法では、x_lとx_rの中点を推定器として用いる。EDE法では、対応する離散的推定器χ_lおよびχ_rの中点を用いる。
• 本アプローチは反復的である：各ステップで区間[x_l, x_r]が改善され、真の拐点pに収束する。
• 3次多項式の場合、x_l, x_r, a, bを組み合わせた補正式を導出することで、pの一貫した推定器が得られる。
• 曲率推定、回帰、平滑化を一切避ける。代わりに、弦と曲線の幾何的および積分的性質にのみ依存する。

実験結果

リサーチクエスチョン

• RQ1 回帰や曲率モデル化を用いずに、弦の幾何的性質と代数的面積に基づいて、凸/凹曲線の拐点を信頼性高く推定できるか？
• RQ2 左/右の弦から導出される面積積分の極値点を用いて、二分法に類似した収束行動を達成する方法は何か？
• RQ3 ESEおよびEDE手法の性能は、特にシグモイド曲線および3次多項式曲線に対して、独立同分布のノイズが存在する状況でいかがようか？
• RQ4 データの非対称性が生じる場合に、ESEの精度を向上させるための3次多項式に対する一貫した補正式を導出可能か？
• RQ5 曲率や平滑化に依存する従来の手法と比較して、本手法の計算効率および精度はどの程度か？

主な発見

• ESEおよびEDE手法は、二分法に類似した収束行動を示し、ノイズのあるデータに対しても高い精度で真の拐点に収束する。
• p=2.5、r=2.0の対称的3次多項式曲線に対して、補正済みESE推定器は2.4933を達成し、真のp=2.5に極めて近い。
• n=500、r=2.0の3次多項式曲線の場合、ESE法はχ_l=4.74およびχ_r=-0.26を生成し、補正推定値として2.4933が得られた。
• ESEおよびEDE手法は高速な収束を示し、標準のIntel Core i5 CPU上でn=10,000のデータ点を処理するFORTRAN実装が7秒未満で完了した。
• ESE法は、データの右側非対称性を示す対称的多項式ではバイアスを示すが、x_l, x_r, a, bを用いた導出された解析的補正式により是正された。
• Rパッケージ「inflection」は、ESEおよびEDE手法の実装を提供しており、FORTRAN, Maple, MATLABでも実装が可能である。