I just modified 3dMEMA so that 3dMEMA could take 1D input files as follows (notice that the output file name has to have an affix of 1D):
3dMEMA -prefix bothLOC.hrf2.Acc.1D -jobs 1 -set A MRI_Subj_1 3dMEMA_S1_coef.1D 3dMEMA_S1_T.1D \
MRI_Subj_2 3dMEMA_S2_coef.1D 3dMEMA_S2_T.1D \
MRI_Subj_3 3dMEMA_S3_coef.1D 3dMEMA_S3_T.1D \
MRI_Subj_4 3dMEMA_S4_coef.1D 3dMEMA_S4_T.1D \
MRI_Subj_5 3dMEMA_S5_coef.1D 3dMEMA_S5_T.1D \
MRI_Subj_6 3dMEMA_S6_coef.1D 3dMEMA_S6_T.1D \
MRI_Subj_7 3dMEMA_S7_coef.1D 3dMEMA_S7_T.1D \
MRI_Subj_8 3dMEMA_S8_coef.1D 3dMEMA_S8_T.1D \
MRI_Subj_9 3dMEMA_S9_coef.1D 3dMEMA_S9_T.1D \
MRI_Subj_10 3dMEMA_S10_coef.1D 3dMEMA_S10_T.1D \
-cio
Alternatively you don't have to put those numbers in 1D files as input, and instead do something like this:
3dMEMA -prefix bothLOC.hrf2.Acc.1D -jobs 1 -set A \
MRI_Subj_1 'R:-2.347991' 'R:-0.946193' \
MRI_Subj_2 'R:-1.729988' 'R:-0.823841' \
MRI_Subj_3 'R:6.018984' 'R:0.619669' \
MRI_Subj_4 'R:-5.100689' 'R:-2.583102' \
MRI_Subj_5 'R:-7.328274' 'R:-0.515742' \
MRI_Subj_6 'R:-6.633659' 'R:-3.660589' \
MRI_Subj_7 'R:-2.653562' 'R:-1.328515' \
MRI_Subj_8 'R:-4.423012' 'R:-1.050487' \
MRI_Subj_9 'R:-14.062942' 'R:-1.254944' \
MRI_Subj_10 'R:-7.640000' 'R:-4.284862' \
-cio
In the output file bothLOC.hrf2.Acc.1D there are four values:
-4.580125
-5.357312
1.9472
9.801461
they are beta, t-stat, tau^2 and Q.
I've submitted the modified program for distribution (probably available tomorrow for you through update). Let me know if it works for you.
Gang
Edited 1 time(s). Last edit at 11/22/2013 02:31PM by Gang.