COGTAG ERROil

EUCLID A : ' ¥ G L E CFCOXF OBLRES

:E'=Q △ ・ Q ・ '

ェ—'=Q △

--½-Q'

T,=L 、 + t J ; ‑ '

＾

Ao=AoTo

こ ^RETURN

F i g . 5 C o r b a l l i s & Traub

型の

F a c t o r m a x

型解法の概要フローチ

ヤート

Fig. 6

サブルーチン

CORBAL

のフローチャート

確認的因子分折における因子的不変性の評価（柴田・辻岡）

R

田

OUS8S

〕

ハウスホルダー法による固有値，固有ベクトルの算出。

⑥ (MOUT

〕

行列の出力。

READY SECTION における REDCAN,REDCA 2 , REDF A 2 , WRITE, WRIT 2 , WR‑

ITE3 の各サブルーチンは，各研究室独自の工夫を必要とするので，これらのプログラムリストは省略する。また，①,Rににおいて使用された KANDAIFACPAC のサプルーチン (VARI‑

MX, PROMAX, PLOTK) については，東村 ( 1 9 7 4 ) を参照いただきたい。

なお，辻岡・柴田 ( 1 9 8 3 ) および柴田・辻岡 ( 1 9 8 3 ) と重複するサプルーチンのプログラムリストは省略する。

‑125‑

横断的研究のための正準型解法のサプルーチン

(Pa tternmax型解法）

￨

SURRDUTINE CONFlECNV,NF,N、M,KP,ITTL,NMA,NMBJ DI MENS I ON AA (30,30), BB (30,30 l ,VA (30, 30) ,VB (30、30),VC(30,30J, VD C 30, 30), VE C:I0,30 J DIMENSION TC30,30l,TVC30、30),TF C 30, 30), TR <30,30)、TE(30,30J, Cf (30, 30), CR C 30, 30) ll IM卜NSI ON VV (60,60) ,Wli (60, 60) ,XX (60,60) ,yy (60,60) DIMENSION DC30l,EC30l,GC60l,Jl(60l IJIMfNSION JTTL(40l,NMAUO),NMU(、JOJ,NAM(60) C OAIA NR1,NR2,NR3/:J0,30,60/ C じ＊事軍•READY SHI ION じ CALL RfOCAN(NN1,Nll2、NV,N「,N,M,AA,U8,VF,Ol CALL wRIIL(Nfll,NR.l,NV,Nl ,AA,BIJ,VV,WW,(1,11,NMA、NMBl C

r.

寡竃事*~flUCNUSTES S[CT IUN C CALL Cl I FF CNR1, NR2, NV ,NF ,AA, BB ,VA, VB, T, TE, E ,NMA ,NMB) C C **事事LVAlUATION SECTION 2 (FACTOR LOADING MATRIX) じ CALL COGT AG CNR1, NR3, NV, NF ,VA, VB,VV ,WW, G, H) CA(L ERRORlNRl ,NV,NF ,VA,VB,VC、VD) CAIL EUCLIDlNRl,NV,NF,VA,VB,VC、VD、[) CALL ANGLUNRl, NV, NF, VA, VB, VC, VD ,VE, D, Fl CALL SOPERM lNR3, NRl, NV, NVV, NF, VV, VA, VB, NAM, NMA, NMB) CALL VARIMXlNH3,NR2,NVV、NF,VV, WW, 0, TV, NAM、0,1,0,ITTU CALL ARNGCG<NR3,NV,Nf,WW,XX,YY,G) CAIL PROMAX(NR3,NR2,NVV,NF,WW

、

VV,GF,CR

、

TF,「R,D, E, IT f L, NAM ,0, I, 0) CALL ARNGCGlNR3,NV,NF,WW,XX,VY,G) CALL ARNGCG(NR3,NV,NF,VV,XX,YY,Gl IF (KP , EQ, 1l ‑CALL Pl.Li「K(NH3, NR2, NVV ,NF ,VV,‑1, 0, 1 .'O,O, 0, !TTL, CR,CF ,NAM) REI URN I・ND

Cliff の解法のサプルーチン

SUIIRUU r r NE Cll FFINRl, NH2, !IV, NI ,AA,

。

U,Vll,Vll,I, IIこ,E,rJM,¥,NMIJ) OIMLNSION MCNll1,1l,BBINll1,ll,VAINRIヽll,VBCNR1,1l、 1 (NR2,1 l, IECNR2,1 l OlMINSJIJN E<ll,NMA(l),NMOII) じ WRllE(6,:SOOl 300 FIJRMAll!HI/SX,'PNOCHU~IES CI.IFF(lNlER) ME1110D'/l

じじ

富事富寡S = A'•R ・・‑‑‑‑> S'寡S‑‑‑‑‑> cl•ll•Q' じ tl(J 10 1=1,NF no 10 .1°1、N「 VA(l,J):0.0 DO 10 K=l、r,v 10 VA(l,J)cVA<f,J)1AA<K,ll•08(K,J) IHI Vi 1=1,NI IHI 15 Jこ1,NF l[(l,.f)clJ.O IJIJ 15 K• 1 ,NF 15 Tf(J,,fl=H(l,J)IVA(K,ll•VACK、J) C C,11.l IIOIISRSCNR2,NF,NF,Nr,1E,E,1) じ

（：

富＊事事p = I'凛IJ‑1/2寓Q'皐(lj寓IJ‑1/2) G IHl ,IJ J 01 ,NI IJII 2(1·.J•.J ,NF r IC I ,.fl ,・IJ. IJ Oil 20 K・.I, NF 20 HC[,.flcff.(i,-fHVA(l,Kl•TCK,.J) flfl 2'.> ,f,.1 ,NF 00 ?'.> I c 1 ,NI 2 5 I LC r , J l・. H. Cf ,,, l / Sil 111 Cr (.))) C C 事事軍事••\•IA ... , UlllHfJ<;IJNAI Sflf.llTIUN MATRIX C VA J C ***寓•T<•flf···-··> fllllflUfillNAL Sill !JI ION MATRIX (VB) じ Oil 1,0 I cJ ,NV flll 40 ,I・1,Nf VA (I,,!) 00.0 ilCI t,O Kこ1,NF 110 VA< [ ,,J)0VA< I ,Jl IAA(I ,Kl<TECK,J) 110 4; fcl,NV /JO 1>5 ,I・1 ,NF VBCl,JJ二o.o 110 1,5 Kこ1、NF 45 Vfl(l,.flcVfl(l,.f)l8H(J,Kl拿T(K,.fl C 咄Iff. (6‑:1221 322「llllMATClfl,5X,'IN「I●~ ・1 A'/l CALI Milli「CNR2、Nf,NF,IE,0,0,0,Ul WRI rru,,、,2,,) 324 FORMA'f<1110/5X,'INTEfl'Tfl'I) CAI.I MOlf I CNR?,NF ,Nf, T ,0,0,0,0) WRI IE<6‑:J12) 312 FORMA.IC1111//5X,'OIHHUGONAL SOLUTION MATRIX (VA)<I) CALL MOUTCNHl,NV,N「,VA,NMA,0,1,0) WR! TE Cfi,314) 314 FORMATCIHO//SX,'ORTHOGONAL SOLUTION MATRIX (VB)'/) CAI.I. MOUl CNR!, NV ,NF ,VB,NMB,0, 1 ,Ol RE TURN END

豆岡汁柿「芹砂禄蛍裕漉」瀕

^16~ffi

%

超行列を構成するためのサブルーチン

SUBROUI IN[ SllPLRM(NR3,NR1,NV,NW

、

NF,VV,VA,Vll

、

NAM,NMA,NMIJ) DIMENSION VVCNR3、!.),VACNR1 ,1 l ,VB CNR! ,1 l ,NAM Cl l ,NMAC!)、NMBCl) DD 10 l=!,NV I I=NV+ I DD 10 J=1,NF VV<I,J)=VACJ,J) VVCI IヽJ)=VB (IヽJ) 10 CONTINUE

