私はaを実行しglm、各係数の標準誤差を抽出したいだけです。インターネットで関数を見ましたが、機能se.coef()せず、戻ります"Error: could not find function "se.coef""。

ggplot2 でモデルをプロットしたいと思います。信頼区間を推定するときに使用したい堅牢な分散共分散行列を推定しました。

ggplot2 に VCOV を使用するように指示することはできますか? あるいは、predict.lm に VCOV マトリックスを使用するように強制することはできますか? ダミーの例:

拡張された VCOV マトリックスを使用して「正しい」予測を取得できれば、ggplot に追加するのは比較的簡単です。

標準エラーをDOSバッチファイルのファイルにパイプするにはどうすればよいですか？標準出力のみを使用して配管する>>と、標準エラーは引き続きコンソールに送られます。

私の問題の詳細：

自動化された毎日のバックアップで、コマンドラインからWinRARを実行しています。そして、次の例では、WinRarの出力をパイプしますが、私が最も必要としているエラー出力はパイプしませんwinraroutput.txt。

問題は時々ファイルが使用されていることです、そしてそれらがアーカイブで見落とされたことを知りたいとき、そして私たちが戻る必要がある場合に備え.txtて各ファイルの隣のファイルにこれを記録し.rarます。不足しているファイルは、プログラムを再インストールすることで簡単に置き換えることができるため、不足していることがわかっている限り、それらを置き換えることは大したことではありません。ですから、必要なときに知っておくとよい情報であり、必要ではありません。

標準エラー出力だけを.txtファイルに出力し、可能であるが必須ではない場合は、通常の出力をコンソールに残すにはどうすればよいですか？

ボーナスポイント：

空白（エラーなし）の場合のファイルの削除方法を教えていただければボーナスポイント！ここで質問：Windowsバッチファイルを使用してファイルが空の場合、ファイルを検出（および削除）するにはどうすればよいですか？。

私は標準誤差のある単純な棒グラフをプロットしようとしていますが、それが私を夢中にさせています。私はいくつかの例を調べて、これまでに行きました：

私はスクリプトの最後の行を通り抜けようとしていますが、かなり長い間それを続けています：

私は本当にこれを乗り越えたいです！

ありがとう、

バズ

RでFama-MacBeth回帰を実行し、標準誤差を計算するパッケージがあるかどうか誰かが知っていますか？sandwichパッケージとそのNewey-West標準誤差を推定する機能、およびクラスタリング用の関数を提供していることを認識しています。しかし、ファーマ・マクベスに関しては何も見ていません。

私の問題はこれです：NAロバストな標準誤差の計算でいくつかの値を取得する必要がある場所を取得します。

クラスターロバストな標準誤差を使用して、固定効果パネル回帰を実行しようとしています。このために、私はp。3はStock/Watson（2006）に続く（後でEconometricaで公開され、アクセスできる人向け）。(M/(M-1)*(N-1)/(N-K) クラスターの数が有限でデータが不均衡であるため、下向きのバイアスに対抗して自由度を修正したいと思います。

同様の問題がStackOverflowの[1、2]の前に投稿され、 CrossValidatedの関連する問題[ 3 ]が投稿されています。

Arai（および最初のリンクの回答）は、関数に次のコードを使用します（以下にデータを追加してコメントを付けます）。

、ここでgcenter、平均からの偏差を計算します（固定効果）。次に、クラスター変数として回帰を続行しDS_CODEます（データに「データ」という名前を付けました）。

計算したい

ただし、分散のuj（clx上記の式を参照）を計算する場合、最初にリグレッサーのいくつかの値を取得し、次に多くのゼロを取得します。この入力ujが分散に使用される場合、NAs結果のみになります。

私のデータ

私のデータは特殊な構造である可能性があり、問題を理解できないため、Hotmailからのリンクとしてすべてを投稿します。その理由は、他のデータ（新井（2011）から取得）では問題が発生しないためです。ご迷惑をおかけして申し訳ございませんが、ご覧いただければ幸いです。このファイルは、純粋にデータを含む5MBの.txtファイルです。

からlmオブジェクトを取得し、標準エラーを抽出したい

関数の概要、名前、係数を知っています。

ただし、要約は、標準エラーに手動でアクセスする唯一の方法のようです。

どうすればseを出力できるかわかりますか？

対応するデータ行の横に実行されている小さな四分位数のプロットを含む、Tufteによって考案されたテーブルの素晴らしいバージョンがあることを覚えている場合：

パッケージNMOFと関数を使用してRにこのようなソリューションを実装しますqTable。これにより、基本的に上記のテーブルが作成され、LaTeXコードとして出力されます。

この視覚化の方法は、提示する情報が少量で、完全なグラフのためにスペースを無駄にしたくない場合に特に便利です。しかし、私は少しハックしたいと思いqTableます。分位数プロットを表示する代わりに、平均の標準誤差を表示したいと思います。私はそのような機能をハッキングするのは得意ではなく、ブルートフォースを使ってそれを行いました。qTable分位数を計算する関数の行を置き換えました。

標準誤差を計算する非常に残忍なもので：

私はそれがおそらく持続不可能なアプローチであることを知っていますが、実際にはある程度機能します-それは標準誤差をプロットしますが、プロットにこのギャップを保ちます：

そして、このギャップをなくしたいと思います。

これは、上記で説明した変更を加えたqTable関数です。

別のプログラムで実行された回帰分析の結果があり、それらが有意であるかどうかを R でテストしたいと考えています。ls.diag() が回帰結果の標準誤差と t 検定を計算することは知っていますが、非常に特殊な入力形式 (つまり、lsfit() の結果) が必要なので、それを使用できるとは思いません。r には、回帰分析の標準誤差と t 検定を計算する関数があり、関数に関連する係数を手動で簡単に与えることができますか?

新しい変数を作成するために使用するために、クラスター ロバスト標準誤差の値を保存する必要があります。

mean コマンドを使用してクラスターの堅牢な標準エラーを取得できますが、stata にはこの値が保存されません。

推定値のクラスター ロバスト標準誤差を計算し、それを使用して新しい変数を作成するためにこの値を保存する方法について何か提案はありますか?