1

R で data.frame コマンドを実行すると (以下に示す - 「マクロ」はモデル内で関心のある変数であることに注意してください)、変数 、fit、se、lowerおよびの出力が得られますupper。を除いて、各出力が私に何を伝えているかを知っていますfit。

> data.frame(effect(c("Macro"), model))
    Macro   fit        se         lower      upper
  1     C   45.30041   5.650558   34.14164   56.45918
  2     R   33.73317   4.394917   25.05406   42.41229

効果コマンドを実行すると (最初は標準偏差が得られると思っていました)、フィットと同じ数値が得られます。

> effect(c("macro"), model)
   Macro effect
   Macro
          C          R
   45.30041   33.73317

これは本当に標準偏差ですか、それとも平均をよりよく表しているのですか? そしてもちろん、私がこれらの潜在的な解釈の両方から完全に外れているという選択肢が常にあります.

4

1 に答える 1

2

fit回帰モデルが与えられた場合の当てはめ値または予測値を表します。カテゴリ予測子の場合、それは平均です：

library(effects) ## to access the effect() function

m1 <- lm(weight ~ group, data = PlantGrowth)

data.frame(effect(c("group"), m1))
  group   fit        se    lower    upper
1  ctrl 5.032 0.1971284 4.627526 5.436474
2  trt1 4.661 0.1971284 4.256526 5.065474
3  trt2 5.526 0.1971284 5.121526 5.930474

# CALCULATE MEANS
aggregate(weight ~ group, data = PlantGrowth, mean)
  group weight
1  ctrl  5.032
2  trt1  4.661
3  trt2  5.526

effect()それが標準偏差を与えると思った理由がわかりません。?effect関数を使用すると何が得られるかを確認してください。

于 2017-02-21T21:34:37.487 に答える