問題タブ [robust]

私は ggplot2 を使用して、負の二項回帰の予測値をプロットしようとしています。1 つはバイナリ変数がオンになっており、もう 1 つはオフになっています。したがって、比較できる 2 つのプロットが 2 つあります。

ここのリンクは、ページの下部でそれを行う方法を示していますが、堅牢な標準誤差を使用して予測値のプロットの周りに陰影を作成できるようにしたいと考えています。predict() 関数からこれを取得する方法がわかりません。このコード例から、堅牢な標準エラーを取得して、プロットされた線の周りを陰にする回避策はありますか?

このサイトのコードを使用して、堅牢な標準エラーを生成します。

私が使用しているモデルはこれです：

そして、私が使用しているデータのサンプルは次のとおりです。

関数 rlm (MASS) は、ロバストな回帰のために M と MM の両方の推定を許可します。ggplot2 で MM ロバスト回帰からスムーザーをプロットしたいのですが、stat_smooth で method = "rlm" を選択すると、自動的に選択される推定方法が M タイプになると思います。

ggplot2 を介して rlm 関数の MM タイプ推定手法を選択する方法はありますか?

これが私のコードです：

rlm サマリー自体で結果を確認しましたが、ggplot2 が (デフォルト?) M 推定を使用していることは確かです。

rlm 関数から MM 推定を使用するにはどうすればよいですか?

よろしくお願いします！

私は現在、Net::FTPSSL::Robust で TLS 暗号化を使用して多数のサーバーのディレクトリをトラバースしています。ファイルをダウンロードしてから、ダウンロード中に削除したいと考えています。フォルダ名を削除したくありません。

問題は、ロバストには get() しかなく、必要なファイルを削除しないことです。これを行う方法が必要であることはわかっています。ここでわかるように、削除は Net::FTPSSL::Robust の一部ではありませんが、Net::FTPSSL の一部であるため、機能しません。

誰でもこれに光を当てることができますか、それともより良い解決策がありますか? Perl はまったく初めてなので、よろしくお願いします。

plm ライブラリを使用して固定効果回帰を実行し、sandwich,lmtest ライブラリを使用して堅牢な標準誤差を計算しています。回帰の実行に問題はありませんが、標準誤差を計算しようとすると、次のエラーが発生する場合があります。

係数または「通常の」標準誤差 (つまり、等分散) の計算に問題はありません。また、二次項を省略しても、ロバストな標準誤差を問題なく計算できます。

何が起こっているのか誰にも分かりますか？計画行列が特異な場合、係数は計算されていないはずなので、標準誤差を計算するときに問題がどこから来ているのかわかりません。

ありがとう！

キーのリストがあります['foo_a','foo_b','foo_c','fnord']

fnordここでの同様の解決策はすべて、テキストに 'sがないことを前提としています。

私は仕事をするこのコードを持っています:

これはもっとエレガントにできるのではないかと私は強く疑っていますが、私の python-fu は今のところ十分に強力ではありません。

いくつかの提案を聞きたいです。

編集します。追加の背景と説明。
これらは、他の開発者が翻訳に使用するためにアプリケーションに追加したキーです。それらには共通の接頭辞が必要ですが、人々は忘れてしまい、他のコードから切り取って貼り付けます。プレフィックスセパレータとしての「_」は単なる慣習です。セパレーターが使用されていると想定しないことをお勧めします。70% は完全に恣意的なしきい値です。「最も普及している」または「優勢」も機能します。
はい、これは python 2.7 であり、引用符内のスペースは単なる視覚的なアーティファクトです。

ロバストな回帰では、Matlab で最小絶対偏差回帰を行うために、どの重みを付ける必要がありますか?

堅牢な回帰用lmrobのライブラリを使用して、R で関数を使用しています。robustbaseとして使用しますrob_reg<-lmrob(y~0+.,dat,method="MM",control=a1)。私が使用する要約を返したいとき、summary(rob_reg)堅牢な回帰が行うことの1つは、データ内の外れ値を識別することです。要約出力の特定の部分から、次のことがわかります。

6508 observations c(49,55,58,77,104,105,106,107,128,134,147,153,...) are outliers with |weight| <= 1.4e-06 ( < 1.6e-06);

すべての外れ値をリストします。この場合は 6508 です (大部分を削除して ... に置き換えました)。これらの外れ値をどうにかして取得し、データから削除する必要があります。私が以前にしたことsummary(rob_reg)$rweightsは、観測のすべての重みを取得し、上記の例の特定の値よりも小さい重みを持つ観測を削除するため に使用することでし1.6e-06た。最初にすべての観測値の重みを取得せずに、外れ値のみのリストを取得する方法はありますか?