RestingStateStats.bat

92 views
Skip to first unread message

Lin Cai

unread,
Jul 26, 2022, 2:35:56 AM7/26/22
to HCP-Users
Dear HCP-experts,

When I run the RestingStateStats.bat in the HCPpipelines folder,
I got an error message as follows.
*******************************************************************************************
Error using load
Unable to read file '/mnt/data2/derivatives/HCPpipeline/sub-K1908061200/MNINonLinear/Results/PA_run-01_epi/PA_run-01_epi_hp0.ica/.fix'. No such file or
directory.

Error in RestingStateStats (line 115)
Inoise=load(noiselist);
*******************************************************************************************

I modified the RestingStateStats.bat file according to my own data structure, e.g., I have 4 runs for resting-state fMRI. 
To my understanding, for each run, this batch file will generate one RestingStateStats folder, right?
Besides, I have successfully obtained the MSMAll files according to the sequence of (1)PreFreeSurfer(2)FreeSurfer(3)PostFreeSurfer(4)GenericfMRIVolume(5)GenericfMRISurface(6)IcaFixProcessing(7)PostFix(8)MSMAll(9)
DedriftAndResample.
Could you please help me find which steps possibly result in the current error?

The following settings are for RestingStateStats.bat.
*******************************************************************************************
#Parameter setting
fMRINames="PA_run-01_epi AP_run-01_epi PA_run-02_epi AP_run-02_epi"
OrigHighPass="0" #Specified in Sigma
Caret7_Command="wb_command"
#RegName="MSMRSNOrig3_d26_DR_DeDrift"
#RegName="MSMAll_2_d41_WRN_DeDrift"
RegName="NONE"

LowResMesh="32"
FinalfMRIResolution="2"
BrainOrdinatesResolution="2"
SmoothingFWHM="2"
OutputProcSTRING="_hp0_clean"
dlabelFile="/mnt/data2/cai_practice/Covid_Project/code/Q1-Q6_RelatedValidation210.CorticalAreas_dil_Final_Final_Areas_Group_Colors.32k_fs_LR.dlabel.nii"
MatlabRunMode="1"
BCMode="NONE" #One of REVERT (revert bias field correction), NONE (don't change biasfield correction), CORRECT (revert original bias field correction and apply new one, requires ??? to be present)
OutSTRING="stats"
WM="$HCPPIPEDIR/global/config/FreeSurferWMRegLut.txt"
CSF="$HCPPIPEDIR/global/config/FreeSurferCSFRegLut.txt" 
*******************************************************************************************

Thanks in advance,
Lin

Lin Cai

unread,
Jul 27, 2022, 4:24:56 AM7/27/22
to hcp-...@humanconnectome.org
Dear HCP-experts,

I tried to run RestingStateStats.bat again by setting fMRINames="run12_PA_AP" which was the name of a concatenated multi-run FIX.  
This time, I came across another problem, despite kind of solving the previous problem. The error message is as follows.

**************************************************************************************************************************
RestingStateStats - About to reshape
RestingStateStats - About to remove fakeNIFTI fileRestingStateStats - Finished fslmaths filtering of cifti input
Error using load
Unable to read file '/mnt/data2/cai_practice/Covid_Project/derivatives/HCPpipeline/sub-K1908061200/MNINonLinear/Results/run12_PA_AP/Movement_Regressors.txt'. No such file or
directory.

Error in RestingStateStats (line 171)
  confounds=load([motionparameters '.txt']);
**************************************************************************************************************************

I noticed that I need the Movement_Regressors.txt file after I checked Line 919 of the script. 
Line 919 in RestingStateStats.sh: motionparameters="${ResultsFolder}/Movement_Regressors" #No .txt

BUT, in my corresponding folder(/MNINonLinear/Results/run12_PA_AP/),
I just found two following related txt files (e.g., Movement_Regressors_demean.txt, Movement_Regressors_hp0_clean.txt),
even though Movement_Regressors.txt indeed exists for each run(e.g., /MNINonLinear/Results/AP_run-01_epi/).

So I am curious about whether I can modify the script to use Movement_Regressors_demean.txt instead?
Or I did something wrong  when I preprocessed the functional data, resulting in incompleteness of output files. 

Thank you very much in advance.

Lin



--
You received this message because you are subscribed to the Google Groups "HCP-Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hcp-users+...@humanconnectome.org.
To view this discussion on the web visit https://groups.google.com/a/humanconnectome.org/d/msgid/hcp-users/7cc9b5f7-1efb-4173-92fa-e64f326af4d8n%40humanconnectome.org.


--
Lin Cai

Project Researcher
University of Tokyo, Center for Evolutionary Cognitive Sciences
Tokyo, Japan

Harms, Michael

unread,
Jul 27, 2022, 11:15:21 AM7/27/22
to hcp-...@humanconnectome.org

 

Hi,

I’m not sure if RestingStateStats has been updated to work with MR-FIX.  Matt is on vacation, so you may have to repost in a couple weeks for a more definitive response.

 

Cheers,

-MH

 

-- 

Michael Harms, Ph.D.

-----------------------------------------------------------

Associate Professor of Psychiatry

Washington University School of Medicine

Department of Psychiatry, Box 8134

660 South Euclid Ave.                        Tel: 314-747-6173

St. Louis, MO  63110                          Email: mha...@wustl.edu

 

From: Lin Cai <bertr...@gmail.com>
Reply-To: "hcp-...@humanconnectome.org" <hcp-...@humanconnectome.org>
Date: Wednesday, July 27, 2022 at 3:25 AM
To: "hcp-...@humanconnectome.org" <hcp-...@humanconnectome.org>
Subject: Re: [hcp-users] RestingStateStats.bat

 

* External Email - Caution *

 


The materials in this message are private and may contain Protected Healthcare Information or other information of a sensitive nature. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this email in error, please immediately notify the sender via telephone or return mail.

Glasser, Matt

unread,
Aug 6, 2022, 10:50:20 PM8/6/22
to hcp-...@humanconnectome.org

The released version of RestingStateStats does not work with MR+FIX.  RestingStateStats overall is in need of a bit of a rewrite.  It needs to partition variance into:

 

*Original

*Detrend Structured Noise

*sICA+FIX Structured Noise

*tICA Structured Noise

*BOLD Signal

*Unstructured Noise

 

Things like MGTR/WM/CSF/Motion/Physio regressors if still provided should be shown as what they do to the data above and beyond the above recommended denoising steps.  The effect of Wishart Filtering might also be useful to show.  If there are specific measures you need now from RestingStateStats, it might be best to have a look at the code and compute them yourself.

 

Matt.

Reply all
Reply to author
Forward
0 new messages