DD 15 I二I,NV IIニNVfJ NAM(I J二NMA(IJ NAM (I IJニNMD(I J 15 CDNTJNUE C NVV•I• V•2 c̲ RE TURN 「ND.

縦断的研究のための正準型解法のサプルーチン

(Factormax型解法とPatternmax型解法）

超行列の 2 つの部分行列の一致性係数を求めるためのサプルーチン

￨ 127ーー

SUBROUTINE ARNGCGCNR3,NV、NF,WW,XX,YY、Gl C C

・軍"

COEFFICIENT OF CONGRUENCE AFTER ARANGE C DIMENSION WりCNR3ヽ1)、XX(NR3,1l、WCNR3,1l,G(1l C 00 50 1=1,NV I l=NVt I 00 50 J=l,NF JJこNF+J XXCl,J)=WW(J,J) 50 XX CI、JJl=WWCll,Jl N=NF•2 00 55 I =1,N 00 55 J=l、N YYCJ,Jl=O.O 00 55 K=l ,NV 55 YY(!,Jl=YYCl,JltXXCK,Jl•XXCK,J) 00 60 1=1,N 60 G(ll=l.0/SQHT<YY(l,l)l 00 65 J=l ,N 00 65 1=1,N 65 YYCl,Jl=G(ll•YYCl,Jl•G<J) WR I TE (6ヽ351) 351 FORHAT(1H1///5X,'FACTOR CONGRUENCE MATRIX CHOIAIIUNl'/l CALL MOUT<NRJ,N,N,YY,0,0,0,0l WHITEC6,352l 352「ORMAT < 1 H0// 5X,'0 I AGONAL or CONGRUENCE• /l 00 67 1=1,NF I l=NF II 67 WRITEC6,353) l,YYCII、I) 353 FORMATC1H ,5X,13、「8.5) C CT=O.O 00 30 1=1,NF I l=Nft I 30 cr‑c11YY(ll,I) f NV •• Nf CT二CT/「NV WRITE<6d55l. CT 355 FORMAT(lHO/SX、'AV=',F10.5) C RETURN ENO

SUIIIIIJII I I NEじUNIr1 INV,NI ,N,M, ISi 1, [!;[ ,,,1,, I I fl ,NMA,NMIJ) OIMI NSIIIN AA(JO,.lOJ ,IIU(JO.:lOl ,VAC.lO,、iOJ,VUC30dOJ、VCU0,30), VIHJ0,30) ,Vf C30,30J OIMINSIDN !(、<U.:lO),rv UO,JO), fl (,l0,30), TROO,JO), TE (30,30), ‑Cf UO,:<OJ,CHC.lOdOJ llJMI NSION RAUO,、!OJ,BB UO, JO>, RAU C.iO, 30l ,RVAC:l0,30J、!JVBCJ0,30), ‑WAO C:IO、.JOJ,WBO (30, .¥0 l ,WA (.HJ, 30) ,WU C:lO,、lO) DIMfN:SIIIN VV(60、f,OJ,WW (60,60), XX c r,0,60l ,YY Cf,0,60) DIM[NSlllN fl(30l,1.(、IOJ,Gl60J,11(60J OIMINSJON Jl!L(40J、NH~C30J,NMHC、lOJ,NAM(60) （； DAIA NNI ,NR2,NN.l/.l0,J0,60/ C C 竃寡X事lll'AIJY SECTION C CAI I. Ill llじA2(NR1 ,Nil? ,NV ,NF, N,M, IJA, llll, R,¥ll,AA, 1111,VE、0) ,; Al I will IE INRI、NIJ.¥,NV、NI,AA,llil,VV,wW,G,11,NMA,NMII) C,¥11 hRIIJ(NRl,NV、[IA,RII, R,¥H, NMA, NM Ill Tl―CISl.1 .l:Q. 1J JIIIN CALI, WG111H(NR1 ,NV,Nf ,AA,118,WAO、WIIO,VC, IJ, E, NMA, NMII J Fl SF ,； ALL WG1 ,,oR < NRl, NV, NF, HA, RB, AA, llll, WAil, wllO, IJVA, [!Vil、vc、VO,VE,O,E, NMA、NMIIJ ENIJ IF C （；事寡宣•l'RIJCHIJSIFS Sf.Cl IIIN C If (ISi 2 .EQ. ll IIILN CAIi. NISSEI CNR1,NIJ2,NV,NF,HAII、WAO,WHO,WA,WIJ ,E, T, TE、NMA、NMIJJ CALL IN IS IN <NRl ,Nil,, NV, NF ,M,1111, VA,Vll, H, J ,NMA,NMU) ELSE CAI.I Cl lrF<NR1,NIJ2,NV,NF,AA,BH,VA、VB,T, TE ,F ,NMA,NHB) CAI I. IN I w I N(NRl ,NRr, NV,Nf ,W,¥O,WIJO,WA,WB, r I., I ,NMA,NMB) END If C じ寡拿軍事IVAIIIAIIUN SI‑CllllN 1 C「ACfOIJSCORE MATHIXJ じ CAI I WARW02 CNR1 ,NV, NF, R,¥H ,WA ,W!l, VC ,VD, NMA、NMH) C （；富富寡宣fVAI IIAJION SECIION 2 (FACl!IR 1.IJAOING MATRIX) C CALL COGJAG<NR1,N¥J:l、NV,NF,VA,VBヽVV,WW,G,IIJ CAI I. 1‑llROR(NRl,NV,Nr,vA,VB,VC,VO) C,¥1.1. tllCLIOCNRl,NV,NI―, VA,VBヽVCヽVO,El CAI.I. ANGI.F CNH1、NV,Nr,VA,VB,VC、VO,VE,IJ,El CAI I. SIJPI.RMCNR3,NR1 ,NV,NVV,Nf ,VV,VA,VB,NAH,NMA、NMB) CAI I. VARJMXCNR3,NR2, NVV,NF ,VV,WW, IJ, lV,NAM,0, 1,0, I ITU CALL ARNGCGCNHJ,NV,Nf ,WW,XX,YY,G) CALL PJJOMAX CNR3, NIJ2, NVV, NF ,WW, VV, CF, CIJ, Tr, l lJ, IJ, E, !TTL, NAM,0, 1,0l CAI I. AIJNGCGCNJl3,NV,NF,WW,XX,YY,GJ CALL ARNGCGCNR3,NV,Nf,VV,XX,YY,Gl IF<KP .EQ. 1l ‑CALL l'LOIKCNR3,NR2、NVV、NF,VV,‑1,0,1.0,0.0, ITTl.,CR,CF、NAM)

RETIIIJN ENI)

吾楔忘四 +91き

^rj竺}が回

+H9

う焙寄

3

栞官︵涼日・汗耳︶

1281

Nesselroade の解法のサプルーチン

SUBNlllll I Nr: NE SSE LC NR1, NR2, NV, NF, llAB、WAO,WBO、WA,WB, E, T, TE, NMA,NMB) DIMENSION IJAIICNR1,1 l、WAO(NR1、1)、WBO(NR1, 1) ,WA (NR1,1) ,WB(NR1,1) IJIMfNSION

l<NR2,1l,!ECNR2,1l

DIMl,NSIUN E(1l,NMA(I)、NMIJC1l じ C 峯事＊事S = WAO'軍RA8•WBO‑‑‑‑‑> S'亭

s

‑‑‑‑‑> Q寓［）事Q' C Oil 10 I二1、N「 IHI 10 Jc1,NV 1,A (I ,,JJ co 0. 0 DO 10 K=1, NV 10 W,~<J,,l)=WACI、Jl1W,¥O <K, I l富IJAB(K,Jl IJO Vi 1‑1 ,Nf no 1~,1=1,Nr Wll CI ,J l一0.0 llll 1'> K=1,NV 1 o Wll C [ ,,IJ二りII(I,,[) 1 WA (I, Kl峯WllO(K,,I J ,; DO ;15 I二1,NF OIJ .i5 ,1°1,NF IE (I,、Il =O. 0 fJ[) ,J, K二1,Ne 35 TECl,Jl=TE<I、,JJ>WBCK,ll寡Wll(K,Jl （； CALL IIOIJS8S CNR2, NF, NF、Nf,IE,E,Tl （： WR!ll.(6,320) 320

FORMAT(1H1//SX,'Ull'

「1¥0(;1¥NALTRANSi llllMATIUN MATRIX (Q)'/) . CAI I. MOU I (NR2、N「,Nf,1,0,0,0,0l

G じ嵐拿＊峯P = P•0-1 /2事〇'•<Q•ll I /21 じ DO

58

I =1, NF DOり8,J、こ1,NI' TE (I ,Jl二0.0 )） 0 58 K=1,NF 58 TECJ,,ll=lf(I、,ll11;8(I ,Kl拿I(K ,Jl 00 60 ,J=1,NF DO 60 1=1、NF 60 TE CIヽJ)clECIヽJl/SQRT<E<Jll C WRI If (6,321 J ;121 I OIJMAl(1HO//SX,'ORTIIOGONAL TRANSFORMATION

MATRIX lP)'/) CALL

MDUT(NR2,NF,NF,TE

、0,0,0,0) C C 寡軍軍寡WA = WAO 寓IA <=Pl C 00 l¥O 1=1,NV DO 80 J=1,NF WA(l ,JJ=O.O on 80 K=1,NF 80 WA< IヽJ>=WA (J ,Jl♦WAO(l ,Kl•TE <K,J) C WRIIE(6,J30l 330

FORMAl(1H0//5X,'WAO

= ORTHOGONAL WEIGHT MATRIX

(SOLUTION)'/)

CALL MOUi CNR1 ,NV,NF

,WA,NMA,0,1,0)

C C n• 軍WR= WBO• TB C=Pl C OU 90 !=1、NV DO 90 Jo1,NF woe r ,JJ=o.o 00 90 K=1,NF 90 WBCl,J)=WBC!,JJIWBO(I

、Kl•TCKヽJ) WRITEC6,331l 331

FORMAT(1H0//5X,

、WBO= ORTHOGONAL WEIGHT MATRIX

CSOLUTJON)'/l

CALL MOUT <NR1, NV, NF、WB,NMBヽ0,1ヽ0) C

RETURN END

N esselroade の変換行列を用いて因子負荷量行列を回転

するサプルーチン

SUUROUIINE INTSTN<NR1、NR2,NV、NF,AA, OB ,VA,VU, TE, T ,NMA, NMU) DIMENSION AACNR1,ll、BBCNR1,I l ,VACNR1 ,! l ,VBCNRI ,!) , TE <NR2

、

ll,TC NR2, I l, NMA <I l

、

NMBCIl SlRIJCTURE ROTATION <AFT[R NESSELROAOE'S

M[THOlll

C C

富重凛富 C DO 10 !=I ,NV DO 10 J=I ,NF VA ClハJ),o. 0 VB<I ,J)=0,0 00 10 K•l,NF VA<l ,J)=VA<I ,J)tAAl I、Kl•H(K,J) 10 VB (IヽJ)=VB (I、J)tOBll,K)•f(K,Jl.

C WRI1E<6,300) 300

FORMAT<1Hl//5X,'AH

= A.• T <ORTJJOGONAL SOI.UTIIIN)'I) CALL MOIi「CNRI、NV,NF ,VA, NMA ,0, 1 ,O) WRllE(6d05) 305

FORMAT(1H0//5X,'OT

= 8• 1 <ORTIIOGONAL SOLU. 「ION)'I) CALL MOUT (NRI

、

NF,VB,NMU,0,1

ヽ0) C

RETURN ENll

Cli

釘の変換行列を用いて因子得点のための重み行列を回

転するサプルーチン

SUBROUTINE INTW「N<NR! ,NR2 ,NV ,NF

,W1¥0,WBO,l1A,WB,

IE, 1, NM,¥,NMB) 0 I卜l[NSlflNWi¥OCNR1,1)、¥jBO<NRI ,!) ,wA (NIii ,1) ,WIJ(Nfll, I)、 TE <NR2、1>,T<NR2、1),NMA (1), NMB C 1) WEIGIIT ROTATION <AFTER CLIFF'S METIIIJOl C

C••••

互耳汁柿「芹助柿器裕曲」瀕

_16~ffi

%

DD 10 1=1,NV DD 10 J=l ,NF

WA<I,J)=O.O WB(l,J)=O.O DO

