問題タブ [confidence-interval]

PCA（主成分分析）プロットの信頼楕円の描画について、以前に質問に回答しました。

最終的に機能したのは、FactoMineRパッケージ（PCA、coord.ellipse、およびplot.PCA関数の組み合わせ）でした。これで信頼楕円を描くことができますが、coord.ellipse関数のlevel.confオプションが実際に何をするのかわかりません。信頼水準が高くなる小さな楕円を期待していましたが、逆のことが起こっています。

ご協力ありがとうございました！

私は次のモデルを持っています

とi を使用してpredict、level=0母集団の平均応答をプロットできます。母集団全体の nlme オブジェクトから 95% 信頼区間/予測バンドを抽出してプロットするにはどうすればよいですか?

Rでは、predict.lmは線形回帰の結果に基づいて予測を計算し、これらの予測の信頼区間を計算することもできます。マニュアルによると、これらの間隔はフィッティングの誤差分散に基づいていますが、係数の誤差間隔には基づいていません。

一方、ロジスティック回帰とポアソン回帰に基づいて予測を計算するpredict.glmには、信頼区間のオプションがありません。そして、ポアソンとロジスティック回帰に意味のある洞察を提供するために、そのような信頼区間をどのように計算できるかを想像するのにも苦労しています。

そのような予測に信頼区間を提供することが意味のある場合はありますか？それらはどのように解釈できますか？そして、これらの場合の仮定は何ですか？

最近のプロジェクトで modFit (R パッケージの一部、FME) を使用しています。モデル予測パラメータの信頼区間を取得することに関心がありますが、その方法がわかりません。modFit ルーチンの要約を見ると、各パラメーターの標準誤差と t 値が示されています。また、パラメータ相関表も提供されます。

この情報を使用してパラメータの信頼区間を抽出するにはどうすればよいですか? タスクを複雑にしすぎているように感じますが、行き詰まっています。どんな助けでも大歓迎です。

これは、私が取り組まなければならないこの要約情報です（質問に混乱を招くだけではないことを願っているので、必要に応じて無視してください）：

RパッケージrandomForestレポートは、フォレスト内の各ツリーの二乗誤差を意味します。ただし、データの各ケースの信頼度を測定する必要があります。randomForestは、単一ツリーの予測を平均することによってケースごとの予測を計算するため、ケースごとの標準誤差、したがって信頼区間も計算できるはずです。これは、出力のrandomForestオブジェクトを使用して実行できますか（ある場合：どのように？）、またはソースコードを掘り下げる必要がありますか？

「CI エリア」の下では、smth を意味します。@cbeleites の回答のグラフにプロットされているものや、Ron Pearson の投稿に示されているもののように.:

ありがとうございました。

正規分布を仮定して、信頼区間を計算したいサンプルデータがあります。

numpy および scipy パッケージを見つけてインストールし、numpy が平均値と標準偏差 (データがリストの numpy.mean(data)) を返すようにしました。サンプル信頼区間の取得に関するアドバイスは大歓迎です。

次のように定義されているいくつかのデータポイントをプロットする必要があります

c（x、y、stdev_x、stdev_y）

たとえば、ポイントとその周囲の1つの等高線を示す例として、95％の信頼限界を表す散布図として。理想的には、ポイントの周りの楕円形にプロットしたいのですが、その方法がわかりません。サンプルを作成してプロットし、stat_density2d（）を追加することを考えていましたが、等高線の数を1に制限する必要があり、その方法を理解できませんでした。

これは、Rを使用した入門統計の演習です。

rmrデータセットを使用して、代謝率と体重をプロットします。線形回帰モデルを関係に適合させます。適合モデルによると、体重70 kgの場合の予測代謝率はどれくらいですか？線の傾きに95％の信頼区間を与えます。

rmrデータ・セットは「ISwR」パッケージに含まれています。次のようになります。

与えられたxで予測されたyを計算する方法を知っていますが、勾配の信頼区間を計算するにはどうすればよいですか？

私はPythonが初めてで、lmfitパッケージを使用して自分の計算をチェックしようとしていますが、(1)次のテスト(および2)のエラーのデータ(sig)のエラーを含める方法がわかりません以下に示す conf_interval2d で取得します):

これは今までは問題なく動作していましたが、上記で質問したように、データ (sig_chla) の誤差を含めて、加重カイ 2 乗を計算できるようにするにはどうすればよいでしょうか?

パート 2: さらに、信頼区間 xs, ys, grid = conf_interval2d(mi, 'a', 'b', 20, 20) をプロットできるように、次を使用しようとすると

次のエラーが表示されます。

* ValueError: インテント (キャッシュ|非表示) | オプションの配列の作成に失敗しました -- 次元を定義する必要がありますが、(0,) を取得しました

また

ルーチン DGESV のパラメータ 4 が正しくありませんでした DGESV の Mac OS BLAS パラメータ エラー、パラメータ #0、(使用不可)、0 です