6

私は 2 方向反復測定計画 (3 x 2) を使用しており、効果サイズ (部分イータ 2 乗) の計算方法を知りたいと考えています。

その中にデータを含むマトリックスがあります（aと呼ばれます）のように（反復測定）

         A.a          A.b           B.a        B.b          C.a           C.b
1        514.0479     483.4246      541.1342   516.4149     595.5404      588.8000
2        569.0741     550.0809      569.7574   599.1509     621.4725      656.8136
3        738.2037     660.3058      812.2970   735.8543     767.0683      738.7920
4        627.1101     638.1338      641.2478   682.7028     694.3569      761.6241
5        599.3417     637.2846      599.4951   632.5684     626.4102      677.2634
6        655.1394     600.9598      729.3096   669.4189     728.8995      716.4605

idata =

    Caps    Lower
       A       a
       A       b
       B       a
       B       b
       C       a
       C       b

カーパッケージで反復測定ANOVAを実行する方法を知っています（タイプ3 SSは私の分野では標準ですが、論理エラーが発生することはわかっています..誰かが私が5歳のように私にそれを説明したいなら、私は大好きですそれを理解するために):

summary(Anova(lm(a ~ 1),
          idata=idata,type=3, 
          idesign=~Caps*Lower)),
    multivariate=FALSE)

私がやりたいことは、要約印刷のこの部分を取ることだと思います：

球状性を仮定した単変量 III 型反復測定 ANOVA

                     SS num Df Error SS den Df        F    Pr(>F)    
(Intercept)     14920141     1   153687      5 485.4072 3.577e-06 ***
Caps            33782        2     8770     10  19.2589  0.000372 ***
Lower           195          1    13887      5   0.0703  0.801451    
Caps:Lower      2481         2      907     10  13.6740  0.001376 **

そして、それを使用して、部分 ETA の二乗を計算します。したがって、私が間違っていなければ、最初の列から SS を取得し、それを各効果の (それ自体 + その行の SS エラー) で割る必要があります。これは正しい方法ですか？もしそうなら、どうすればいいですか？サマリー出力から値を参照する方法がわかりません。

4

1 に答える 1

8

部分イータ二乗は、パッケージ内のetasq関数で計算できますheplots

library(car)
mod <- Anova(lm(a ~ 1),
idata = idata,
type = 3,
idesign = ~Caps*Lower)

mod

library(heplots)
etasq(mod, anova = TRUE)

あなたは計算について尋ねているので:
From ?etasq: '単変量線形モデルの場合、古典的なη^2 = SSH / SSTおよび部分η^2 = SSH / (SSH + SSE). これらは一方向設計では同一です。

class例のように を持つモデルの η^2 の計算のコードを調べたい場合は、 を使用できますgetS3method(f = "etasq", class = "Anova.mlm")。

于 2013-08-06T23:45:52.277 に答える