10 K=l,NF WA (I ,J)=WA(! ,J)tWAO< I ,K)凜TE(K,J) 10

WBC!,J)=WB(l,JHWBO(l,K)

竃l(K,J) じ WRITE (6,300) 30.0

FDRMAl<lHl//SX,'WA

= WAO富TA(CLIFF)'/) CALL MOUT <NR1 ,NV ,NF

,WA,NMA,0,1,0)

C WRITEC6,305) 305.

FORMAl<1HD//5X,'WB

= WBD• TB (CLIFF>'ll CALL MOUT <NR1,NV,NF ,WB,NMB,0,1ヽ0) RETURN END

ー 129￨

縦断的研究のための Corballis 型解法のサプルーチン

（基準段定型） SIJBROUT I NE CONF2D <NV, NF ,N、M,KP,IT ,EP, ITll,NMA,NMBl DIMENSION AA(30,30)

、

BB(30, 30> ,VA (30, 3D), VB (30ヽ30l,VC(30,30l

、

‑VDC30,30l ,VEC30、30) DIMENSION T Cl 5,15), TV(! 5,15>, TF<I 5,15), TR(! 5,15l, TE(l 5,1 5)、 ‑TGC15,15l,THl15,15l,XL<15ヽ15) DIMENSION RAC30,30l

、

RB<30,30l

、

RAB(30, 30), RVA 130, 30> ,RVB 130

、

30), ‑WA(30,30l,WBC30,30l DIMENSION VVC60,60l,WW(60,60) OJ MENS JON 0(30) ,E(30l ,G(60) ,HC60l DI MENS I ON I TTL C40 l, NMA (30) ,NMB < 30), NAM C 60) DATA NRI、NR.2,NR3/30, 1 5, 60/ （；竃事軍富READY SEC I I UN CALL REDF A2 CNR1, NR2 ,NV ,NF ,N,M, RA, RB, RAB,AA, BB, TV, TF、TR,VE,Dl CALL XLM I NT<NR2、NFヽXU CALL WRITECNR!

、

NR3,NV,NF,AA

、

BB,VV,WW,G

、

H,NMA,NMB) CALL WRl13CNR1,NV,RA,RB、RAB,NMA,NMBl C n• 富PROCRUSTESSECT I ON CALL CORBAL INRI ,NR2,NV ,NF, IT ,EP,RAB,AA,BB,VA,VB,VC,VDヽXLヽT、TE, D, E, NMA) C 富竃事•EVALUATION SECTION I <FACTOR SCORE MATRIX)

