History of AFNI updates  

|
October 03, 2015 04:28PM
Hi,

I need some help fixing my script for make_random_timing.py. Based on the help file for make_random_timing.py, the how to #03 example, as well as the the first example in this thread I am trying to come up with a script that includes 100 iterations and outputs the normalized standard deviation for each trial. I have made many efforts but my scripting skills are not good enough to set up the extraction of the LC sums correctly.

Thanks,
Maria

set iterations = 100
set seed = 1234567
set outdir = stim_results
set LCfile = NSD_sums # file to store norm. std. dev. sums in # set pattern = LC # search pattern for LC[0], say
set pattern = 'norm. std.' # search pattern for normalized stdev vals

foreach iter (`count -digits 3 1 $iterations`)

# create empty LC file
echo -n "" > $LCfile

echo -n "iteration: 000"

make_random_timing.py -num_runs 2 -run_time 240 -num_stim 2 -num_reps 15 -show_timing_stats -stim_dur 2 -max_consec 3 -min_rest 2.0 -save_3dd_cmd sample.3dd.command.${iter} -stim_labels first second -t_gran 7.0 -make_3dd_contrasts -verb 4 -seed ${seed} -prefix stim_times.${iter} \
>& out.mrt.$iter
tcsh cmd.3dd.$iter >& out.3dD.$iter

# save the sum of the 2 LC values
set nums = ( `awk -F= '/'"$pattern"'/ {print $2}' out.3dD.${iter}` )

# make a quick ccalc command
set sstr = $nums[1]
foreach num ( $nums[2-1] )
set sstr = "$sstr + $num"
end
set num_sum = `ccalc -expr "$sstr"`

echo -n "$num_sum = $sstr : " >> $LCfile
echo "iteration $iter, seed $seed" >> $LCfile

echo -n "$iter"
end
Subject Author Posted

make_random_timing.py

Kristina November 08, 2012 03:35PM

Re: make_random_timing.py

rick reynolds November 08, 2012 10:26PM

Re: make_random_timing.py

Kristina November 13, 2012 10:18AM

Re: make_random_timing.py

rick reynolds November 14, 2012 02:35PM

Re: make_random_timing.py

Gabe Castillo November 09, 2012 02:28PM

Re: make_random_timing.py

rick reynolds November 10, 2012 02:21AM

Re: make_random_timing.py

maria October 03, 2015 04:28PM

Re: make_random_timing.py

rick reynolds October 03, 2015 09:10PM

Re: make_random_timing.py

maria October 03, 2015 09:53PM

Re: make_random_timing.py

rick reynolds October 03, 2015 11:41PM

Re: make_random_timing.py

maria October 04, 2015 02:21AM

Re: make_random_timing.py

rick reynolds October 04, 2015 10:35PM