0

結果変数と、カテゴリ変数であるx3 つの説明変数があります。a, b, c私の例aでは、8 つのレベル、b4 つのレベル、c35 のレベルがありますが、3 つの変数のすべての組み合わせに観測値があるわけではありません (これはおそらく重要ではありません)。

Stata で次の加法的 ANOVA モデルを実行すると、

anova x a b c
adjust, by(a b) gen(y)

次に、変数およびxによって調整された変数の予測を取得します。調整コマンドは、結果ウィンドウに次の表を出力し、予測が調整された変数を生成します。aby

      |                 b                 
    a |       2        4        8       16
----------+-----------------------------------
   50 | .016655  .018487                  
   75 | .008286  .011237                  
  100 | .005937  .006677  .012467         
  150 | .001905  .004038  .009454         
  200 | .001774  .003107  .007592  .010081
  400 |          .004982  .006853  .009342
  800 |                   .002126   .00521
 1000 |                   .002732  .005221
----------------------------------------------
 Key:  Linear Prediction

私の問題は、変数yには と の各組み合わせの値があるのに対し、上の表にはa, bとの各組み合わせのc値しかないことです。テーブルから結果を保存するにはどうすればよいですか? 表の値と の値の関係は何ですか?aby

前もって感謝します。

更新：私はこれを見つけましたhelp adjust：

推定コマンドで使用されているが、by() 変数リストにも調整変数リストにも含まれていない変数は、観測ごとに現在の値のままになります。ここで、adjust は、by() オプションの変数によって定義された各グループ内のこれらの未指定の変数の平均を代用して、推定された予測の平均 (または対応する確率または指数化された予測) を表示します。

これは私のデータにも当てはまります。たとえば、a=75とのb=2場合c、値は 12、13、14、15、16 になります。y対応する値c=14(平均) は、まさに表に表示されているものです。しかし、値の平均がそれが取る値でない場合はどうなるでしょうか?

4

1 に答える 1

1

adjustこれは、とに関する Stefan Hansen のコメントへの返信marginsです。

一般的に、いいえ。ただし、すべてはモデルと、指定されたもの以外の共変量があるかどうかによって異なります。しかし、次の結果を考慮してください。

sysuse auto, clear 
anova mpg foreign rep78
adjust, by(foreign rep78)
margins foreign#rep78

ここで、結果は一致します。

私はmargins初歩的な用途を超えて流暢ではないので、より複雑な質問は他の人が処理する必要があります.

于 2013-01-16T12:05:17.580 に答える