・c

CALL WG140R<NR1 ,NV,NF、RA,RB,VA,VB,WA、WB,RVA,RVB,VC,VD,VE,D,E, NHA,NMBJ CALL WARWB2 <NRI, NV, NF, RAB,WA,WB,VC,VD,NMA, NMBl CALL WEIGTN(NR1,NR2,NV,NF,WA,WB

、 vc 、

TV,TF) CALL WARWB2 INRI, NV, NF ,RAB,WA,WB、VC,VO,NMA,NMB) C 冨竃寡畢EVALUATION .SECTION 2 <FACTOR LOADING MATRIX> CALL COGTAGCNRl,NR3,NV,NF、VA,VB,VV,WW,G,Hl CALL ERROR <NRI、NV,NF ,VA,VB,VC,VOl CALL EUCLJD<NR1,NV,NF,VA、VB,VC,VD,EJ CALL ANGLE<NRl,NV,NF,VA,VB、VC,VO,VE,D,El CALL CFCDNFCNR2,NF,T,TF,TH,TGl CALL OBLRES <NRI, NR2, NR3, NV, NF

、

VA,VB,VC,VO,TV

、

TF,TR

、

TE,TH,TG, ‑VV,WW,G,11,NMA、NMB,NAH,ITTL,3,KPl

