Hello, i am trying to do a dual regression on HPC, but it’s not working. Each time it seems to work but I don’t have any output.
Here is my code:
module load GCC/8.3.0
module load OpenMPI/3.1.4
module load FSL/6.0.3-Python-3.7.4
\# The melodic output
MELODIC_IC="/home/users/l/lemarec4/data_ICA/GroupICA_T0_T1.ica/melodic_IC.nii.gz"
\# Output directory for dual regression results
RESULTS_DIR="/home/users/l/lemarec4/Results_dualreg"
\# Must be the SAME list (same order) used for MELODIC
INPUT_LIST="/home/users/l/lemarec4/data_ICA/inputs_4_ica.txt"
\# Perform thresholded dual regression to obtain unbiased timeseries for connectomics analyses (e.g., with FSLnets)
THR_FLAG="--thr"
mkdir -p "${RESULTS_DIR}"
\# ---------------------------
\# 3) Run Dual Regression (Stages 1 & 2 only)
\# ---------------------------
\# Arguments: <groupICs> <des_norm> <design.mat> <design.con> <n_perm> <outdir> <subjects...>
\# Here: "0 -1 0" means: no group-level stage 3 here (stages 1 & 2 only)
echo "Running dual_regression on inputs from: ${INPUT_LIST}"
srun dual_regression "${MELODIC_IC}" 0 -1 0 ${THR_FLAG} "${RESULTS_DIR}" \`cat "${INPUT_LIST}"\`
echo "Done. Main outputs:"
echo " - ${RESULTS_DIR}/dr_stage1/ (time courses)"
echo " - ${RESULTS_DIR}/dr_stage2_icXXXX.nii.gz (subject spatial maps per IC)"
Can you help me ?
Thank you so much !