問題タブ [standard-error]

過分散データを使用してポアソン回帰を実行しようとしているため、ハブホワイトの堅牢な標準誤差を使用する必要があると思います。

ただし、glmfit にはそのためのオプションがありません。そして、私が理解していることから、robustfitは線形回帰専用です。私はこれで正しいですか？Matlab で堅牢な標準誤差を使用してポアソン回帰を行う方法はありますか?

まず、この質問がすでに回答されている場合はお詫び申し上げます。これは、stackoverflow に関する私の最初の投稿です。

この作業例のように構成された株式市場データのパネルがあります。

基本的に、x1 を従属変数として、リグレッサーとして x3 を使用し、x2 のラグとリードを使用して時系列回帰を推定したいと考えています。Newey West の標準誤差を使用して、サンプルの各株に対してこれを実行したいと考えています。

ここでいくつかの投稿に行った後、私は次のことを思い付くことができました：

私は時系列で作業していて、ラグ演算子とリード演算子を使用したいので、zooorxts形式を使用する必要があることを読みました (使用したい関数とzoo互換性があるため、形式を選択しました):dyn$lm

zooパッケージにはすべての数値変数が必要なため、Stock 列を数値に変換することに注意してください。

ただし、これにより次のエラー メッセージが表示されます。

最終的には各株式の回帰を別々に推定したいので、これは問題ではないと考えました。これで問題は解決します。

このオブジェクトに対して、各株式の myLm の回帰を推定したいと考えています。

パッケージを使用しplyrて、関数を在庫ごとに分割しようとします。

残念ながら、これにより別のエラーが発生します。

ここで何が間違っていますか？

私はRに比較的慣れていないので、棒グラフにエラーバーを追加する方法を見つけようとしています. 簡単な例を使用すると、誤差範囲を追加したいと考えている 2 年間の細菌の有病率データがあります。まず、x 値と y 値、および 95% 信頼区間の標準誤差を含むデータ フレームを作成します。

次に、エラー バーの上限と下限を設定します。

次に、グラフを p として設定します。

次に、バーをグラフに追加します。

バーの幅を選択します。

次に、エラー バーを追加してみます。

エラーバーを追加すると、グラフがバーからラインに変わります。エラーバーは含まれていますが、データを適切に表すには棒グラフが必要です。どんな助けでも大歓迎です！

2 つのカテゴリ因子 (「生息地」と「地域」) と 1 つの連続変数 (T) があります。「生息地」には 2 つのレベルがあり、「場所」には 8 つのレベルがあります。デフォルトのウィスカーを SE を表すように変更し、中央値を各箱ひげ図の平均値に変更したいと考えています。これを行う方法と、プロット時に両方のカテゴリ要因を考慮に入れる方法はありますか? よろしくお願いします。

これは、boxplot ggplot のデフォルト設定で行ったことであり、中央値間隔で 1 番目と 3 番目の四分位数を示しています。

ここで非常によく似た質問が今まで回答されていないため、この質問を投稿しています。

私は、患者のコホート全体の平均 +/- SEM を、xyplot()すべての患者の値を表す上にプロットするように依頼されました。使用されるデータは、手術を受ける患者の術中の心血管所見を表しています。

これは私のdata.frame電話ですdf

xyplot変数用に最初に作成したものはCI次のようになります

次のようにして、コホート全体の平均 (黒い線) を追加することができました。

今、平均の上/下の線として平均に+/- SEを追加したいのですが、これを行う方法がどこにも見つかりません。

私ができることはlatticeExtra、以下のように黄土線 +/- SE を追加するパッケージを使用することですが、それは私が探している正しい数学関数ではありません。2 つの違いを説明するために、平均線をそこに残しました。

SO とインターネットで広範な検索を実行しましたが、これを行うための適切な機能を見つけることができませんでした。

助けていただければ幸いです。ありがとう。

R で、調査パッケージの svytotal 関数を使用して、ある郡から他のいくつかの郡に移動する人の数の推定値と標準誤差を生成し、反対方向への移動に同じ関数を使用すると、どうすればよいでしょうか。正味の移行と対応する標準誤差を計算しますか? IPUMSデータから一方向の移行を計算する方法は次のとおりです。

nowpumaは、人が現在住んでいる場所の名前を持つ要素です。leftは、その人物が以前の郡を離れたと識別された場合に 1 に等しくなります。結果は次のようになります。

私は次のようなことを望んでいます：

Stata から R へのロジット回帰を複製しようとしています。Stata では、オプション「ロバスト」を使用して、ロバストな標準誤差 (異分散一貫性のある標準誤差) を取得します。Stata からまったく同じ係数を複製できますが、パッケージ「サンドイッチ」では同じ堅牢な標準誤差を得ることができません。

OLS 線形回帰の例をいくつか試しました。R と Stata のサンドイッチ推定量は、OLS に対して同じ堅牢な標準誤差を与えてくれるようです。Stata が非線形回帰 (私の場合はロジット回帰) のサンドイッチ推定量を計算する方法を知っている人はいますか?

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

付属のコード: in R:

ステータス:

Linux で NetBeans 7.4 から stdout と stderr を取得するにはどうすればよいですか? コマンドライン (Ubuntu 14) から NetBeans を実行すると、bin/netbeans スクリプトで確認した内容によると、NetBeans はすべてを /dev/null にリダイレクトするように見えます。stdout と stderr が端末に送信されるのを確認したいと思います。出力をファイルにパイプするように bin/netbeans スクリプトを変更しようとさえしましたが、stdout と stderr のすべてが得られるわけではありません (ただし、私は bash スクリプトの専門家ではありません)。「netbeans --help」で利用できる貧弱なヘルプでは、stdout/stderr の手がかりが得られず、Google で何も見つけられませんでした。

Swing コンポーネントの「デザイン」ビューに切り替えると、NetBeans IDE 自体が例外をスローし、問題を追跡するために必要な完全なスタックトレースが得られないため、stdout と stderr が必要です (1 つの中にいくつかの静的な初期化コードがあります)。問題を引き起こしている可能性が高いクラスですが、実行時にうまく機能し、コードに問題は見られません)。Windowsでnetbeans.exeを実行すると、cmd.exeウィンドウにすべてのstderrとstdoutが表示されますが、Ubuntu 14でこれを行う必要があります。ありがとう。