RETURN END と

Traub の解法 Corballis のサプルーチン

CORBAL<NRI ,NR2,NV,NF、IT,EPヽRAB、AA,BB,VA,VB,VC,VD,XL,SUBROUTINE

、

NMA>T

、

D,E <NRI, I l, BB <NR1, I l ,VA <NR1, I l

、

VB<NRl,1l

、

D !MENS ION RAB (NRI ,I l ,AA VC<NRl,!l,VD<NR!,1l ‑ DI MENS ION XUNR2 ,I), T <NR2,1 l, TE<NR2,1 l DIMENSION D<1l,E(1l,NMA<ll

C C

nu A(A'Al‑1 & B<B'Bl‑1 C DO 5 1=1,NF DO 5 J=!,NF

T,E(l,J>=0.0 DO

5 K=l ,NV

5 lECl,J>=TECl,J)tAACK,J)軍AACK,J> CALL MATINVCNR2,NF,TE,0,El 00 10 I=! ,NV 00 10 J=1,NF

VACJ,J)=O.O 00

10 K=I ,NF 10 VACl,Jl=VACl,J)tAACl,K)富TECK,Jl じ DO 12 1=1,NF 00 12 J=l,NF TEC I ,Jl=O.O DO 12 K=l,NV 12 TECl,Jl=TECJ,J)t88CK,l>•BBCK、J) CALL MA1INVCNR2,NF、TE,0,E> 00 14 1=1,NV 00 14 J=l ,NF VB (I、J)=0,0 00 14 K=l、NF 14 VUCl,J>=VBCl,J)tBBCl,K)富TECK,J) C C ・重••ITERA110N SYCLE C L = CA'A>‑1 A'CRAB‑UAB)BCB'B)‑1 C UAR = OIAGCRAB‑ALB') C

