問題タブ [cdf]

まず、教授が間違った質問をしたかどうかはわかりません。とにかく、F(x)~U(0,1)CDFを生成しようとしましたF(x)=1-(1+x)exp(-x)(この CDF についてはx=g(F(x))、手計算では計算できませんでした)。そしてF(x)、質問が望むものを達成するために根を計算します。

0からまでINFのルート範囲uniroot()は論外だからです。したがって、ニュートン法を使用して記述します。

次に、私のコードは次のようになります。

ただし、F(x)が小さすぎてニュートン法の 1 ステップで結果が 0 未満になると、エラーが発生します。さらに、次のようにコードを修正しました。

ゼロに近い結果を拒否したため、明らかにコードが間違っています。

したがって、私の質問は、私のコードを修正して正しく計算できるかどうか、そうでない場合は、それを行う別の方法があるかどうかです。どんな援助でも大歓迎です。

累積分布関数 CDF を意味する x_i よりも高い値 X を取得する確率を取得したいと考えています。P(X>=x_i)。このコードを使用して、Matlab で実行しようとしました。

データが列ベクトル p1 にあると仮定しましょう。

今、CDFを計算したいのですが、

私が大丈夫なのか、それとも何か間違ったことをしているのか、誰かが私を助けることができるかどうか疑問に思っていますか?

取得した到着間隔の累積分布関数をプロットしたいのですが、データセットは次のとおりです。

他のトピックで見つけた多くの提案を試しましたが、うまくいきませんでした。

もう1つは次のとおりです。

cdf がどのように機能するか正確にわからないため、データに適用できるように、あなたの助けと説明の一部を使用できます..

ありがとう

タイムスタンプが秒単位の長いファイルを読み取り、numpy または scipy を使用して CDF をプロットする必要があります。numpy で試してみましたが、出力が本来あるべきものではないようです。以下のコード: 提案を歓迎します。

次のタイプの理論的意思決定モデルがあります。

\rhat = \r - \rbar = 平均からの偏差

\rhat<0 の場合: 1/(1+e^((|\rhat|^\delta)/t))

\rhat>0 の場合: 1-1/(1+e^((|\rhat|^\gamma)/t))

したがって、モデルは \gamma と \delta によってキャプチャされたほぼ 0 の非対称性を持つ S 字型の曲線です。曲線の傾きは t によって決まります。

\rhat のデータがあり、上記の S カーブをデータの CDF に当てはめて、\delta、\gamma、および t の推定値を取得したいと考えています。私は通常、この問題を二項選択問題と考えていますが、私のデータはそのタイプではありません。それは、分布パラメータを見つけることに沿ったものです。Mathematica の HistogramList コマンドなどを使用して、この推定の経験的 cdf データを回復する必要がありますか?