問題タブ [quantile]

scipy のいくつかの方法を使用して、一連の x (分位) および y (値) サンプルから分位関数(逆 CDF)を補間しようとしています。これは分位関数であるため、値が繰り返されることがあります。たとえば、CDF は最終的に 1 で平坦になるため、x=1 はいくつかの増加する y 値に対して繰り返されます。

一部の補間方法は関数用に構築されているようで、x の繰り返しには満足していません。x の繰り返しは導関数であると仮定する人さえいます (例: Krogh )。

どうすればこれを回避できますか？

私はmatplotlibを初めて使用し、次の情報を使用してプロットを作成したいと考えています:

1. 約 200 個の可変長ベクトルの中央値を結ぶ線 (入力)
2. これらのベクトルの対応する分位数を結ぶ線。
3. 対応するスプレッド (最大点と最小点) を結ぶ線。

つまり、基本的には連続箱ひげ図のようなものです。

ありがとう！

実際の値を含む列を持つデータ フレームがあります。

これらの値を相続人のサイズに応じて分類する追加の列が必要です。たとえば、値が最大値のグループの最小値のグループに属しているかどうかを知りたいです。これら 2 つのグループの要素数を同じにしたいと考えています。

例えば。次の値がある場合:

1、2、3 を 1 に、4、40、50 を 2 にマッピングしたいと思います。データ フレームで簡単に行う方法はありますか。

上記の例では、2 つのグループのみを使用しました。しかし、私はそれを柔軟に保ちたいと思います。たとえば、3 つのグループの場合、1 と 2 を 1 に、3 と 4 を 2 に、40 と 50 を 3 にマッピングしたいと思います。

経験的データを使用して次の積分を再作成しようとしています:

ここで、F、G は cdfs であり、それらの逆数は分位関数です。

これが私のコードです：

これで私はそれを期待しeqces(u,u) = 0ますu = np.random.uniform(0,1,50)が、これは一般的には当てはまりません。私が何か間違ったことをしているのか、それとも代替案を提案しているのか、誰にもわかりますか?

Edit

このコードは、いくつかの分析結果でうまく機能するようです。

2 つの患者タイプ (患者 1 の 33) と (患者 2 の 99) で構成される 132 x 107 のデータセットがあります。

外れ値を探しているので、次のコマンドを使用して、データセットで pca を実行し、最初の 4 つのコンポーネントの qqplot を実行しました

次を使用して2番目のコンポーネントのqqplotを実行すると:

次のグラフは、2 つの明確な外れ値を示しています。左下隅の赤い点は患者 1 です。

これらのポイントを削除できるように、データ内のこれらのポイントのインデックスを計算する簡単な方法はありますか?

パッケージを使用してquantreg、R で次の分位点回帰を実行しています。

を介して係数と要約統計量を抽出します。

私が得た結果は、

ただし、R^2/調整済み R^2 を知りたいです。これは、summary()-コマンドが単純な OLS 回帰に対して提供しているように見えますが、分位回帰の場合はそうではありません。

それらを抽出する方法を知っている人はいますか？