LL=O WRITEC6,311)

311 FORMAH1Hl///5X,'ITER1110N EE') C 1 CONT I NUE LL=Lli 1 C C ‑‑‑‑‑> UAB ITERATION SYCLE C 00 15 1=1,NV 00 15 J=l,NF VCCJ,Jl=O.O 00 15 K=l,NF 15 VCC!,J)=VC(l,J>tAACJ,K)富XLCK,Jl C 00 20 1=1、NV DO 20 J=l、NV VDCl,J>=O,O 00 20 K=l,NF 20 VD(l,J)=VOCJ,J)tVCCJ,K)竃BBCJ,K) C DO 25 1=1,NV 25 OCl)=RABCI,ll‑VDCJヽI) C DO 27 1=1、NV 00 27 J=l ,NV VCCJ,Jl=RABCl,J) 27 IFCJ .EG. J) VCCl,J)=RABCl,J)‑O(J) じ WRITEC6,312l LL 312 FORMATC1H0///5X,'DIAGCRAB+ MATRIX)',5X,'LL =',13) WRITEC6,314l CVCCJ,J>,1=1,NV> 314 FORMATC1H0,5X、12F8.3) C C ‑‑‑‑‑> L‑MATRIX ITERATION SYCLE C 00 40 1=1,NF DO 40 J=l,NV VOCl,Jl=O.O

審挺菩四+9宝 ‑n 苫︱}か因+きう冷辛 8雫甫︵泄王・汗耳︶

￨ 130 1

00 40 K=l ,NV 1,0 VOCl,Jl=VDCIハJltVACK, I l•VC(K,Jl C OD 55 1=1 ,NF DO~5 J=t ,NF VCCIハJ>=O.0 00 55 K=l ,NV 55 VCCl,Jl=VCCIヽJltVDCI ,Kl•VB CK,Jl C C ‑‑‑‑‑>・SUM OF SQURE UK‑I) ‑UK) C DD 60 !=!,NF DO 60 J=t、NF 60 TECl,J>=XLCl,J)‑VCCI、J) r, DO 65 1=1,N「 DU 65 J=I、NF VDCIハJl=O.O DO 65 K=l、NF 65 VDCl,Jl=VDCI、JltTE CK, I>•TE CK,J > C

EE=O.O DO

70 1=1、NF 70 EF.=EEtVD CI, ll C WRITEC6,316l 316 FORMATC1H0/5X,'SUM DF SQURE <RESIDUAL)') WRITEC6,318l LL,EE 318 FORMATC1H0,5X,l5,2X、FI0.6) C DD 75 1=1,NF DO 75 J=t,NF 75 XLCl,Jl=VCCl,Jl

C C

‑‑‑‑‑> ITERATION RETURN OR END C IFCEE .LT. EP) GO TO 90 IFCLL .GE. Ill GO TO 90 C GD TO 1 C C

****

llERATION SYCLE END C 90 CDNT I NUE C C **軍•SOLUTION T = L CL'U‑1/2 C WRITEC6,334l LL 334 FORMATCIH1//5X,'L‑MATRIX CLL=',15,'l'/l CALL MOUTCNR2、NF,NF,XL,O、0,0,0) C DO 78 I=! ,NF DO 78 J=t ,NF TECl,Jl=O.O DO 78 K=t ,NF 78 TECl,Jl=TECIヽJltXL(K, ll•XL(K,Jl C CALL HOUS8SCNR2,NF、NF,NF,TE,D,Tl C WRITEC6,335) COCl>,1=1,NF)

