AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
July 10, 2015 07:30AM
Hi,

I am having problem setting up a 3DLME with between-subject variables. I know my R-package is working because I can run 3DLME with another data set I have.

Here is my code:

3dLME -prefix 150709_78s \
-mask /Masks/amygdala.mask+tlrc \
-jobs 2 \
-model "Gene*Age*Cond" \
-qVars Age \
-ranEff '~1' -SS_type 3 \
-dataTable \
Subj Gene Age Cond InputFile \
s1 WT 22.6 F imput/stats.s1+tlrc'[F#0_Coef]' \
s1 WT 22.6 C imput/stats.s1+tlrc'[C#0_Coef]' \
s2 WT 23.7 F imput/stats.s2+tlrc'[F#0_Coef]' \
s2 WT 23.7 C imput/stats.s2+tlrc'[C#0_Coef]' \
s3 MT 21 F imput/stats.s3+tlrc'[F#0_Coef]' \
s3 MT 21 C imput/stats.s3+tlrc'[C#0_Coef]' \
s4 WT 23.9 F imput/stats.s4+tlrc'[F#0_Coef]' \
s4 WT 23.9 C imput/stats.s4+tlrc'[C#0_Coef]' \
s5 MT 23.1 F imput/stats.s5+tlrc'[F#0_Coef]' \
s5 MT 23.1 C imput/stats.s5+tlrc'[C#0_Coef]' \

.... and so on. I have 70+ subjects in my data table.

When I run the code, the contingency table looks correct. But after the end of data structure info, it gives this error.


Reading input files now...

If the program hangs here for more than, for example, half an hour,
kill the process because the model specification or the GLT coding
is likely inappropriate.

~~~~~~~~~~~~~~~~~~~ Model test failed ~~~~~~~~~~~~~~~~~~~
Possible reasons:

0) Make sure that R packages nlme and phia have been installed. See the 3dLME
help documentation for more details.

1) Inappropriate model specification with options -model, or -qVars.

2) In correct specifications in general linear test coding with -gltCode.

3) Mistakes in data table. Check the data structure shown above, and verify
whether there are any inconsistencies.

4) Inconsistent variable names which are case sensitive. For example, factor
named Group in model specification and then listed as group in the table hader
would cause grief for 3dLME.

** Error:
Quitting due to model test failure...

Can you see any problem with my code?

Also even if I simplify the model and run, say, "Cond" for the model, the script still fails. I almost suspect that there is some problem with the input files. And if so, how can I go about verifying that? I have tried opening my input files with Afni GUI and 3dinfo. They looked fine when I did that.

Thank you so much
Subject Author Posted

3DLME model failure

ace327 July 10, 2015 07:30AM

Re: 3DLME model failure

gang July 10, 2015 10:05AM

Re: 3DLME model failure

ace327 July 10, 2015 10:46AM