問題タブ [confidence-interval]

ポアソン回帰 ( glm()with family='poisson') を使用し、次の形式で回答を提供したいと考えています。

「各予測値には、確率 0.95 の区間 [低、高] からの誤差があります」。

1. することは可能ですか？はいの場合、どのように？
2. いくつかの方法があります。

1) 残差を計算し、低い値と高い値を取得します。どのディストリビューションを使用すればよいですか?

2) 関数を使用しますpredict.glm(model, type = 'link', se.fit=TRUE)。正規分布の信頼区間を計算します。使用する

低い値と高い値を把握します。

2番目のケースを実装しました。問題は、上記の形式でどのように回答できるかです。

R で ccf を使用して、2 つの時系列データ セットの相互相関を計算したいと考えていました。私の質問は、プロット内の相関係数のいずれかが青い破線の外にある場合、手動で見ずにどうすれば知ることができるでしょうか? 何万もの時系列データを扱う必要があるためです。前もって感謝します！

Rの二項比率の区間推定の範囲確率を計算する一般的な関数を作成しようとしています。これは、さまざまな事前分布のWald、Clopper-Pearson、HPD区間などのさまざまな信頼区間法で行う予定です。

理想的には、間隔を計算するために使用する必要があるメソッドを引数として取ることができる1つの関数が必要です。それでは私の質問：関数を引数として別の関数に含めるにはどうすればよいですか？

例として、Exact Clopper-Pearson間隔の場合、次の関数があります。

ここで、Exact（x、n）は、適切な間隔を計算するための関数にすぎません。私はを頂きたい

そのため、間隔計算の方法ごとに個別のカバレッジ関数ではなく、カバレッジ確率を計算するための1つの関数があります。これを行うための標準的な方法はありますか？説明が見つかりませんでした。

ありがとう、ジェームズ

ターゲット画像で実際に顔認識を行うときの信頼レベルを収集する方法を見つけようとしています。いくつかの例を検索しましたが、実装方法を確認できるものは見つかりませんでした。すべての助けに感謝します、ありがとう。

極値分析（リターンレベルの計算）を実行するために使用している風データがあります。パッケージ「evd」、「extRemes」、「ismev」でRを使用しています。

ある期間Tのリターンレベル（RL）を推定するために、GEV、ガンベル、ワイブル分布を近似しています。GEV
とガンベルの場合、extRemes :: return.level（）を使用してRLと信頼区間を取得できます。働き。

いくつかのコード：

しかし、ワイブルの場合、CIを生成するためにいくつかの助けが必要です。

分位点回帰の信頼区間を構築する方法があります。

デフォルトは何ですか？5%ですか？10% など、別のものを見つける方法は?

関数への入力である多くのデータセットがあります。データはデータ テーブルに保存され、関数出力の信頼区間を計算しています。ただし、すべての入力データが同じで、次のようなエラーが発生する場合があります。すべての値が等しい場合の 0 または NA のような任意の値への間隔)? 例えば：

戻り値：

変化

戻り値:「t のすべての値は 100 に等しい \n 信頼区間を計算できません」 これは他のエラーにつながります。

結果を次のようにしたいと思います。

HH ライブラリの ci.plot を使用して、単純な線形回帰問題に対する信頼区間プロットを作成しました。

そして、それは本当にうまく機能します。ただし、プロットに追加したい。まず、時間内に簡単に順序を追えるように、ポイントに参加したいと思います。私は試した：

しかし、それは好きではありません:

...等

ci.plot に機能を追加する方法はありますか?

アップデート

私はちょうど実行してしまった：

図は次のようになります (破線は conf. int. ではなく予測間隔であることに注意してください。

私がまだできないことは、カラーマップの凡例を表示することです (Matlab では、これはコマンド: colorbar になります)。アイデアは、データのわずかなヒステリシスを強調することです。1995 年から 2048 年にかけて、色は青からピンクへと変化します。

カラーバーを表示する方法を知っている人はいますか?

一部のデータ

ここで要求されたのは、操作するデータのサブセットです。

タラバガニの成熟時のサイズを予測するために使用しているロジスティック回帰モデルがありますが、ブート パッケージを使用したブートストラップ用のコードの設定に問題があります。これは私が持っているものです:

カニの 50% が成熟するサイズに関心があるため、統計 CL=98.97 をブートストラップしたいと思いますが、関数を設定してその統計を指定する方法がわかりません。一般的なブートストラップ関数は言うまでもありません。私の 95% CI を取得するには、どんな助けでも大歓迎です! ありがとう！

ACGTE という文字を含むさまざまな長さの一連のn=400シーケンスがあります。たとえば、A の後にCがある確率は次のとおりです。

そして、これは一連の経験的シーケンスから計算できるため、

仮定:

次に、遷移行列を取得します。

しかし、ファットの信頼区間を計算することに興味があります。どうすればそれを実行できるかについて何か考えがありますか?