SET OLANG = English. set mprint off. DEFINE GAPIM_I (x = !charend('/') !default(X)/y = !charend('/')!default(Outcome) /groupid = !tokens(1) !default(group) /xn = !tokens(1) !default("'X'") /yn = !tokens(1) !default("'Outcome'")/groupn = !tokens(1) !default("'group'") /upper_lab=!tokens(1) !default("'Highs'")/lower_lab =!tokens(1) !default("'Lows'") /alpha=!tokens(1) !default(.05)/malpha=!tokens(1) !default(.10) /Xmiss=!tokens(1) !default(1)/ofile=!tokens(1) !default("'gapim_I.txt'") /directory=!tokens(1) !default("'c:\'") /clist = !cmdend). cd !directory. Comment to do list check to see if complete model runs X variable check make sure dich are effect coded recode continuous X variables more output for describing best model descriptives for X', I, and I'. Comment save original dataset. dataset name mydata. execute. dataset copy gapimdata. Comment Compute group composition variables. dataset activate gapimdata. aggregate outfile=agg.sav/break=!groupid /Xnn=n(!x) /Gidn=n(!groupid) /Xmean=mean(!x) /Xsd=sd(!x). Execute. get file="agg.sav". dataset name agg. Compute xmist=Gidn-Xnn. select if Xnn>2. Execute. dataset activate gapimdata. execute. MATCH FILES /table agg/file */ by !groupid. execute. Compute sy=1. if (sysmis(Xnn)) sy=0. if (!Xmiss =1 and Gidn>Xnn) sy=0. if (sysmis(!groupid)) sy=0. COMPUTE filter_$=(sy=1). VARIABLE LABEL filter_$ 'sy=1 (FILTER)'. VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'. FORMAT filter_$ (f1.0). FILTER BY filter_$. EXECUTE. Comment Compute the group variance and the actor and others similarity variables. SORT CASES BY !groupid. Compute mxmbxr = 1. if (!groupid=lag(!groupid))mxmbxr= lag(mxmbxr)+1. EXECUTE. compute Xvar=Xsd*Xsd. compute Xprime=(Xnn*Xmean-!x)/(Xnn - 1). compute I=!x*Xprime. If (xnn>2) Iprime=(Xnn*Xmean**2 - Xvar - 2*I)/(Xnn - 2). compute Xcon = !x - Xprime. compute XconI = I - Iprime. If (xnn>2) Xdiv = (2*I*(Xnn-1)+Iprime*(Xnn-1)*(Xnn-2))/(Xnn*(Xnn-1)). execute. OMS /SELECT TABLES /IF COMMANDS= ['Descriptives'] SUBTYPES=['Descriptive Statistics'] /DESTINATION FORMAT=SAV OUTFILE='file4.sav' /COLUMNS DIMNAMES=['Variables' 'Labels' 'Statistics']. dataset activate gapimdata. SET Printback=On. DESCRIPTIVES VARIABLES=!x !y /MISSING LISTWISE /STATISTICS=MEAN STDDEV MIN MAX. SET Printback=Off. OMSEND. get file="file4.sav". dataset name file4. flip. execute. dataset close file4. erase file="file4.sav". dataset name filenox. IF ($casenum=14) xmin=lag(var001,9). IF ($casenum=14) xmax=lag(var001,8). IF ($casenum=14) xmean=lag(var001,7). IF ($casenum=14) xsd=lag(var001,6). IF ($casenum=14) ymin=lag(var001,4). IF ($casenum=14) ymax=lag(var001,3). IF ($casenum=14) ymean=lag(var001,2). IF ($casenum=14) ysd=lag(var001,1). IF ($casenum=14) nxn=var001. Compute id=1. execute. compute ncn=$casenum. execute. select if NCN=14. MATCH FILES /FILE * /DROP ncn var001 case_lbl. execute. OMS /SELECT Tables /IF COMMANDS=['Frequencies'] Subtypes= [ 'Frequencies'] /DESTINATION FORMAT=SAV OUTFILE='freqn.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. dataset activate agg. SET Printback=On. FREQUENCIES VARIABLES=xnn /ORDER=ANALYSIS. SET Printback=off. OMSEND. Get file='freqn.sav'. Dataset name freqn. if (var2 ='Total') var2= '0'. Compute npg=number(var2,f10.0). If (npg=0)npg=$sysmis. Execute. OMS /SELECT TABLES /IF COMMANDS= ['Descriptives'] SUBTYPES=['Descriptive Statistics'] /DESTINATION FORMAT=SAV OUTFILE='file5.sav' /COLUMNS DIMNAMES=['Variables' 'Labels' 'Statistics']. SET Printback=On. DESCRIPTIVES VARIABLES=frequency npg /MISSING LISTWISE /STATISTICS= MIN MAX. OMSEND. get file="file5.sav". dataset name file5. Compute id=1. execute. MATCH FILES /FILE file5 /DROP command_ subtype_ label_ frequency_N npg_N. execute. OMS /SELECT TABLES /IF COMMANDS= ['Descriptives'] SUBTYPES=['Descriptive Statistics'] /DESTINATION FORMAT=SAV OUTFILE='filexm.sav' /COLUMNS DIMNAMES=['Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate agg. Compute xmist=0. If (Gidn>Xnn)xmist=1. DESCRIPTIVES VARIABLES=xmist /MISSING LISTWISE /STATISTICS= SUM. OMSEND. get file="filexm.sav". dataset name filexm. execute. Compute id=1. MATCH FILES /FILE * /DROP command_ subtype_ label_ ValidNListwise_N. execute. Comment Xmist_Sum and Xmist_N. OMS /SELECT Tables /IF COMMANDS=['Frequencies'] Subtypes= [ 'Frequencies'] /DESTINATION FORMAT=SAV OUTFILE='freqx.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. dataset activate agg. SET Printback=On. FREQUENCIES VARIABLES=Xmean /ORDER=ANALYSIS. SET Printback=off. OMSEND. Get file='freqx.sav'. Dataset name freqx. If (var2="Total") var2="5555". Compute nhom=number(var2,f10.3). Select if nhom=1 or nhom=-1 or nhom=5555. Compute id= 1. MATCH FILES /FILE * /DROP command_ subtype_ label_ var1 var2 percent validpercent cumulativepercent. execute. CASESTOVARS /id=id. Compute nhoml=0. Compute nhomu=0. If (nhom.1=-1) nhoml=frequency.1. If (nhom.2=-1) nhoml=frequency.2. If (nhom.1=1) nhomu=frequency.1. If (nhom.2=1) nhomu=frequency.2. dataset close agg. erase file="agg.sav". execute. Comment Get Total sample size. OMS /SELECT TABLES /IF COMMANDS= ['Descriptives'] SUBTYPES=['Descriptive Statistics'] /DESTINATION FORMAT=SAV OUTFILE='filet.sav' /COLUMNS DIMNAMES=['Variables' 'Labels' 'Statistics']. dataset activate gapimdata. SET Printback=On. DESCRIPTIVES VARIABLES=Xnn /MISSING LISTWISE /STATISTICS=MEAN. SET Printback=Off. OMSEND. get file="filet.sav". dataset name filet. flip. execute. dataset name filen. IF ($casenum=6) nxnt=var001. Compute id=1. compute ncn=$casenum. execute. select if NCN=6. MATCH FILES /FILE * /DROP ncn var001 case_lbl. execute. Comment Get X sample size. OMS /SELECT TABLES /IF COMMANDS= ['Descriptives'] SUBTYPES=['Descriptive Statistics'] /DESTINATION FORMAT=SAV OUTFILE='filex.sav' /COLUMNS DIMNAMES=['Variables' 'Labels' 'Statistics']. dataset activate gapimdata. SET Printback=On. DESCRIPTIVES VARIABLES=!x /MISSING LISTWISE /STATISTICS=MEAN. SET Printback=Off. OMSEND. get file="filex.sav". dataset name filex. flip. execute. dataset name fileq. IF ($casenum=6) nxnx=var001. Compute id=1. compute ncn=$casenum. execute. select if NCN=6. MATCH FILES /FILE * /DROP ncn var001 case_lbl. execute. Comment Get Y sample size. OMS /SELECT TABLES /IF COMMANDS= ['Descriptives'] SUBTYPES=['Descriptive Statistics'] /DESTINATION FORMAT=SAV OUTFILE='filey.sav' /COLUMNS DIMNAMES=['Variables' 'Labels' 'Statistics']. dataset activate gapimdata. SET Printback=On. DESCRIPTIVES VARIABLES=!y /MISSING LISTWISE /STATISTICS=MEAN. SET Printback=Off. OMSEND. get file="filey.sav". dataset name filey. flip. execute. dataset name filek. IF ($casenum=6) nxny=var001. Compute id=1. compute ncn=$casenum. execute. select if NCN=6. MATCH FILES /FILE * /DROP ncn var001 case_lbl. execute. Comment Empty Model Model 0 REML. dataset activate gapimdata. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var0.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. MIXED !y WITH !x Xprime I Iprime !clist /FIXED = !clist /METHOD = reml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var0.sav'. dataset name var0. If ($casenum=2) emptyevar=lag(estimate,1). If ($casenum=2) icc0=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. Comment Empty Model Model 0 ML. dataset activate gapimdata. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdevx.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. MIXED !y WITH !x Xprime I Iprime !clist /FIXED = !clist /METHOD = ml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. OMSEND. SET Printback=on. get file='fdevx.sav'. dataset name fdevx. compute emptydev=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . execute. Comment Main Effects Model Model 1 REML. dataset activate gapimdata. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var1.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe1.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Model Dimension'] /DESTINATION FORMAT=SAV OUTFILE='gcount.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. MIXED !y WITH !x Xprime I Iprime !clist /FIXED = !x Xprime !clist /METHOD = reml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='gcount.sav'. dataset name gcount. compute id=1. Compute ngroups=-1. if (var3='Number of Subjects')ngroups=Number(var4,F10.0). select if ngroups>0. MATCH FILES /FILE * /DROP var1 var2 var3 var4 command_ subtype_ label_ . execute. get file='var1.sav'. dataset name var1. If ($casenum=2) evar1=lag(estimate,1). If ($casenum=2) icc1=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe1.sav". dataset name pe1. If ($casenum=3) xpr1se=Std.Error. If ($casenum=3) xpr1p=sig. If ($casenum=3) xpr1=estimate. If ($casenum=3) x1se=lag(Std.Error,1). If ($casenum=3) x1p=lag(sig,1). If ($casenum=3) x1=lag(estimate,1). If ($casenum=3) int1se=lag(Std.Error,2). If ($casenum=3) int1p=lag(sig,2). If ($casenum=3) int1=lag(estimate,2). execute. compute ncn=$casenum. EXECUTE. select if NCN=3. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Main Effects Model Model 1 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev1.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime !clist /FIXED = !x Xprime !clist /METHOD = ml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev1.sav'. dataset name fdev1. compute dev1=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . execute. Comment Actor Only Model 2 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var2.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe2.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean !clist /FIXED = !x !clist /METHOD = reml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var2.sav'. dataset name var2. If ($casenum=2) evar2=lag(estimate,1). If ($casenum=2) icc2=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe2.sav". dataset name pe2. If ($casenum=2) x2se=Std.Error. If ($casenum=2) x2p=sig. If ($casenum=2) x2=estimate. If ($casenum=2) int2se=lag(Std.Error,1). If ($casenum=2) int2p=lag(sig,1). If ($casenum=2) int2=lag(estimate,1). execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Actor Only Model 2 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev2.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. dataset activate gapimdata. SET Printback=On. MIXED !y WITH !x Xprime I Iprime Xmean !clist /FIXED = !x !clist /METHOD = ml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev2.sav'. dataset name fdev2. compute dev2=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Others Only Model 3 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var3.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe3.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean !clist /FIXED = Xprime !clist /METHOD = reml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var3.sav'. dataset name var3. If ($casenum=2) evar3=lag(estimate,1). If ($casenum=2) icc3=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe3.sav". dataset name pe3. If ($casenum=2) xpr3se=Std.Error. If ($casenum=2) xpr3p=sig. If ($casenum=2) xpr3=estimate. If ($casenum=2) int3se=lag(Std.Error,1). If ($casenum=2) int3p=lag(sig,1). If ($casenum=2) int3=lag(estimate,1). execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Others Only Model 3 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev3.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean !clist /FIXED = Xprime !clist /METHOD = ml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev3.sav'. dataset name fdev3. compute dev3=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Group Model 4 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var4.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe4.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean !clist /FIXED = Xmean !clist /METHOD = reml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var4.sav'. dataset name var4. If ($casenum=2) evar4=lag(estimate,1). If ($casenum=2) icc4=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe4.sav". dataset name pe4. If ($casenum=2) x4se=Std.Error. If ($casenum=2) x4p=sig. If ($casenum=2) x4=estimate. If ($casenum=2) int4se=lag(Std.Error,1). If ($casenum=2) int4p=lag(sig,1). If ($casenum=2) int4=lag(estimate,1). execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Group Model 4 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev4.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean !clist /FIXED = Xmean !clist /METHOD = ml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev4.sav'. dataset name fdev4. compute dev4=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Contrast Model 5 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var5.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe5.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean Xcon !clist /FIXED = Xcon !clist /METHOD = reml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var5.sav'. dataset name var5. If ($casenum=2) evar5=lag(estimate,1). If ($casenum=2) icc5=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe5.sav". dataset name pe5. If ($casenum=2) x5se=Std.Error. If ($casenum=2) x5p=sig. If ($casenum=2) x5=estimate. If ($casenum=2) int5se=lag(Std.Error,1). If ($casenum=2) int5p=lag(sig,1). If ($casenum=2) int5=lag(estimate,1). execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Contrast Model 5 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev5.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xmean Xcon !clist /FIXED = Xcon !clist /METHOD = ml /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev5.sav'. dataset name fdev5. compute dev5=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Complete Model Model 6 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var6.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe6.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime !clist /METHOD = reml /FIXED = !x Xprime I Iprime !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var6.sav'. dataset name var6. If ($casenum=2) evar6=lag(estimate,1). If ($casenum=2) icc6=estimate. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe6.sav". dataset name pe6. If ($casenum=5) ipr6se=Std.Error. If ($casenum=5) ipr6p=sig. If ($casenum=5) ipr6=estimate. If ($casenum=5) i6se=lag(Std.Error,1). If ($casenum=5) i6p=lag(sig,1). If ($casenum=5) i6=lag(estimate,1). If ($casenum=5) xpr6se=lag(Std.Error,2). If ($casenum=5) xpr6p=lag(sig,2). If ($casenum=5) xpr6=lag(estimate,2). If ($casenum=5) x6se=lag(Std.Error,3). If ($casenum=5) x6p=lag(sig,3). If ($casenum=5) x6=lag(estimate,3). If ($casenum=5) int6se=lag(Std.Error,4). If ($casenum=5) int6p=lag(sig,4). If ($casenum=5) int6=lag(estimate,4). execute. compute ncn=$casenum. EXECUTE. select if NCN=5. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Complete Model Model 6 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev6.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime !clist /METHOD = ml /FIXED = !x Xprime I Iprime !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev6.sav'. dataset name fdev6. compute dev6=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Person Fit Model Model 7 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var7.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe7.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime !clist /METHOD = reml /FIXED = !x Xprime I !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var7.sav'. dataset name var7. If ($casenum=2) evar7=lag(estimate,1). If ($casenum=2) icc7=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe7.sav". dataset name pe7. If ($casenum=4) i7se=Std.Error. If ($casenum=4) i7p=sig. If ($casenum=4) i7=estimate. If ($casenum=4) xpr7se=lag(Std.Error,1). If ($casenum=4) xpr7p=lag(sig,1). If ($casenum=4) xpr7=lag(estimate,1). If ($casenum=4) x7se=lag(Std.Error,2). If ($casenum=4) x7p=lag(sig,2). If ($casenum=4) x7=lag(estimate,2). If ($casenum=4) int7se=lag(Std.Error,3). If ($casenum=4) int7p=lag(sig,3). If ($casenum=4) int7=lag(estimate,3). execute. compute ncn=$casenum. EXECUTE. select if NCN=4. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Person Fit Model Model 7 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev7.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime !clist /METHOD = ml /FIXED = !x Xprime I !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev7.sav'. dataset name fdev7. compute dev7=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Diversity Model 8 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var8.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe8.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xdiv !clist /METHOD = reml /FIXED = !x Xprime Xdiv !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var8.sav'. dataset name var8. If ($casenum=2) evar8=lag(estimate,1). If ($casenum=2) icc8=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe8.sav". dataset name pe8. If ($casenum=4) i8se=Std.Error. If ($casenum=4) i8p=sig. If ($casenum=4) i8=estimate. If ($casenum=4) xpr8se=lag(Std.Error,1). If ($casenum=4) xpr8p=lag(sig,1). If ($casenum=4) xpr8=lag(estimate,1). If ($casenum=4) x8se=lag(Std.Error,2). If ($casenum=4) x8p=lag(sig,2). If ($casenum=4) x8=lag(estimate,2). If ($casenum=4) int8se=lag(Std.Error,3). If ($casenum=4) int8p=lag(sig,3). If ($casenum=4) int8=lag(estimate,3). execute. compute ncn=$casenum. EXECUTE. select if NCN=4. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Diversity Model 8 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev8.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime Xdiv !clist /METHOD = ml /FIXED = !x Xprime Xdiv !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev8.sav'. dataset name fdev8. compute dev8=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Interaction Contrast Model Model 9 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var9.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe9.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI !clist /METHOD = reml /FIXED = !x Xprime XconI !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var9.sav'. dataset name var9. If ($casenum=2) evar9=lag(estimate,1). If ($casenum=2) icc9=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe9.sav". dataset name pe9. If ($casenum=4) i9se=Std.Error. If ($casenum=4) i9p=sig. If ($casenum=4) i9=estimate. If ($casenum=4) xpr9se=lag(Std.Error,1). If ($casenum=4) xpr9p=lag(sig,1). If ($casenum=4) xpr9=lag(estimate,1). If ($casenum=4) x9se=lag(Std.Error,2). If ($casenum=4) x9p=lag(sig,2). If ($casenum=4) x9=lag(estimate,2). If ($casenum=4) int9se=lag(Std.Error,3). If ($casenum=4) int9p=lag(sig,3). If ($casenum=4) int9=lag(estimate,3). execute. compute ncn=$casenum. EXECUTE. select if NCN=4. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Interaction Contrast Model Model 9 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev9.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI !clist /METHOD = ml /FIXED = !x Xprime XconI !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev9.sav'. dataset name fdev9. compute dev9=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Actor Model Model 10 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var10.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe10.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI !clist /METHOD = reml /FIXED = !x I !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var10.sav'. dataset name var10. If ($casenum=2) evar10=lag(estimate,1). If ($casenum=2) icc10=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe10.sav". dataset name pe10. If ($casenum=3) i10se=Std.Error. If ($casenum=3) i10p=sig. If ($casenum=3) i10=estimate. If ($casenum=3) x10se=lag(Std.Error,1). If ($casenum=3) x10p=lag(sig,1). If ($casenum=3) x10=lag(estimate,1). If ($casenum=3) int10se=lag(Std.Error,2). If ($casenum=3) int10p=lag(sig,2). If ($casenum=3) int10=lag(estimate,2). execute. compute ncn=$casenum. EXECUTE. select if NCN=3. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Actor Model Model 10 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev10.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI !clist /METHOD = ml /FIXED = !x I !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev10.sav'. dataset name fdev10. compute dev10=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Others Model, Main Effects and Interaction Model 11 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var11.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe11.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI !clist /METHOD = reml /FIXED = Xprime Iprime !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var11.sav'. dataset name var11. If ($casenum=2) evar11=lag(estimate,1). If ($casenum=2) icc11=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe11.sav". dataset name pe11. If ($casenum=3) ipr11se=Std.Error. If ($casenum=3) ipr11p=sig. If ($casenum=3) ipr11=estimate. If ($casenum=3) xpr11se=lag(Std.Error,1). If ($casenum=3) xpr11p=lag(sig,1). If ($casenum=3) xpr11=lag(estimate,1). If ($casenum=3) int11se=lag(Std.Error,2). If ($casenum=3) int11p=lag(sig,2). If ($casenum=3) int11=lag(estimate,2). execute. compute ncn=$casenum. EXECUTE. select if NCN=3. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Others Model, Main Effects and Interaction Model 11 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev11.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI !clist /METHOD = ml /FIXED = Xprime Iprime !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev11.sav'. dataset name fdev11. compute dev11=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Group Model Model 12 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var12.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe12.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI Xdiv Xmean !clist /METHOD = reml /FIXED = Xmean Xdiv !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var12.sav'. dataset name var12. If ($casenum=2) evar12=lag(estimate,1). If ($casenum=2) icc12=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe12.sav". dataset name pe12. If ($casenum=3) i12se=Std.Error. If ($casenum=3) i12p=sig. If ($casenum=3) i12=estimate. If ($casenum=3) x12se=lag(Std.Error,1). If ($casenum=3) x12p=lag(sig,1). If ($casenum=3) x12=lag(estimate,1). If ($casenum=3) int12se=lag(Std.Error,2). If ($casenum=3) int12p=lag(sig,2). If ($casenum=3) int12=lag(estimate,2). execute. compute ncn=$casenum. EXECUTE. select if NCN=3. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Group Model Model 12 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev12.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI Xdiv Xmean !clist /METHOD = ml /FIXED = Xmean Xdiv !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev12.sav'. dataset name fdev12. compute dev12=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . Comment Constrast Model, Main Effects and Interaction, Model 13 REML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [ ' Covariance Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='var13.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Parameter Estimates'] /DESTINATION FORMAT=SAV OUTFILE='pe13.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI Xcon !clist /METHOD = reml /FIXED = Xcon XconI !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='var13.sav'. dataset name var13. If ($casenum=2) evar13=lag(estimate,1). If ($casenum=2) icc13=estimate. execute. compute ncn=$casenum. EXECUTE. select if NCN=2. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 WaldZ estimate sig lowerbound upperbound ncn command_ subtype_ label_ . execute. get file="pe13.sav". dataset name pe13. If ($casenum=3) i13se=Std.Error. If ($casenum=3) i13p=sig. If ($casenum=3) i13=estimate. If ($casenum=3) x13se=lag(Std.Error,1). If ($casenum=3) x13p=lag(sig,1). If ($casenum=3) x13=lag(estimate,1). If ($casenum=3) int13se=lag(Std.Error,2). If ($casenum=3) int13p=lag(sig,2). If ($casenum=3) int13=lag(estimate,2). execute. compute ncn=$casenum. EXECUTE. select if NCN=3. Compute id=1. MATCH FILES /FILE * /DROP estimate sig DF lowerbound upperbound var1 ncn command_ subtype_ label_ t. execute. Comment Constrast Model, Main Effects and Interaction, Model 13 ML. OMS /SELECT Tables /IF COMMANDS=[' Mixed'] Subtypes= [' Information Criteria'] /DESTINATION FORMAT=SAV OUTFILE='fdev13.sav' /COLUMNS DIMNAMES=[ 'Variables' ]. SET Printback=On. dataset activate gapimdata. MIXED !y WITH !x Xprime I Iprime XconI Xcon !clist /METHOD = ml /FIXED = Xcon XconI !clist /PRINT = SOLUTION TESTCOV /REPEATED=mxmbxr | SUBJECT(!groupid) COVTYPE(CSR). execute. SET Printback=on. OMSEND. get file='fdev13.sav'. dataset name fdev13. compute dev13=var2. compute ncn=$casenum. EXECUTE. select if NCN=1. Compute id=1. MATCH FILES /FILE * /DROP var1 var2 ncn command_ subtype_ label_ . OMS /SELECT Tables /IF COMMANDS=['Regression'] SUBTYPES=['Model Summary'] /DESTINATION FORMAT=SAV OUTFILE='xemx.sav'. dataset activate gapimdata. REGRESSION /MISSING LISTWISE /STATISTICS CHANGE /DEPENDENT !Y /METHOD=ENTER !X Xprime I Iprime /METHOD=ENTER !X Xprime I Iprime !CList. OMSEND. get file="xemx.sav". dataset name xemx. execute. COMPUTE ID=1. Compute ncov=df1. MATCH FILES /FILE * /DROP FChange VAR1 df1 df2 Sig.Fchange Command_ subtype_ label_ RSquareChange. compute ncn=$casenum. EXECUTE. select if NCN=2. execute. Comment Compute Correlations. dataset activate gapimdata. OMS /SELECT Tables /IF COMMANDS=[' Correlations'] Subtypes= [' Correlations'] /DESTINATION FORMAT=SAV OUTFILE='corr.sav' /COLUMNS DIMNAMES=[ 'Variables' 'Labels' 'Statistics']. SET Printback=On. CORRELATIONS /VARIABLES= !x Xprime I Iprime !Y /MISSING=Listwise /Print=Nosig. execute. SET Printback=On. OMSEND. execute. get file="corr.sav". dataset name corr. comment execute. flip. execute. dataset name corr2. IF ($casenum=48) r12=lag(var001,43). IF ($casenum=48) r13=lag(var001,41). IF ($casenum=48) r14=lag(var001,32). IF ($casenum=48) r23=lag(var001,17). IF ($casenum=48) r24=lag(var001,15). IF ($casenum=48) r34=lag(var001,13). compute ncn=$casenum. execute. select if NCN=48. Compute id=1. comment MATCH FILES /FILE * /DROP ncn case_lbl. dataset activate filenox. MATCH FILES /file corr2/file gcount/file file5/file fdevx/file var0/file fdev1/file var1/File=pe1 /file fdev2/file var2/File=pe2/file fdev3/file var3/File=pe3/file fdev3/file var3/File=pe3 /file fdev4/file var4/File=pe4/file fdev5/file var5/File=pe5 /File=pe6/file fdev6/file var6/File=pe7/file fdev7/file var7 /FILE=* BY id. MATCH FILES /File=pe8/file fdev8/file var8/File=pe9/file fdev9/file var9 /File=pe10/file fdev10/file var10/File=pe11/file fdev11/file var11 /File=pe12/file fdev12/file var12/File=pe13/file fdev13/file var13 /File=filek/File=fileq/File=filen/File=freqx/File=filexm/file xemx /FILE=* BY id. dataset close filen. dataset close filek. dataset close fileq. String ess BestMod BestSMod (a50). String pv1 pv2 pv3 pv4 (a1). String row1 row2 row3 rr1 rr2 rr3 rr4 rr5 rr6 rr7 rr8 rr9 rr0 rr10 rr11 rr12 rr13 tto ttu tte tti tty tmp (a109). String rrr1 rrr2 rrr3 rrr4 rrr5 rrr6 rrr7 rrr8 rrr9 rrr10 rrr11 rrr12 rrr13 (a79). String sss1 sss2 sss3 sss4 sss5 sss6 sss7 sss8 sss9 sss10 sss11 sss12 sss13 (a79). String Modname0 Modname1 Modname2 Modname3 Modname4 Modname5 Modname6 Modname7 Modname8 Modname9 Modname10 Modname11 Modname12 Modname13 (a79). String txt0 txt1 txt2 txt3 txt4 txt5 txt6 txt7 tta txt8 txt9 wrn (a1699). Compute xmeann=(xmean-xmin)/(xmax-xmin). Compute zzz =(xmax-xmin)*sqrt(nxn*xmeann*(1-xmeann)/(nxn-1)). Compute xoox=abs(zzz-xsd)/(zzz+xsd). Compute dix=0. If (xoox < .0000001) dix=1 Compute ymeann=(ymean-ymin)/(ymax-ymin). Compute zzz =(ymax-ymin)*sqrt(nxn*ymeann*(1-ymeann)/(nxn-1)). Compute xoox=abs(zzz-ysd)/(zzz+ysd). Compute diy=0. If (xoox < .0000001) diy=1. Compute nx1=trunc((npg_minimum-1)/2). Compute nx2 =npg_minimum - nx1-1. Compute rxr = (nx1*(nx1-1)/2+nx2*(nx2-1)/2-nx1*nx2)/((nx1+nx2)*(nx1+nx2-1)/2). Compute df1=2. Compute df2=1. Compute df3=1. Compute df4=1. Compute df5=1. Compute df6=4. Compute df7=3. Compute df8=3. Compute df9=3. Compute df10=2. Compute df11=2. Compute df12=2. Compute df13=2. Compute SABIC0=emptydev + 3*ln((Ngroups + 2)/24). Compute SABIC1=dev1 + 5*ln((Ngroups + 2)/24). Compute SABIC2=dev2 + 4*ln((Ngroups + 2)/24). Compute SABIC3=dev3 + 4*ln((Ngroups + 2)/24). Compute SABIC4=dev4 + 4*ln((Ngroups + 2)/24). Compute SABIC5=dev5 + 4*ln((Ngroups + 2)/24). Compute SABIC6=dev6 + 7*ln((Ngroups + 2)/24). Compute SABIC7=dev7 + 6*ln((Ngroups + 2)/24). Compute SABIC8=dev8 + 6*ln((Ngroups + 2)/24). Compute SABIC9=dev9 + 6*ln((Ngroups + 2)/24). Compute SABIC10=dev10+ 5*ln((Ngroups + 2)/24). Compute SABIC11=dev11+ 5*ln((Ngroups + 2)/24). Compute SABIC12=dev12+ 5*ln((Ngroups + 2)/24). Compute SABIC13=dev13+ 5*ln((Ngroups + 2)/24). DO REPEAT MulR=MulR1 TO MulR13/evar=evar1 evar2 evar3 evar4 evar5 evar6 evar7 evar8 evar9 evar10 evar11 evar12 evar13. Compute MulR=0. If (evarBestR) MBest=1. If (MulR1>BestR) BestMod="Main Effects Only". If (MulR1>BestR) BestR =MulR1. If (MulR2>BestR) MBest=2. If (MulR2>BestR) BestMod="Actor Only". If (MulR2>BestR) BestR =MulR2. If (MulR3>BestR) MBest=3. If (MulR3>BestR) BestMod="Others Only". If (MulR3>BestR) BestR =MulR3. If (MulR4>BestR) MBest=4. If (MulR4>BestR) BestMod="Main EffectsNorm". If (MulR4>BestR) BestR =MulR4. If (MulR5>BestR) MBest=5. If (MulR5>BestR) BestMod="Main Effects Contrast". If (MulR5>BestR) BestR =MulR5. If (MulR6>BestR) MBest=6. If (MulR6>BestR) BestMod="Complete". If (MulR6>BestR) BestR =MulR6. If (MulR7>BestR) MBest=7. If (MulR7>BestR) BestMod="Person Fit". If (MulR7>BestR) BestR =MulR7. If (MulR8>BestR) MBest=8. If (MulR8>BestR) BestMod="Diversity". If (MulR8>BestR) BestR =MulR8. If (MulR9>BestR) MBest=9. If (MulR9>BestR) BestMod="Interaction Contrast". If (MulR9>BestR) BestR =MulR9. If (MulR10>BestR) MBest=10. If (MulR10>BestR) BestMod="Actor Only (Main Effects and Interaction)". If (MulR10>BestR) BestR =MulR10. If (MulR11>BestR) MBest=11. If (MulR11>BestR) BestMod="Others Only (Main Effects and Interaction)". If (MulR11>BestR) BestR =MulR11. If (MulR12>BestR) MBest=12. If (MulR12>BestR) BestMod="Norm (Main Effects and Interaction)". If (MulR12>BestR) BestR =MulR12. If (MulR13>BestR) MBest=13. If (MulR13>BestR) BestMod="Contrast (Main Effects and Interaction)". If (MulR13>BestR) BestR =MulR13. Compute MBestS=0. Compute BestS =SABIC0. Compute BestSMod="Empty". If (SABIC11) txt1=Concat(Rtrim(txt1)," There are ",Rtrim(tto)," covariates which are all controlled in the analyses."). If (ncov =1) txt1=Concat(Rtrim(txt1)," There is one covariate which is controlled in the analyses."). If (!xmiss=0) txt1=Concat(Rtrim(txt1)," A ",!groupn," containing missing data on ",!xn,", is retained as long as the ",!groupn," has at least three members."). If (!xmiss=1) txt1=Concat(Rtrim(txt1)," Any ",!groupn," which contains missing data on ",!xn," is dropped from the analysis."). Comment Descriptives. Compute tto=Ltrim(String(nxn,F15.0)). Compute tta=Ltrim(String(ngroups,F15.0)). Compute txt2=Concat(" For the analysis, there are ",Rtrim(tto)," persons in ",Rtrim(tta)," ",!groupn,"s."). Compute ttu=Ltrim(String(npg_minimum,F15.0)). Compute tto=Ltrim(String(npg_maximum,F15.0)). If (npg_n=1)txt2=Concat(Rtrim(txt2)," The ",!groupn," size for all ",!groupn,"s is ",Rtrim(ttu),"."). If (npg_n=2)txt2=Concat(Rtrim(txt2)," The two ",!groupn," sizes are ",Rtrim(ttu)," and ",Rtrim(tto),"."). If (npg_n>2)txt2=Concat(Rtrim(txt2)," The ",!groupn," sizes range from ",Rtrim(ttu)," to ",Rtrim(tto),"."). Compute tto=Ltrim(String(xmist_sum,F15.0)). If (!xmiss=1 and xmist_sum>1)txt2=Concat(Rtrim(txt2)," A total of ",Rtrim(tto)," ",!groupn,"s were lost due to missing data on ",!xn,"."). If (!xmiss=1 and xmist_sum=1)txt2=Concat(Rtrim(txt2)," One ",!groupn," was lost due to missing data on ",!xn,"."). If (!Xmiss=1) ymiss = xmist_n - nxn - xmist_sum. If (!Xmiss=0) ymiss = xmist_n - nxn . Compute tto=Ltrim(String(ymiss,F15.0)). If (!xmiss=1 and ymiss>1)txt2=Concat(Rtrim(txt2),"A total of ",Rtrim(tto)," ",!groupn,"s were lost due to missing data on ",!yn,"."). If (!xmiss=1 and ymiss=1)txt2=Concat(Rtrim(txt2),"A total of ",Rtrim(tto)," ",!groupn," was lost due to missing data on ",!yn,"."). If (!xmiss=0 and ymiss>1)txt2=Concat(Rtrim(txt2),"A total of ",Rtrim(tto)," ",!groupn," were lost due to missing data on either ",!yn," or ",!xn,"."). If (!xmiss=0 and ymiss=1)txt2=Concat(Rtrim(txt2),"A total of ",Rtrim(tto)," ",!groupn," was lost due to missing data on either ",!yn,"or ",!xn,"."). Compute tto=Ltrim(String(nhoml,F15.0)). Compute ttu=Ltrim(String(nhomu,F15.0)). Compute txt2=Concat(Rtrim(txt2)," There are ",Rtrim(tto)," all-",!lower_lab," ",!groupn,"s and ",Rtrim(ttu)," all-",!upper_lab," ",!groupn,"s."). Compute ttu=Ltrim(String(nxnt-nxnx,F15.0)). If (nxnt-nxnx>1) txt2=Concat(Rtrim(txt2)," There are ",Rtrim(ttu)," observations missing for ",!xn,"."). If (nxnt-nxnx=1) txt2=Concat(Rtrim(txt2)," There is ",Rtrim(ttu)," one observation missing for ",!xn,"."). Compute ttu=Ltrim(String(nxnt-nxny,F15.0)). If (nxnt-nxny>1) txt2=Concat(Rtrim(txt2)," There are ",Rtrim(ttu)," observations missing for ",!yn,"."). If (nxnt-nxny=1) txt2=Concat(Rtrim(txt2)," There is ",Rtrim(ttu)," observation missing for ",!yn,"."). Compute xhigh=nxn*(xmean+1)/2. Compute xlow=nxn -xhigh. Compute ttu=Ltrim(String(xlow,F15.0)). Compute tto=Ltrim(String(100*xlow/nxn,F15.1)). Compute txt2=Concat(Rtrim(txt2)," There are ",Rtrim(ttu)," (",Rtrim(tto),"%) ",!lower_lab," and"). Compute ttu=Ltrim(String(xhigh,F15.0)). Compute tto=Ltrim(String(100*xhigh/nxn,F15.1)). Compute txt2=Concat(Rtrim(txt2)," ",Rtrim(ttu)," (",Rtrim(tto),"%) ",!upper_lab,"."). Compute txt2=Concat(Rtrim(txt2)," The means and standard deviations of ",!xn," and ",!yn," are presented in Table 1."). Compute tto=Ltrim(String(x6,F15.3)). Compute ttu=Concat(" (p = ",Ltrim(String(x6p,F15.3)),")"). If (x6p< .001) ttu=" (p < .001)". Compute txt3=Concat(Rtrim(txt3)," The effect of the actor's ",!xn," is ",Rtrim(tto),Rtrim(ttu),"." ). Compute tto=Ltrim(String(2*x6,F15.3)). If (x6<0) tto=Ltrim(String(-2*x6,F15.3)). Compute ttu ="lower". if (x6>0) ttu="higher". Compute txt3=Concat(Rtrim(txt3)," ",!upper_lab," score ",Rtrim(tto)," units ",Rtrim(ttu)," on ",!yn," than ",!lower_lab,"."). Compute tto=Ltrim(String(xpr6,F15.3)). Compute ttu=Concat(" (p = ",Ltrim(String(xpr6p,F15.3)),")"). If (xpr6p< .001) ttu=" (p < .001)". Compute txt3=Concat(Rtrim(txt3)," The effect of the ",!xn," of the other ",!groupn," members is ",Rtrim(tto),Rtrim(ttu),"." ). Compute tto=Ltrim(String(2*xpr6,F15.3)). If (xpr6<0) tto=Ltrim(String(-2*xpr6,F15.3)). Compute ttu ="lower". If (xpr6>0) ttu="higher". Compute txt3=Concat(Rtrim(txt3)," An actor, all of whose other ",!groupn," members are ",!upper_lab,", scores ",Rtrim(tto)," units ",Rtrim(ttu)," on ",!yn," than a member all of whose other members are ",!lower_lab,"."). Compute tto=Ltrim(String(i6,F15.3)). If (i6<0) tto=Ltrim(String(-i6,F15.3)). Compute ttu=Concat(" (p = ",Ltrim(String(i6p,F15.3)),")"). If (i6p< .001) ttu=" (p < .001)". Compute txt3=Concat(Rtrim(txt3)," The effect of the actor's similarity to others in the ",!groupn," on ",!xn," is ",Rtrim(tto),Rtrim(ttu),"." ). Compute tto=Ltrim(String(2*i6,F15.3)). Compute ttu ="lower". if (i6>0) ttu="higher". Compute txt3=Concat(Rtrim(txt3)," An actor, who is totally similar to the other ",!groupn," members on ",!xn,", scores ",Rtrim(tto)," units ",Rtrim(ttu)," on ",!yn," than an actor who is totally dissimilar."). Compute tto=Ltrim(String(ipr6,F15.3)). Compute ttu=Concat(" (p = ",Ltrim(String(ipr6p,F15.3)),")"). If (ipr6p< .001) ttu=" (p < .001)". Compute txt3=Concat(Rtrim(txt3)," The effect of the others' similarity in ",!xn," is ",Rtrim(tto),Rtrim(ttu),"." ). Compute tto=Ltrim(String((1-rxr)*ipr6,F15.3)). If (ipr6<0) tto=Ltrim(String(-(1-rxr)*ipr6,F15.3)). Compute ttu ="lower". if (ipr6>0) ttu="higher". Compute tta=Ltrim(String(nx1,F15.0)). Compute tte=Ltrim(String(nx2,F15.0)). Compute tty=Ltrim(String(npg_minimum,F15.0)). Compute txt3=Concat(Rtrim(txt3)," In a ",Rtrim(tty),"-person ",!groupn,", a person with others who are completely homogeneous, either all"). Compute txt3=Concat(Rtrim(txt3)," ",!upper_lab," or ",!lower_lab,", scores ",Rtrim(tto)," units ",Rtrim(ttu)," on ",!yn," than a person with ",Rtrim(tte)," ",!upper_lab). If (nx1=nx2)txt3=Concat(Rtrim(txt3)," and ",Rtrim(tta)," ",!lower_lab," others."). If (nx1<>nx2)txt3=Concat(Rtrim(txt3)," (or ",!lower_lab,") and ",Rtrim(tta)," ",!lower_lab," (or ",!upper_lab,") others."). Compute tto=Ltrim(String(icc6,F15.3)). If (icc6<0) txt4=Concat(" There is a negative correlation between ",!groupn," members' scores, r = ",Rtrim(tto),"."). Compute tto=Ltrim(String(icc6*evar6,F15.3)). If (icc6>=0)txt4=Concat(" The variance due to ",!groupn,"s is equal to ",Rtrim(tto)," and"). Compute tto=Ltrim(String(icc6,F15.3)). If (icc6>=0)txt4=Concat(Rtrim(txt4)," the proportion of variance due to ",!groupn," or intraclass correlation is equal to ",Rtrim(tto),"."). Compute tto=Ltrim(String(MulR6,F15.3)). Compute txt4=Concat(Rtrim(txt4)," The multiple correlation explained by the four GAPIM-I fixed effects is equal to ",Rtrim(tto),"."). Compute tto=Ltrim(String(emptydev-dev6,F15.3)). Compute pchi=1-CDF.CHISQ(emptydev-dev6,4). Compute ttu=Concat(" (p = ",Ltrim(String(pchi,F15.3)),")"). If (pchi< .001) ttu=" (p < .001)". Compute tty="". if (ncov >0) tty= "or covariates". Compute txt4=Concat(Rtrim(txt4)," The chi square test with four degrees of freedom that compares the Complete Model to the Empty Model (a model with no predictors,",Rtrim(tty),") equals ",Rtrim(tto),Rtrim(ttu),"."). If (pchi!alpha)txt4=Concat(Rtrim(txt4)," Because this chi square test is not statistically significant, we do not have evidence that any of the terms of the Complete Model are non-zero."). Compute tto=Ltrim(String(dev1-dev6,F15.3)). Compute pchi=1-CDF.CHISQ(dev1-dev6,4). Compute ttu=Concat(" (p = ",Ltrim(String(pchi,F15.3)),")"). If (pchi< .001) ttu=" (p < .001)". Compute txt4=Concat(Rtrim(txt4)," The chi square test with two degrees of freedom that compares the Complete Model to the Main Effects Model equals ",Rtrim(tto),Rtrim(ttu),"."). If (pchi!alpha)txt4=Concat(Rtrim(txt4)," Because this chi square test is not statistically significant, we do not have evidence that the interaction effects of the Complete Model are non-zero."). Compute tto=Ltrim(String(SABIC6,F15.2)). Compute ttu=Ltrim(String(SABIC0,F15.2)). Compute txt4=Concat(Rtrim(txt4)," The sample size adjusted BIC for the Complete Model is ",Rtrim(tto)," whereas the value for the Empty Model is ",Rtrim(ttu),"."). If (SABIC6SABIC0) txt4=Concat(Rtrim(txt4)," Because the index is larger, the Complete Model is a not a better fitting model than the Empty Model."). Compute tto=Ltrim(String(BestR ,F15.3)). Compute txt5=Concat(" The best fitting model in terms of highest multiple correlation is the ",Rtrim(BestMod)," Model with a multiple R of ",Rtrim(tto),"."). Compute tto=Ltrim(String(BestS ,F15.3)). Compute tty="". If (MBest=MBestS) tty=" also". Compute txt5=Concat(Rtrim(txt5)," (The best fitting model in terms of lowest SABIC is",Rtrim(tty)," the ",Rtrim(BestSMod)," Model with an SABIC of ",Rtrim(tto),".)"). If (MBest=0)txt5=Concat(Rtrim(txt5)," This implies that there are no effects of ",!xn," ",!groupn," composition."). If (MBest=1)txt5=Concat(Rtrim(txt5)," This implies that there an effect of one's own ",!xn," and an effect due other members' ",!xn,"."). If (MBest=2)txt5=Concat(Rtrim(txt5)," This implies that there is only an effect due to one's own ",!xn,"."). If (MBest=3)txt5=Concat(Rtrim(txt5)," This implies that there is only an effect due to the ",!xn," of the other ",!groupn," members."). If (MBest=4)txt5=Concat(Rtrim(txt5)," This implies that there is an effects due to the ",!xn," of the ",!groupn,"."). If (MBest=5)txt5=Concat(Rtrim(txt5)," This implies that there is a contrast effect between one's own ",!xn," and the ",!xn," of the other members' of the ",!groupn). If (MBest=5)txt5=Concat(Rtrim(txt5),". What matters is how different the actor is on",!xn," from the others in the ",!groupn,"."). If (MBest=6)txt5=Concat(Rtrim(txt5)," The Complete Model is discussed above."). If (MBest=7)txt5=Concat(Rtrim(txt5)," This implies that what is important is the match on ",!xn," between the actor and the others in the ",!groupn,"."). If (MBest=8)txt5=Concat(Rtrim(txt5)," This implies that there is an effect due to the ",!groupn,"'s diversity in ",!xn,"."). If (MBest=9)txt5=Concat(Rtrim(txt5)," This implies that there is effect if the actor is different from the others in the ",!groupn," on ",!xn," and the others are similar to each other on ",!xn,"."). If (MBest=10)txt5=Concat(Rtrim(txt5)," This implies is an effect of own ",!xn," and an effect of the actor' similarity to other ",!groupn," members."). If (MBest=11)txt5=Concat(Rtrim(txt5)," This implies is an effect of other ",!groupn," members",!xn," and an effect of the others' similarity to each other."). If (MBest=12)txt5=Concat(Rtrim(txt5)," This implies that the ",!groupn," mean and overall ",!groupn," diversity are important."). If (MBest=13)txt5=Concat(Rtrim(txt5)," This implies that there is a contrast effect both for the effect of ",!xn," and the effects of similarity."). Compute II=0. DO REPEAT df=df1 TO df13/dev= dev1 dev2 dev3 dev4 dev5 dev6 dev7 dev8 dev9 dev10 dev11 dev12 dev13/SABIC=SABIC1 SABIC2 SABIC3 SABIC4 SABIC5 SABIC6 SABIC7 SABIC8 SABIC9 SABIC10 SABIC11 SABIC12 SABIC13. Compute II=II+1. If (II=Mbest) chisqM=emptydev-dev. If (II=MBest) dfM=df. If (II=MBest) SABICM=SABIC. END REPEAT. If (MBest>0) pchi=1-CDF.CHISQ(chisqM,dfM). Compute ttu=Ltrim(String(chisqM,F15.3)). Compute tto=Ltrim(String(dfM,F15.0)). Compute tte=Concat(" (p = ",Ltrim(String(pchi,F15.3)),")"). If (pchi< .001) tte=" (p < .001)". If (Mbest>0) txt5=Concat(Rtrim(txt5)," The chi square test comparing this model to the Empty Model is equal to ",Rtrim(ttu)," with ",Rtrim(tto)," degrees of freedom",Rtrim(tte),"."). Compute ttu=Ltrim(String(SABICM,F15.3)). If (Mbest>0) txt5=Concat(Rtrim(txt5)," The sample size adjusted BIC is equal to ",Rtrim(ttu),"."). Compute tto=Ltrim(String(xmean,F15.3)). Compute ttu=Ltrim(String(xsd,F15.3)). Compute tte=CHAR.LPAD(Rtrim(tto),30-length(!xn)). Compute tmp=CHAR.LPAD(Rtrim(ttu),18). Compute row1=Concat(!xn,Rtrim(tte),Rtrim(tmp)). Compute tto=Ltrim(String(ymean,F15.3)). Compute ttu=Ltrim(String(ysd,F15.3)). Compute tte=CHAR.LPAD(Rtrim(tto),30-length(!yn)). Compute tmp=CHAR.LPAD(Rtrim(ttu),18). Compute row2=Concat(!yn,Rtrim(tte),Rtrim(tmp)). Comment Output 14 Model results. Compute ttu="Empty". Compute tto="0". Compute tto=CHAR.LPAD(Rtrim(tto),38-length(Rtrim(ttu))). Compute tta="0". Compute tta=CHAR.LPAD(Rtrim(tta),13). Compute tti="0". Compute tti=CHAR.LPAD(Rtrim(tti),13). Compute tte="0". Compute tte=CHAR.LPAD(Rtrim(tte),14). Compute tmp=Ltrim(String(SABIC0,F15.3)). Compute tmp=CHAR.LPAD(Rtrim(tmp),15). Compute tty=Ltrim(String(MulR0,F15.3)). Compute tty=CHAR.LPAD(Rtrim(tty),12). Compute rr0=Concat(Rtrim(ttu),Rtrim(tto),Rtrim(tta),Rtrim(tti),Rtrim(tte),Rtrim(tmp),Rtrim(tty)). Compute ttu="Main Effects". Compute pv1=" ". If (x1p1) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). Compute tte=Ltrim(String(xmax,F15.2)). If (xmax>1)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". The maximum value of ",!xn," (",Rtrim(tte),") is greater than one which is not ordinarily allowable. "). If (xmin<-1) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). Compute tte=Ltrim(String(xmin,F15.2)). If (xmin<-1)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". The minimum value of ",!xn," (",Rtrim(tte),") is less than minus one which is not ordinarily allowable. "). If (diy = 1) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). If (diy =1 )wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". The outcome variable ",!yn," is a dichotomy while the analysis assumes it is not. "). If (nxnt-nxnx>0 and !xmiss=0) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). If (nxnt-nxnx>0 and !xmiss=0)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". There are one or more missing observations for the variable ",!xn,", but ",!groupn,"s containing those missing cases are still retained in the analysis."). If (npg_n>1) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). If (npg_n>1)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". Because group sizes vary, the chi square difference tests involving diversity, contrast, and group effects are approximate."). If (abs(xlow/nxn-.5)>.25) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). If (abs(xlow/nxn-.5)>.25)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". The variable ",!xn," is skewed which leads to co-linearity. "). If (r12*r12 >.36 or r13*r13 >.36 or r14*r14 >.36 or r23*r23 >.36 or r24*r24 >.36 or r34*r34 >.36 ) wnum=wnum+1. Compute mco=0. If (r12*r12 >.36 or r13*r13 >.36 or r14*r14 >.36 or r23*r23 >.36 or r24*r24 >.36 or r34*r34 >.36 ) mco=1. Compute tmp=Ltrim(String(wnum,F15.0)). If (mco=1)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". There is multicolinearity:"). Compute tto=Ltrim(String(r12,F15.3)). If (mco=1 and r12*r12 >.36)wrn=Concat(Rtrim(wrn)," The correlation between X and X' is ",Rtrim(tto),"."). Compute tto=Ltrim(String(r13,F15.3)). If (mco=1 and r13*r13 >.36)wrn=Concat(Rtrim(wrn)," The correlation between X and I is ",Rtrim(tto),"."). Compute tto=Ltrim(String(r23,F15.3)). If (mco=1 and r23*r23 >.36)wrn=Concat(Rtrim(wrn)," The correlation between X' and I is ",Rtrim(tto),"."). Compute tto=Ltrim(String(r14,F15.3)). If (mco=1 and r14*r14 >.36)wrn=Concat(Rtrim(wrn)," The correlation between X and I' is ",Rtrim(tto),"."). Compute tto=Ltrim(String(r24,F15.3)). If (mco=1 and r24*r24 >.36)wrn=Concat(Rtrim(wrn)," The correlation between X' and I' is ",Rtrim(tto),"."). Compute tto=Ltrim(String(r34,F15.3)). If (mco=1 and r34*r34 >.36)wrn=Concat(Rtrim(wrn)," The correlation between I and I' is ",Rtrim(tto),"."). If (ngroups<10) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). Compute tto=Ltrim(String(ngroups,F15.0)). If (ngroups<10)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". There are only ",Rtrim(tto)," groups, probably too few for a GAPIM analysis. "). If (ngroups>9 and ngroups<20) wnum=wnum+1. Compute tmp=Ltrim(String(wnum,F15.0)). Compute tto=Ltrim(String(ngroups,F15.0)). If (ngroups>9 and ngroups<20)wrn=Concat(Rtrim(wrn)," ",Rtrim(tmp),". There are only ",Rtrim(tto)," groups, probably a marginal number for a GAPIM analysis. "). If (wnum=0) wrn=" ". If (wnum>1) wrn=Concat("WARNINGS:",Rtrim(wrn)). If (wnum=1) wrn=Concat("WARNING:",Rtrim(wrn)). ALTER TYPE txt0 txt1 txt2 txt3 txt4 txt5 txt6 txt7 txt8 row1 row2 rr0 rr1 rr2 rr3 rr4 rr5 rr6 rr7 rr8 rr9 rr10 rr11 rr12 rr13 rrr1 rrr2 rrr3 rrr4 rrr5 rrr6 rrr7 rrr8 rrr9 rrr10 rrr11 rrr12 rrr13 sss1 sss2 sss3 sss4 sss5 sss6 sss7 sss8 sss9 sss10 sss11 sss12 sss13 wrn (A=AMIN) /print none. Execute. WRITE OUTFILE=!ofile TABLE /wrn/ /"GROUP ACTOR-PARTNER INTERDEPENDENCE MODEL"/ /txt0//txt1 //"RESULTS"/ /"Descriptives"//txt2 / /"Complete Model"//txt3 //txt4//"The Best Fitting Model"//txt5/ / /' Table 1: Descriptive Statistics' //'Variable Mean Standard Deviation' /'--------------------------------------------------------' /row1/row2 / //' Table 2: Model Estimates and Fit' //"Model X X' I I' SABIC R" /'-----------------------------------------------------------------------------------------------------------' /rr0/rr1/rr2/rr3/rr4/rr5 /rr6/rr7/rr8/rr9/"Constraints on both Main Effects and Interactions" /rr10/rr11/rr12/rr13 / //' Table 3: Chi Square Tests' //"Model Tested Chi Square df p" /'------------------------------------------------------' /'EMPTY MODEL AS THE COMPARISON MODEL (significant means a better fitting model)' /rrr1/rrr2/rrr3/rrr4/rrr5 /rrr6/rrr7/rrr8/rrr9/"Constraints on both Main Effects and Interactions" /rrr10/rrr11/rrr12/rrr13 /'MAIN EFFECTS MODEL AS THE COMPARISON MODEL (significant means a worse fitting model)' /sss1/sss2/sss3/sss4 /'COMPLETE MODEL AS THE COMPARISON MODEL (significant means a worse fitting model)' /sss5/sss6/sss7/"Constraints on both Main Effects and Interactions" / sss8/ sss9/ sss10/ sss11/sss12/sss13. dataset activate gapimdata. SAVE OUTFILE='gapimdata.sav' /DROP=Gidn xmist sy filter_$ mxmbxr Xvar /COMPRESSED. execute. dataset close gapimdata. dataset activate mydata. dataset close filenox. dataset close freqn. erase file='freqn.sav'. dataset close freqx. erase file='freqx.sav'. dataset close file5. erase file='file5.sav'. dataset close gcount. erase file='gcount.sav'. dataset close corr2. dataset close corr. erase file='corr.sav'. dataset close filet. erase file="filet.sav". dataset close filex. erase file="filex.sav". dataset close filey. erase file="filey.sav". dataset close filexm. erase file="filexm.sav". dataset close var0. dataset close var1. dataset close var2. dataset close var3. dataset close var4. dataset close var5. dataset close var6. dataset close var7. dataset close var8. dataset close var9. dataset close var10. dataset close var11. dataset close var12. dataset close var13. dataset close fdevx. dataset close fdev1. dataset close fdev2. dataset close fdev3. dataset close fdev4. dataset close fdev5. dataset close fdev6. dataset close fdev7. dataset close fdev8. dataset close fdev9. dataset close fdev10. dataset close fdev11. dataset close fdev12. dataset close fdev13. dataset close pe1. dataset close pe2. dataset close pe3. dataset close pe4. dataset close pe5. dataset close pe6. dataset close pe7. dataset close pe8. dataset close pe9. dataset close pe10. dataset close pe11. dataset close pe12. dataset close pe13. dataset close xemx. erase file='var0.sav'. erase file='var1.sav'. erase file='var2.sav'. erase file='var3.sav'. erase file='var4.sav'. erase file='var5.sav'. erase file='var6.sav'. erase file='var7.sav'. erase file='var8.sav'. erase file='var9.sav'. erase file='var10.sav'. erase file='var11.sav'. erase file='var12.sav'. erase file='var13.sav'. erase file='fdevx.sav'. erase file='fdev1.sav'. erase file='fdev2.sav'. erase file='fdev3.sav'. erase file='fdev4.sav'. erase file='fdev5.sav'. erase file='fdev6.sav'. erase file='fdev7.sav'. erase file='fdev8.sav'. erase file='fdev9.sav'. erase file='fdev10.sav'. erase file='fdev11.sav'. erase file='fdev12.sav'. erase file='fdev13.sav'. erase file='pe1.sav'. erase file='pe2.sav'. erase file='pe3.sav'. erase file='pe4.sav'. erase file='pe5.sav'. erase file='pe6.sav'. erase file='pe7.sav'. erase file='pe8.sav'. erase file='pe9.sav'. erase file='pe10.sav'. erase file='pe11.sav'. erase file='pe12.sav'. erase file='pe13.sav'. erase file='xemx.sav'. execute. !ENDDEFINE.