Dear AFNIsts:
Recently I began trying to run 3dMEMA, but I didn't get very far before encountering a problem that I couldn't figure out how to solve. I tried running 3dMEMA using a command where script.F-R.R contains the input arguments and tempout.F-R.R contains the output and the full error message. Here's the command:
Rscript script.F-R.R | & tee tempout.F-R.R
The contents of the input and output files are pasted below. I have the latest version of AFNI with R installed, running on Linux CentOS 5.2.
Any help would be greatly appreciated.
Thanks in advance,
Will
script.F-R.R:
source("/usr/local/afni/3dMEMA.R")
mema.F-R
3
1
3
F-R
22
4014
3818
4007
3804
3787
3734
3696
3810
3418
3649
3602
3800
3903
4072
3343
4081
3959
3989
3846
4009
4095
4120
F
4014.F.bl5+tlrc.BRIK
4014.F.t.bl5+tlrc.BRIK
3818.F.bl5+tlrc.BRIK
3818.F.t.bl5+tlrc.BRIK
4007.F.bl5+tlrc.BRIK
4007.F.t.bl5+tlrc.BRIK
3804.F.bl5+tlrc.BRIK
3804.F.t.bl5+tlrc.BRIK
3787.F.bl5+tlrc.BRIK
3787.F.t.bl5+tlrc.BRIK
3734.F.bl5+tlrc.BRIK
3734.F.t.bl5+tlrc.BRIK
3696.F.bl5+tlrc.BRIK
3696.F.t.bl5+tlrc.BRIK
3810.F.bl5+tlrc.BRIK
3810.F.t.bl5+tlrc.BRIK
3418.F.bl5+tlrc.BRIK
3418.F.t.bl5+tlrc.BRIK
3649.F.bl5+tlrc.BRIK
3649.F.t.bl5+tlrc.BRIK
3602.F.bl5+tlrc.BRIK
3602.F.t.bl5+tlrc.BRIK
3800.F.bl5+tlrc.BRIK
3800.F.t.bl5+tlrc.BRIK
3903.F.bl5+tlrc.BRIK
3903.F.t.bl5+tlrc.BRIK
4072.F.bl5+tlrc.BRIK
4072.F.t.bl5+tlrc.BRIK
3343.F.bl5+tlrc.BRIK
3343.F.t.bl5+tlrc.BRIK
4081.F.bl5+tlrc.BRIK
4081.F.t.bl5+tlrc.BRIK
3959.F.bl5+tlrc.BRIK
3959.F.t.bl5+tlrc.BRIK
3989.F.bl5+tlrc.BRIK
3989.F.t.bl5+tlrc.BRIK
3846.F.bl5+tlrc.BRIK
3846.F.t.bl5+tlrc.BRIK
4009.F.bl5+tlrc.BRIK
4009.F.t.bl5+tlrc.BRIK
4095.F.bl5+tlrc.BRIK
4095.F.t.bl5+tlrc.BRIK
4120.F.bl5+tlrc.BRIK
4120.F.t.bl5+tlrc.BRIK
R
4014.R.bl5+tlrc.BRIK
4014.R.t.bl5+tlrc.BRIK
3818.R.bl5+tlrc.BRIK
3818.R.t.bl5+tlrc.BRIK
4007.R.bl5+tlrc.BRIK
4007.R.t.bl5+tlrc.BRIK
3804.R.bl5+tlrc.BRIK
3804.R.t.bl5+tlrc.BRIK
3787.R.bl5+tlrc.BRIK
3787.R.t.bl5+tlrc.BRIK
3734.R.bl5+tlrc.BRIK
3734.R.t.bl5+tlrc.BRIK
3696.R.bl5+tlrc.BRIK
3696.R.t.bl5+tlrc.BRIK
3810.R.bl5+tlrc.BRIK
3810.R.t.bl5+tlrc.BRIK
3418.R.bl5+tlrc.BRIK
3418.R.t.bl5+tlrc.BRIK
3649.R.bl5+tlrc.BRIK
3649.R.t.bl5+tlrc.BRIK
3602.R.bl5+tlrc.BRIK
3602.R.t.bl5+tlrc.BRIK
3800.R.bl5+tlrc.BRIK
3800.R.t.bl5+tlrc.BRIK
3903.R.bl5+tlrc.BRIK
3903.R.t.bl5+tlrc.BRIK
4072.R.bl5+tlrc.BRIK
4072.R.t.bl5+tlrc.BRIK
3343.R.bl5+tlrc.BRIK
3343.R.t.bl5+tlrc.BRIK
4081.R.bl5+tlrc.BRIK
4081.R.t.bl5+tlrc.BRIK
3959.R.bl5+tlrc.BRIK
3959.R.t.bl5+tlrc.BRIK
3989.R.bl5+tlrc.BRIK
3989.R.t.bl5+tlrc.BRIK
3846.R.bl5+tlrc.BRIK
3846.R.t.bl5+tlrc.BRIK
4009.R.bl5+tlrc.BRIK
4009.R.t.bl5+tlrc.BRIK
4095.R.bl5+tlrc.BRIK
4095.R.t.bl5+tlrc.BRIK
4120.R.bl5+tlrc.BRIK
4120.R.t.bl5+tlrc.BRIK
12
1
1
~/data/masks/EZML.mask.333+tlrc.BRIK
0
0
tempout.F-R.R:
[1] "#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
[1] " ================== Welcome to 3dMetaAna.R ================== "
[1] "AFNI Mixed-Effects Meta-Analysis Modeling Package!"
[1] "#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
[1] "Version 0.1.1, Sept. 9, 2009"
[1] "Author: Gang Chen (gangchen@mail.nih.gov)"
[1] "Website - [
afni.nimh.nih.gov];
[1] "SSCC/NIMH, National Institutes of Health, Bethesda MD 20892"
[1] "#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++"
[1] "################################################################"
[1] "Please consider citing the following if this program is useful for you:"
Gang Chen et al., hopefully something coming soon...
[
afni.nimh.nih.gov]
[1] "################################################################"
[1] "Use CNTL-C on Unix or ESC on GUI version of R to stop at any moment."
[1] "On a multi-processor machine, parallel computing will speed up the program significantly."
[1] "Choose 1 for a single-processor computer."
[1] "-----------------"
[1] "The following types of group analysis are currently available:"
[1] "1: one condition with one group;"
[1] "2: one condition across 2 groups with homoskedasticity (same variability);"
[1] "3: two conditions with one group;"
[1] "4: one condition across 2 groups with heteroskedasticity (different variability)."
[1] "Since the contrast between the 2 conditions will be the 1st minus the 2nd, choose"
[1] "an appropriate order between the 2 conditions to get the desirable contrast."
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
[1] "-----------------"
Error in if (as.integer(size) == size) { :
missing value where TRUE/FALSE needed
Calls: source -> eval.with.vis -> eval.with.vis -> lapply -> FUN
Execution halted
