afni/stats/

[giscus[bot]]

afni/stats/

https://liaopan.github.io/afni/stats/

[LL2421]

我怎么导入make_Timings.sh脚本，输出的txt文件是空的，1D文件是4个星号。

[LL2421]

此的脚本可以提取：#!/bin/bash

检查并创建被试列表

if [ ! -f subjList.txt ]; then
ls -d sub-* > subjList.txt 2>/dev/null
fi

检查是否找到被试

if [ ! -s subjList.txt ]; then
echo "错误: 没有找到任何被试目录（sub-*）"
exit 1
fi

echo "总共有 $(cat subjList.txt | wc -l) 个被试"

循环处理每个被试

for subj in cat subjList.txt; do
echo "正在处理: $subj"

# 检查被试目录是否存在
if [ ! -d "$subj/func" ]; then
    echo "  警告: $subj/func 目录不存在"
    continue
fi

# 进入func目录
cd "$subj/func" || { echo "无法进入 $subj/func 目录"; cd ../..; continue; }

# 处理run1
run1_file="${subj}_task-flanker_run-1_events.tsv"
if [ -f "$run1_file" ]; then
    echo "  处理 $run1_file"
    # 提取incongruent_correct事件（跳过标题行）
    tail -n +2 "$run1_file" | awk -F'\t' '$3=="incongruent_correct" {print $1, $2, 1}' > incongruent_run1.txt
    # 提取congruent_correct事件
    tail -n +2 "$run1_file" | awk -F'\t' '$3=="congruent_correct" {print $1, $2, 1}' > congruent_run1.txt
    
    echo "  incongruent_run1.txt: $(wc -l < incongruent_run1.txt) 行"
    echo "  congruent_run1.txt: $(wc -l < congruent_run1.txt) 行"
else
    echo "  错误: 找不到 $run1_file"
    # 返回上级目录
    cd ../..
    continue
fi

# 处理run2
run2_file="${subj}_task-flanker_run-2_events.tsv"
if [ -f "$run2_file" ]; then
    echo "  处理 $run2_file"
    tail -n +2 "$run2_file" | awk -F'\t' '$3=="incongruent_correct" {print $1, $2, 1}' > incongruent_run2.txt
    tail -n +2 "$run2_file" | awk -F'\t' '$3=="congruent_correct" {print $1, $2, 1}' > congruent_run2.txt
    
    echo "  incongruent_run2.txt: $(wc -l < incongruent_run2.txt) 行"
    echo "  congruent_run2.txt: $(wc -l < congruent_run2.txt) 行"
else
    echo "  警告: 找不到 $run2_file"
fi

# 转换为AFNI格式
echo "  转换为AFNI格式..."

# 检查是否有congruent的txt文件
if [ -f "congruent_run1.txt" ] || [ -f "congruent_run2.txt" ]; then
    echo "  转换congruent文件"
    timing_tool.py -fsl_timing_files congruent_run*.txt -write_timing congruent.1D
    echo "  congruent.1D 内容:"
    cat congruent.1D
else
    echo "  注意: 没有congruent的txt文件"
    echo "****" > congruent.1D
fi

# 检查是否有incongruent的txt文件
if [ -f "incongruent_run1.txt" ] || [ -f "incongruent_run2.txt" ]; then
    echo "  转换incongruent文件"
    timing_tool.py -fsl_timing_files incongruent_run*.txt -write_timing incongruent.1D
    echo "  incongruent.1D 内容:"
    cat incongruent.1D
else
    echo "  注意: 没有incongruent的txt文件"
    echo "****" > incongruent.1D
fi

# 返回到原始目录
cd ../..
echo ""

done

echo "处理完成！"

显示统计信息

echo ""
echo "====== 统计信息 ======"
for subj in cat subjList.txt; do
if [ -d "$subj/func" ]; then
echo "$subj:"
for cond in congruent incongruent; do
if [ -f "$subj/func/$cond.1D" ]; then
echo " $cond.1D: $(grep -v '^$' "$subj/func/$cond.1D" | wc -l) 行"
else
echo " $cond.1D: 不存在"
fi
done
fi
done