335 FORMATC1 HOl/5X,'OIAGONAL ELEMENT'/5X,IOFIO. 5/l 00 80 I=! ,HF 80 0 Cl l =SQRTCO CI ll WRITEC&,33&) <D<ll、I=!,NF) 33& FORMATCIHOl/5X,'SQRTCOIAGl = ORTHOGONAL FC'/5Xヽ!OF10.5/l 00 81 !=!,NF 81 0(1)=!.0/0(1) C 00 82 1=1,NF 00 82 J=l ,HF TECl,Jl=O.O 00 82 K=l ,NF 82 TE Cl ,Jl=TE (l ,Jl♦T Cl ,Kl•O CKl•TCJ ,Kl C 00 84 1=1,NF DO 84 J=l ,NF TC! ヽJl=0‑0 DO 84.K=l,NF 84 TCl,Jl=TCl,Jl+XLCl,Kl•TECK,Jl C ‑ WR !TE (6,320) 320 FORMATC1 Hl //5X,'ORTHOGONAL TRANS FORMAT ION MATRIX CSOLUTIONl'/l CALL MOUT CNR2,NF ,NF, T ,0,0,0,0l PROCRUSTES ROTATION <FACTOR LOADING MATRIX) C C

****

C DO 86 I=! ,NV DO 86 J=I ,NF VACIハJl=O. 0 DO 86 K=I ,NF 86 VACl,Jl=VACl,Jl+AACl,Kl•TCK,Jl C DO 88 l=l,NY DD 88 J=I、NF 88 VB<I,Jl=BB<I,Jl C WRITEC6,325> 325 FDRMAT<1H0///5X,'PRDCRUSTES SDLUT(DN (VA)'/l CALL MDUT<NR1,NV,NF、VA,NMA,゜ヽ1,0) C REfURN END

涯淵汁柿「＃鼎柿器裕潮」滅

16~ffi

%

確認的因子分析における因子的不変性の評価（柴田・辻岡）

〔参考文献〕

浅野長一郎

1 9 7 1 因子分析法通論共立出版

B r o w n e , M. W. 1 9 6 7 On o b l i q u e p r o c r u s t e s r o t a t i o n . P s y c h o m e t r i k a , 3 2 , 1 2 5 ‑ 1 3 2 .

C a t t e l l , R . B . 1 9 4 4 " P a r a l l e l p r o p o r t i o n a l p r o f i l e s " and o t h e r p r i n c i p l e s f o r d e t e r m i n i n g t h e c h o i c e o f f a c t o r s by r o t a t i o n . P s y c h o m e t r i k a , 9 , 2 6 7 ‑ 2 8 3 .

C l i f f , N . 1 9 6 6 O r t h o g o n a l r o t a t i o n t o c o n g r u e n c e . P s y c h o m e t r i k a , 3 1 , 3 3 ‑ 4 2 .

C o r b a l l i s , M. C . & T r a u b , R . E . 1 9 7 0 L o n g i t u d i n a l f a c t o r a n a l y s i s , P s y c h o m e t r i k a , 3 5 , 7 9 ‑ 9 8 . E v a n s , G . T . 1 9 7 1 T r a n s f o r m a t i o n o f f a c t o r m a t r i c e s t o a c h i e v e c o n g r u e n c e , B r i t . ] . m a t h .

s t a t i s t . P s y c h o / . , 2 4 , 2 2 ‑ 4 8 .

G r e e n , B . F . 1 9 5 2 The o r t h o g o n a l a p p r o x i m a t i o n o f an o b l i q u e s t r u c t u r e i n f a c t o r a n a l y s i s . P s y c h o m e t r i k a ,

17,

4 2 9 ‑ 4 4 0 .

Harman, H . H . 1 9 6 0 Modern F o c t o r A n a l y s i s . U n i v . o f C h i c a g o P r e s s .

H e n r y s s o n , S . 1 9 5 7 A P P i i c a b i l i t y of f a c t o r a n a l y s i s i n t h e b e h a v i o r a l s c i e n c e s : A m e t h o d o l o g i c a l s t u d y , S t o c k h o l m : A l m q v i s t and W i k s e l l .

H o r s t , P . 1 9 6 1 G e n e r a l i z e d c a n o n i c a l c o r r e l a t i o n s and t h e i r a p p l i c a t i o n s t o e x p e r i m e n t a l d a t a . ] . c l i n . P s y c h o l , Monograph S u p p l e m e n t , N o . 1 4 , 3 3 1 ‑ 3 4 7 .

東村高良 1 9 7 4 多変量解析フ゜ログラムの使用法 ( 1 )KANDA! FAC‑PACK 関西大学電子計算機室，電子計算機利用研修会テキスト， N o .4 , 3 1 ‑ 4 3 .

J o r e s k o g , K . G . 1 9 6 9 A g e n e r a l a p p r o a c h t o c o n f i r m a t o r y maximum l i k e l i h o o d f a c t o r a n a l y s i s . P s y c h o m e t r i k a , 3 4 , 1 8 3 ‑ 2 0 2 .

M e r e d i t h , W. 1 9 6 4 N o t e s on t h e f a c t o r i a l i n v a r i a n c e . P s y c h o m e t r i k a , 3 4 , 1 8 3 ‑ 2 0 2 .

M o s i e r , C . I . 1 9 3 9 D e t t e r m i n 1 n g a s i m p l e s t r u c t u r e when l o a d i n g s f o r c e r t a i n t e s t s a r e known.

P s y c h o m e t r i k a 1 4 , 1 4 9 ‑ 1 6 2 .

N e s s e l r o a d e , J . R . 1 9 7 2 N o t e on t h e " L o n g i t u d i n a l f a c t o r a n a l y s i s " M o d e l . P s y c h o m e t r i k a , 3 7 , 1 8 7 ‑ 1 9 1 .

N e s s e l r o a d e , J . R . & B a l t e s , P . B . 1 9 7 0 Ona dilemma o f c o m p a r a t i v e f a c t o r a n a l y s i s : A s t u d y o f f a c t o r matching b a s e d on random d a t a . E d u c . p s y c h o l . m e a s m t . , 3 0 , 9 3 5 ‑ 9 4 8 .

N e s s e l r o a d e , J . R . , B a l t e s , P . B . & L a b o u v i e , E . W. 1 9 7 1 E v a l u a t i n g f a c t o r i n v a r i a n c e i n o b l i q u e s p a c e : B a s e l i n e d a t a g e n e r a t e d from random n u m b e r s . M u l t i v a r i a t e B e h a v i o r a l R e s e a r c h , 6 , 2 3 3 ‑ 2 4 1 .

