Hi everyone. This one is probably for Gang,
I have 2 quick questions about 3dLME syntax.
Currently I have 2 within-subjects fixed variables - cs (csm, csp, cst) and scan (scan1, scan2) - and I’d like to investigate the difference between each cs within the same scan and between the same cs between the two scans. So I’d like t-test between csm & csp, csm & cst, and csp & cst during scan1 and then during scan2 while also getting t-tests between the csm in scan 1 and the csm in scan2 (and also the csp & cst). I have set up my 3dLME call as follows with gltCode 1-6 for the companions within the same scans the gltCode 7-9 for the comparisons between the same cs but different scans. While I think the first 1-6 are right looking at the help page for 3dLME, the contrasts 7-9 threw me off because there were no examples of 2 within-subjects variables so I just want to ensure this will get me the contrasts I am looking for.
3dLME -prefix LME_test1 -jobs 24 \ -model "cs*scan+sex+stais+stait" \ -qvars "sex,stais,stait" \ -ranEff '~1+CS" \ -SS_type 3 \ -num_glt 9 \ -gltLabel 1 'csm-csp1' -gltCode 1 'cs : 1*csm -1*csp scan : 1*scan1' \ -gltLabel 2 'csp-cst1' -gltCode 2 'cs : 1*csp -1*cst scan : 1*scan1' \ -gltLabel 3 'csm-cst1' -gltCode 3 'cs : 1*csm -1*cst scan : 1*scan1' \ -gltLabel 4 'csm-csp2' -gltCode 4 'cs : 1*csm -1*csp scan : 1*scan2' \ -gltLabel 5 'csp-cst2' -gltCode 5 'cs : 1*csp -1*cst scan : 1*scan2' \ -gltLabel 6 'csm-cst2' -gltCode 6 'cs : 1*csm -1*cst scan : 1*scan2' \ -gltLabel 7 'csm1-csm2' -gltCode 7 'cs : 1*csm scan : 1*scan1 -1*scan2' \ -gltLabel 8 'csp1-csp2' -gltCode 8 'cs : 1*csp scan : 1*scan1 -1*scan2' \ -gltLabel 9 'cst1-csmt' -gltCode 9 'cs : 1*cst scan : 1*scan1 -1*scan2' \ -num_glf 2 \ -glfLabel 1 'csm-csp1' -glfCode 1 'scan : 1*scan1 cs : 1*csm -1*csp & 1*csm -1*cst' \ -glfLabel 2 'csm-csp2' -glfCode 2 'scan : 1*scan2 cs : 1*csm -1*csp & 1*csm -1*cst' \ -dataTable @table.txt
My last question - I have read in another post you can’t have something like a ‘sex’ variable in your model as I have above because it is not quantitative. is this true? The ‘condition’ codes in the examples are not quantitative either but then these are not covariates. Since I am interested in controlling for sex in my model is there a way to include this kind of information or no?
EDIT: I thought of one more question in the meantime. For the input files should these be the appropriate sub-briks of the stats dataset (like in the 3dANOVAs) or something else like the REML datasets?