「統計学関連なんでもあり」の過去ログ--- 042

No.10226　反復測定を多変量分散分析？　　【しち】　2009/06/29(Mon) 13:12

すみません，反復測定分散分析のデザインで多変量分散分析することについて，アドバイスをいただければと思います．

今，被験者間要因1（水準2）×被験者内要因1（水準3）という実験計画で1変量を集めたデータの解析を，SPSS16.0で試みております．
このとき，反復測定分散分析で解析したところ，「多変量検定」の結果が出力されました．
教科書を読んだところ，反復測定デザインの解析について被験者内要因の各水準の反応を変量とみなして多変量分散分析することもある，と書いてありましたので，その解析だと推測しました．
しかし，被験者間要因1（水準2）の実験計画で3変量を集めたデータとして，多変量分散分析してみると，出力される結果が違うことが分かりました．

こうした点について疑問を持ったのですが，反復測定分散分析で出力される「多変量検定」が，どのようなデザインの解析か，教えていただけないでしょうか．

No.10228　Re: 反復測定を多変量分散分析？　　【にゃんちゅう】　2009/06/30(Tue) 13:01

次の分析をしてみると同じ結果になりましたが，どう違うと言っているのでしょうか？


*********************************************.
title 森・吉田 1990 例題3・2・3 p107.
subtitle 1要因に対応がなく，1要因に対応がある場合.
data list free/sub a b1 b2 b3 b4.
begin data
 1 1 3 4 6 5
 2 1 3 3 6 7
 3 1 1 4 6 8
 4 1 3 5 4 7
 5 1 5 7 8 9
 6 2 3 2 3 2
 7 2 5 6 2 3
 8 2 2 3 3 3
 9 2 4 6 6 4
10 2 6 4 5 6
end data.
manova b1 b2 b3 b4 by a(1,2)
   /wsfactor=b(4)
   /design
.
GLM b1 b2 b3 b4 BY a
  /WSFACTOR=factor1 4 Polynomial 
  /METHOD=SSTYPE(3)
  /PRINT=DESCRIPTIVE 
  /CRITERIA=ALPHA(.05)
  /WSDESIGN=factor1 
  /DESIGN=a.
****************************************.

manova
* * * * * * * * * * * * * * * * * A n a l y s i s   o f   V a r i a n c e -- Design   1 * * * * * * * * * * * * * * * * *

 EFFECT .. a BY B
 Multivariate Tests of Significance (S = 1, M = 1/2, N = 2 )

 Test Name      Value        Exact F     Hypoth. DF     Error DF      Sig. of F
 Pillais         .74620       5.88019           3.00       6.00             .032
 Hotellings     2.94010       5.88019           3.00       6.00             .032
 Wilks           .25380       5.88019           3.00       6.00             .032
 Roys            .74620
 Note.. F statistics are exact.

 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

* * * * * * * * * * * * * * * * * A n a l y s i s   o f   V a r i a n c e -- Design   1 * * * * * * * * * * * * * * * * *

EFFECT .. B
 Multivariate Tests of Significance (S = 1, M = 1/2, N = 2 )

 Test Name    Value          Exact F  Hypoth. DF      Error DF        Sig. of F

 Pillais      .67370          4.12934         3.00        6.00             .066
 Hotellings  2.06467          4.12934         3.00        6.00             .066
 Wilks        .32630          4.12934         3.00        6.00             .066
 Roys         .67370
 Note.. F statistics are exact.

-------------------------------------------------------
GLM
			Multivariate Tests(b)
Effect		          Value	F	Hypothesis df	Error df Sig.
factor1	Pillai's Trace	         .674	4.129a	3.000	6.000	.066
	Wilks' Lambda	         .326	4.129a	3.000	6.000	.066
	Hotelling's Trace	2.065	4.129a	3.000	6.000	.066
	Roy's Largest Root	2.065	4.129a	3.000	6.000	.066
factor1 * a	Pillai's Trace 	.746	5.880a	3.000	6.000	.032
	Wilks' Lambda	　　　　 .254	5.880a	3.000	6.000	.032
	Hotelling's Trace	2.940	5.880a	3.000	6.000	.032
	Roy's Largest Root	2.940	5.880a	3.000	6.000	.032
a. Exact statistic
b. Design: Intercept + a 
 Within Subjects Design: factor1

あれ，MANOVAではRoy's Largest Rootの結果がちゃんとでていない。

No.10229　Re: 反復測定を多変量分散分析？　　【しち】　2009/06/30(Tue) 15:01

レスありがとうございます．

私はGUIのほうから多変量分散分析したのですが，試しにご提案のデータを解析しようとすると，以下のようになりました．

GLM b1 b2 b3 b4 BY a
/METHOD=SSTYPE(3)
/INTERCEPT=INCLUDE
/CRITERIA=ALPHA(.05)
/DESIGN= a.

これだと「a BY B」が含まれないからダメなのだということが分かりました．
ただやはり，変量でありながら因子になることがよく理解できないのですが…．

No.10230　Re: 反復測定を多変量分散分析？　　【にゃんちゅう】　2009/06/30(Tue) 15:06

WSFACTORおよびWSDESIGNの指定がないのが問題なのです。