P i n n e a u , S . R . & Newhouse, A . 1 9 6 4 Measures o f i n v a r i a n c e and c o m p a r a b i l i t y i n f a c t o r a n a l y s i s f o r f i x e d v a r i a b l e s . P s y c h o m e t r i k a , 2 9 , 2 7 1 ‑ 2 8 1 .

芝祐順 1 9 7 1 因子スコアの推定式高木貞二編現代心理学と数量化東京大学出版会芝祐順 1 9 7 9 因子分析法（第 2 版）東京大学出版会

柴田満 1 9 8 2 因子分析におけるプロクラステス諸解法の数値例比較（資料）関西大学大学院人間科学，

2 0 , 2 1 1 ‑ 2 2 5 .

柴田満因子分析におけるフ゜ロクラステス諸解法の位置づけ心理学評論（投稿予定）

柴田満・辻岡美延 1 9 8 3 確認的因子分析のための総合確認システム ‑Factormax 型解法と総合評

価システム—関西大学社会学部紀要， 15

( 1 ) , 1 4 5 ‑ 1 8 6 .

Schonemann, P . H . 1 9 6 6 A g e n e r a l i z e d s o l u t i o n o f t h e o r t h o g o n a l p r o c r u s t e s p r o b l e m . P s y c h o m e ‑ t r i k a , 3 1 , 1 ‑ 1 0 .

辻岡美延 1 9 7 5

確認的因子分析における因子構造の確認化ー~因子得点のための評価システム—•関西

大学社会学部紀要， 7( 1 ) , 9 8 ‑ 1 0 6 .

‑131‑

ドキュメント内確認的因子分析における因子的不変性の評価 (ページ 35-43)

EUCLID A : ' ¥ G L E CFCOXF OBLRES

:E'=Q △ ・ Q ・ '

--½-Q'

T,=L 、 + t J ; ‑ '

＾

Ao=AoTo

こ RETURN

F i g . 5 C o r b a l l i s & Traub

F a c t o r m a x

Fig. 6

CORBAL

R

OUS8S

ハウスホルダー法による固有値，固有ベクトルの算出。

⑥ (MOUT

行列の出力。

READY SECTION における REDCAN,REDCA 2 , REDF A 2 , WRITE, WRIT 2 , WR‑

ITE3 の各サブルーチンは，各研究室独自の工夫を必要とするので，これらのプログラムリスト は省略する。また，①,Rににおいて使用された KANDAIFACPAC のサプルーチン (VARI‑

MX, PROMAX, PLOTK) については，東村 ( 1 9 7 4 ) を参照いただきたい。

なお，辻岡・柴田 ( 1 9 8 3 ) および柴田・辻岡 ( 1 9 8 3 ) と重複するサプルーチンのプログラムリ ストは省略する。

‑125‑

横断的研究のための正準型解法のサプルーチン

r.

、

、

Cliff の解法のサプルーチン

。

（：

16~ffi

%

超行列を構成するためのサブルーチン

、

、

縦断的研究のための正準型解法のサプルーチン

超行列の 2 つの部分行列の一致性係数を求めるためのサプルーチン

・軍"

吾楔忘四 +91き

+H9

3

Nesselroade の解法のサプルーチン

s

58

N esselroade の変換行列を用いて因子負荷量行列を回転

するサプルーチン

、

、

、

、

Cli

釘の変換行列を用いて因子得点のための重み行列を回

転するサプルーチン

C••••

16~ffi

%

縦断的研究のための Corballis 型解法のサプルーチン

、

、

、

、

、

、

、

、

・c

、 vc 、

、

、

、

Traub の解法 Corballis のサプルーチン

、

、

、

、

、

****

****

16~ffi

%

確認的因子分析における因子的不変性の評価（柴田・辻岡）

こ ^RETURN

ITE3 の各サブルーチンは，各研究室独自の工夫を必要とするので，これらのプログラムリストは省略する。また，①,Rににおいて使用された KANDAIFACPAC のサプルーチン (VARI‑

なお，辻岡・柴田 ( 1 9 8 3 ) および柴田・辻岡 ( 1 9 8 3 ) と重複するサプルーチンのプログラムリストは省略する。

^16~ffi

_16~ffi

〔参考文献〕

東村高良 1 9 7 4 多変量解析フ゜ログラムの使用法 ( 1 )KANDA! FAC‑PACK 関西大学電子計算機室，電子計算機利用研修会テキスト， N o .4 , 3 1 ‑ 4 3 .

芝祐順 1 9 7 1 因子スコアの推定式高木貞二編現代心理学と数量化東京大学出版会芝祐順 1 9 7 9 因子分析法（第 2 版）東京大学出版会

柴田満 1 9 8 2 因子分析におけるプロクラステス諸解法の数値例比較（資料）関西大学大学院人間科学，

柴田満因子分析におけるフ゜ロクラステス諸解法の位置づけ心理学評論（投稿予定）

柴田満・辻岡美延 1 9 8 3 確認的因子分析のための総合確認システム ‑Factormax 型解法と総合評